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CHAPTER  I 


INTRODUCTION 


1.1  Summary  of  Thesis 

The  green  line  (5577A)  is  a  bright,  persistent  component  of  the  visible 
airglow.  It  is  produced  by  an  electric  quadrupole  transition  from  the  meta¬ 
stable  second  excited  state  (^Sq)  to  the  first  excited  state  (ID2)  of  atomic 
oxygen.  These  two  excited  states  all  lie  in  the  same  electron  shell  of  the 
atom  and  have  the  same  electron  configuration  as  the  ground  state  of 
Is22s22p'^,  which  is  the  ^P2,i,o- 

This  emission  is  present  in  both  the  daytime  and  night  airglow  and 
in  the  aurora,  and  despite  a  long  history  of  study  it  is  still  not  fully 
understood.  The  emission  in  the  dayglow  and  the  nightglow  is  relatively 
homogeneous  spatially  and  global  in  coverage.  In  the  aurora,  the  emission 
is  much  brighter  than  the  airglow,  highly  structured  and  very  localized 
being  restricted  to  higher  latitudes. 

The  structure  of  the  5577A  emission  with  altitude  and  the  chemistry 
responsible  for  the  production  of  the  emission  are  complex.  The  vertical 
structure  for  the  emission  has  two  distinct  layers  in  the  airglow  each  with 
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its  own  set  of  production  and  loss  mechanisms.  The  chemistry  for  either  of 
these  layers  is  not  completely  known.  The  auroral  emission  is  not 
understood  either  since  it  overlaps  the  upper  and  lower  layer  altitudes  and 
it  tends  to  contain  some  parts  of  the  chemistry  of  both  layers  as  sources  and 
losses. 

This  thesis  addresses  the  upper  layer  emission  in  the  dayglow.  The 
chemistry  that  has  been  proposed  as  responsible  for  this  layer  is 
investigated  and  a  photochemical  model  which  describes  this  emission  has 
been  developed.  The  model  is  compared  to  observations  from  the 
Atmosphere  Explorer  satellites  to  determine  the  relative  importance  of  the 
individual  reactions.  For  some  of  the  reactions,  the  properties  of  the  rate 
coefficients,  branching  ratios,  and  cross  sections  are  not  well  known  or 
there  are  large  discrepancies  between  different  studies.  The  model  is  used 
to  determine  the  best  parameters  for  some  of  these  reactions  based  on 
observations.  The  model  is  then  applied  to  the  twilight  airglow. 

The  remainder  of  this  chapter  gives  an  extensive  account  of  the 
history  of  the  discovery  of  the  greenline  emission  and  the  determination  of 
the  responsible  transition.  The  various  proposed  chemical  reactions  for 
sources  and  losses  are  discussed  for  each  of  the  dayglow  layers  and  the 
auroral  case.  In  addition,  a  summary  of  the  observational  techniques  used 
to  study  this  emission  is  presented. 

The  Atmosphere  Explorer  (AE)  satellite  program  is  discussed  in 
chapter  two.  The  Visible  Airglow  Experiment  (VAE),  which  was  the 
instrument  on  the  AE  spacecraft  which  measured  atmospheric  emissions, 
is  described  in  detail,  together  with  an  account  of  the  various  data  analysis 
techniques  that  are  applicable  to  the  data.  Other  instruments  on  the 
satellites,  which  provided  measurements  of  the  neutral  and  ionized 


3 


atmospheric  constituents  that  are  needed  in  the  various  chemical 
reactions,  are  also  discussed. 

Details  of  the  photochemical  model  are  given  in  Chapter  3.  This 
chapter  presents  the  various  atmospheric  models  that  are  required  by  the 
photochemical  model,  followed  by  a  description  of  the  photochemical  code 
and  the  specific  reactions  that  are  employed  by  it. 

The  next  three  chapters  present  the  results  derived  from  the 
photochemical  model.  Chapter  4  details  the  calculation  of  the  total  rate 
coefficient  for  a  reaction  that  is  a  major  source  of  the  0(^S)  at  lower 
altitudes.  In  Chapter  5,  the  model  is  used  with  the  results  from  Chapter 
four  to  determine  the  effective  branching  ratio  for  the  energy  exchange 
reaction  and  the  magnitude  of  the  electron  impact  cross  section  of  atomic 
oxygen  that  best  represents  the  observed  data.  With  these  parameters 
determined.  Chapter  5  also  describes  the  application  of  the  model  to  the 
twilight  5577A  emission.  Conclusions  and  suggestions  for  further  work  are 
given  in  Chapter  six. 


1J2  Introduction  to  the  5577A  01  Emission 

The  5577A  green  line  is  produced  by  an  electronic  transition  in 
atomic  oxygen  from  the  to  the  state.  This  is  an  electric  quadrupole 
transition  which  has  a  radiative  lifetime  of  about  0.75  sec  in  the  upper  state 
(Kemahan  and  Pang,  1975).  This  emission  is  predominant  in  the  aurora 
and  is  one  of  the  brightest  emissions  in  the  visible  nightglow.  Figure  1.1  is 
an  energy  diagram  for  atomic  oxygen  showing  the  01  transitions. 


Figure  1.1  Energy  level  diagram  of  atomic  oxygen,  showing 
potential  and  observed  transitions.  Solid  arrows  are 
observed  transitions,  dashed  arrows  have  not  been 
observed.  Wavelengths  are  in  nanometers  except  as 
noted.  (After  Rundle,  1971). 
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In  the  atmosphere,  the  5577A  dayglow  emission  occurs  in  two 
distinct  layers.  The  lower  layer  is  at  about  95  km,  while  the  upper  layer 
peaks  at  about  150  km.  These  two  layers  are  shown  in  Figure  1.2.  The 
green  line  has  been  studied  more  than  any  other  airglow  emission,  and 
previous  investigations  have  been  reviewed  by  Bates  (1978)  and  Torr  and 
Torr  (1982). 

The  present  investigation  concentrates  on  the  upper  layer  which 
peaks  in  the  dayglow  at  approximately  150  km  and  has  a  complex 
photochemical  processes  for  production.  The  majority  of  the  past  studies 
were  focused  on  the  auroral  regions  and  in  the  nightglow  (Sharp  et  al., 
1979;  Sharp  and  Torr,  1979;  Rees,  1984;  Yee  and  Killeen,  1986;  McDade  and 
Llewellyn,  1986,  Gerikova  and  Shepherd,  1988).  A  study  of  the  5577  A 
airglow  including  modeling  of  the  dayglow  was  carried  out  by  Frederick  et 
al.  (1976)  using  observations  from  the  Visible  Airglow  Experiment  (VAE)  on 

board  the  Atmosphere  Explorer  satellites.  Since  then,  many  new  laboratory 
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and  observational  results  have  improved  our  understanding  of  the  5577A 
dayglow  and  the  complex  chemical  processes  which  produce  the  0(lS) 
metastable  state.  This  section  will  give  a  history  of  the  5577A  emission  and 
discuss  the  physical  processes  and  morphology  of  the  emission.  Also  the 
different  techniques  that  can  be  used  to  measure  the  emission  will  be 
described. 

1.2.1  History 

o 

The  green  line  at  5577A  has  been  studied  more  than  any  other 
emission.  As  the  brightest  emission  in  the  visible  nightglow  and  the 
dominant  color  visible  in  the  aurora,  it  attracts  much  attention.  Despite 
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Figure  1.2  The  5577A  volume  emission  rate  profile  by  Wallace 
and  McElroy  (1966).  Low  altitudes  are  corrected  for 
Rayleigh  scattering  as  indicated  by  squares  and  triangles. 
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extensive  studies,  the  major  processes  responsible  for  the  emission  are  not 
completely  resolved.  In  this  section,  an  account  of  the  various  suggested 
mechanisms  for  this  emission  will  be  given  along  with  an  historical 
perspective.  Since  the  emission  is  found  in  all  parts  of  the  sky  (aurora, 
nightglow,  and  dayglow),  no  attempt  is  made  here  to  distinguish  the 
emission  in  the  dayglow  alone.  This  will  be  done  in  later  chapters. 

Studies  of  the  "auroral”  green  line  date  back  to  1868  when  Anders.  J. 
Angstrom  (1814-1874)  began  making  quantitative  spectral  studies  of  the 
aurora.  Before  this  only  visual  color  studies  of  the  aurora  were  performed. 

o 

Angstrom  (1868)  determined  the  wavelength  of  the  green  emission  in  the 
aurora  as  5567 A  but  was  unable  to  determine  its  source.  It  was  suggested 
the  green  line  emission  would  be  present  even  though  the  aurora  was  not 
visible  (Angstrom,  1868;  Vogel,  1872,  and  Wiechert,  1902). 

It  was  noted  by  Campbell  (1895)  and  Path  (1908)  that  the  "auroral" 
green  line  emission  was  present  in  the  entire  sky  at  all  times.  This  was 
confirmed  by  Slipher  (1919)  using  a  large  number  of  green  line 
spectrograms.  This  was  termed  the  permanent  aurora  by  Yntema  (1909)  or 
the  nonpolar  aurora  by  Rayleigh  (1924).  The  terms  airglow  and  Nightglow 
were  not  introduced  until  the  early  1950's  (Elvey,  1950;  Roach  and  Pettit, 
1951). 

Early  spectroscopic  work  by  Vegard  (1913)  placed  the  green  line 
emission  near  5577k.  Using  photographs  of  the  green  line  interference 

fringes  of  the  night  sky  from  a  Fabry-Perot  etalon,  Babcock  determined  the 
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wavelength  of  the  emission  to  be  5577. 350±0. 005  A  with  a  maximum  width 
of  0.35  A.  The  results  of  these  measurements  led  to  the  conclusion  that  the 
auroral  and  nightglow  green  line  emission  had  a  common  source. 
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Sorting  out  the  contributions  to  the  night  sky  from  various  terrestrial 
and  extra-terrestrial  sources  was  a  significant  problem,  as  pointed  out  by 
Fabry  (1919,1921).  These  sources  included  the  stellar  background, 
scattering  of  sunlight  by  the  interplanetary  gas  and  dust,  and  the 
atmospheric  emission.  Polarization  of  the  night  sky  was  found  to  be  almost 
non-existant  (Rayleigh,  1919;  Babcock,  1919).  In  the  plane  of  the  azimuth  of 
the  sun,  Dufay  (1928,  1929)  found  a  polarization  of  2-4%,  which  led  to  the 
conclusion  that  about  15%  of  the  light  in  the  night  sky  was  from  the  zodical 
light. 

During  his  investigations  of  the  spectrum  of  the  night  sky,  Rayleigh 
(1922a,  1922b,  1923,  1928,  1931)  concluded  that  the  terrestrial  component  of 
the  night  sky  was  quantitatively  different  from  the  aurora.  The  night  sky 
lacked  the  negative  bands  of  nitrogen  seen  in  the  aurora,  and  the  green  line 
showed  no  consistent  enhancement  toward  the  auroral  regions.  The  first 
absolute  measurements  of  green  line  brightness  were  reported  by  Rayleigh 
(1930).  These  measurements  in  terms  of  surface  brightness  led  to  the 
photometric  unit  named  after  him  (Hunten  et  al.,  1956). 

Even  though  the  bright  green  line  emission  was  predominant  in  all 
parts  of  the  sky,  its  source  was  not  known.  There  were  several  sources 
proposed  and  then  discarded  in  the  process  of  determining  the  mechanism 
for  the  emission.  At  the  turn  of  the  century,  it  was  suggested  that  krypton 
gas  was  responsible  for  the  green  line.  Krypton  gas  does  have  a  line  at 
5570. sA  along  with  several  others  close  to  prominent  nitrogen  lines,  and  it 
was  proposed  that  krypton  could  explain  the  entire  auroral  spectrum 
(Ramsay,  1919). 

An  unknown  gas  called  geocoronium  was  suggested  by  Wegener 
(1911)  as  the  source  of  the  green  line.  This  gas  was  supposed  to  be  lighter 
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than  nitrogen,  originated  in  the  solar  atmosphere  and  was  transported  to 
the  earth  (Henriksen  and  Egeland,  1988).  This  hypothesis  did  not  last  when 
the  implied  characteristics  of  the  gas  were  compared  to  the  aurora.  If  the 
gas  was  lighter  than  nitrogen,  then  the  fall  off  of  the  green  line  with  height 
would  be  smaller  than  that  for  the  nitrogen  emissions  (Vegard,  1923a). 
Also,  Vegard  pointed  out  that  there  was  no  place  in  the  periodic  table  for 
such  an  element. 

In  the  early  1920's,  it  was  suggested  by  Vegard  (1923a,  1923b,  1924) 
that  the  source  of  the  green  line  emission  was  solid  nitrogen.  He  found  that 
bombarding  frozen  nitrogen  crystals  with  accelerated  electrons  produced  a 
series  of  emissions  coinciding  with  an  auroral  spectrum  (Vegard,  1924). 
Also,  he  found  a  broad  feature  around  the  5577A  emission.  This  broad 
feature  was  found  to  consist  of  emissions  at  5556A,  5572 A,  56 10  A  and  5654A 
but  a  line  at  5577A  was  not  observed  (McClennan  and  Shrum,  1924). 
Further  experiments  by  McClennan  and  Shrum  (1925)  and  McClennan  and 
McLeod  (1927)  found  that  the  green  line  was  due  to  atomic  oxygen  via  a 
previously  unknown  transition. 

The  nature  of  this  unknown  transition  was  first  hinted  at  by  Bowen 

(1927,  1928)  through  his  studies  of  the  spectra  of  planetary  nebula.  He 

invoked  the  idea  of  forbidden  transitions  and  in  the  following  quote  set  up  a 

key  premise  in  the  development  of  thermospheric  photochemistry. 

In  the  spectra  of  gaseous  nebulae  several  very  strong  lines  are 
formed  which  have  not  been  duplicated  in  any  terrestrial 
source.  Many  lines  of  evidence  point  to  the  fact  that  these  lines 
are  emitted  by  an  element  of  low  atomic  weight.  Since  the 
spectra  of  light  elements  as  excited  in  terrestrial  sources  are 
well  known,  this  leads  to  the  conclusion  that  there  must  be 
some  condition,  presumably  low  density,  which  exists  in  the 
nebulae,  that  causes  additional  lines  to  be  emitted. 
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From  studies  of  the  Zeeman  pattern  it  was  concluded  that  the  5577A 
emission  was  from  a  forbidden  and  previously  unknown  (ID- IS)  transition 
of  atomic  oxygen  (McLennan,  McLeod,  and  McQuarrie,  1927;  McLennan, 
McLeod,  and  Ruedy,  1928;  McLennan,  1928;  and  Sommer,  1928).  Frerichs 
and  Campbell  (1930)  confirmed  the  identification  of  the  green  line  as 
arising  from  an  electric  quadrupole  transition.  Frerichs  (1930)  used  the 
ultraviolet  spectrum  of  oxygen  to  predict  the  energy  of  the  ID2  and  ^Sq  states. 
With  the  identification  of  the  source  and  transition  for  the  green  line 
emission,  the  next  issue  was  to  determine  the  mechanisms  that  produce 
0(lS)  in  the  upper  atmosphere. 

1.2.2  Physical  Processes 

After  the  identification  of  the  source  of  the  5577A  emission  as  being 
due  to  the  (ID-IS)  transition  of  atomic  oxygen.  Chapman  (1931)  made  the 
first  attempted  theoretical  explanation  for  the  feature  using  a  process  first 
proposed  by  Frank  in  1929  (as  cited  by  Bates,  1981).  Chapman  noted  that 
McLennan  (1930)  and  Rayleigh  (1931)  observed  that  the  emission  at  5577A  is 
essentially  constant  throughout  the  night.  This  meant  that  the  energy 
involved  is  comparable  to  the  estimate  of  the  energy  needed  to  ionize  the 
upper  layer.  Chapman  (1931)  argued  that  to  maintain  the  emission  at  a 
constant  level  throughout  the  night,  the  number  of  transitions  producing 
5577A  photons  must  be  small  compared  to  the  number  of  precursor 
particles. 

Chapman  looked  to  the  reservoir  of  energy  from  the  dissociation  of  O2 
as  the  source  of  0(^8).  He  first  noted  that  the  O2  association  energy  of  5.03 
eV  was  sufficient  to  provide  the  excitation  energy  of  4.2  eV  for  0(^8).  It  is 
clear  from  his  description,  although  he  never  explicitly  wrote  the  reaction. 
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that  the  three-body  reaction  that  has  become  known  as  the  Chapman 
mechanism  (Equation  [1.1])  was  responsible,  i.e. 

0('P)  -»■  0('P)  0('P)  — O2  +  OCS)  ,  [1.1] 

where  Yj  is  the  reaction  rate  coefficient. 

Based  on  the  estimated  density  for  the  0(3p),  the  height  of  the 
emission  due  to  the  Chapman  mechanism  was  placed  at  about  100  km. 
This  altitude  was  not  supported  by  many  observations  using  the  van  Rhijn 
method  (van  Rhijn,  1921;  see  Chamberlain,  1961),  which  placed  the  altitude 
much  higher.  The  van  Rhijn  method  makes  use  of  the  fact  that  the  precise 
changes  in  intensity  with  zenith  angle  depend  on  the  height  of  the  emitting 
layer.  The  method  makes  the  assumptions  of  an  optically  and  linearly  thin 
emitting  layer  that  is  homogeneous  and  spherically  symmetric.  There  is 
no  provision  for  extinction  in  the  lower  atmosphere  for  the  emission  in  the 
method.  The  invalidity  of  any  of  these  assumptions  cause  the  van  Rhijn 
method  to  give  widely  varying  results. 

The  van  Rhijn  method  was  shown  later  to  be  unreliable  for  height 
determinations,  but  the  matter  was  not  resolved  for  many  years.  Phillips 
(1956)  used  a  Fabry-Perot  interferometer  to  determine  the  Doppler  width  of 

0 

the  5577 A  emission,  and  found  it  would  correspond  to  a  temperature  of  155 
to  231  K,  which  is  consistent  with  altitudes  of  about  100  km.  This  was 
followed  by  rocket  observations  that  placed  the  peak  of  the  layer  between  90 
and  102  km  (Berg  et  al.,  1956;  Koomen  et  al.,  1956;  others  summarized  by 
Nardi,  1991). 

The  rocket  measurements  of  Wallace  and  McElroy  (1966)  established 
that  the  green  line  actually  had  two  distinct  layers.  In  addition  to  the  layer 
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at  about  95  km,  there  is  an  upper  F-region  layer.  Figure  1.2  shows  the 
results  of  their  rocket  flight.  The  rest  of  this  section  will  discuss  these  two 
layers  separately  beginning  with  the  lower  layer.  For  completeness,  a 
discussion  of  auroral  mechanisms  will  be  presented  last. 

Lower-altitude  component.  The  Chapman  mechanism  as  the 
source  of  0(^8)  was  not  seriously  questioned  for  some  time.  The  estimation 
for  the  rate  coefficient,  yi,  is  difficult  to  make  by  either  experiment  or 
theory.  Also,  it  was  believed  that  only  a  small  number  of  the  associations 
formed  would  need  to  produce  0(^8)  to  give  the  measured  intensity.  Nicolet 
(1959)  deduced  that  only  one  in  one  thousand  associations  need  produce  a 
5577A  photon  to  account  for  the  measurements  of  Rayleigh  (1930).  At  this 
time  quenching  of  0(^8)  was  thought  to  be  unimportant  due  to  the 
investigation  by  Kvifte  and  Vegard  (1947)  that  gave  the  rate,  P2.  for 
equation 


0(‘8)  +  O2  — ^  OCP)  +  O2  ,  [1.2] 

to  be  4.0  X  10  * cm3  3-1  for 

0('8)  +  OCP)  — ^  OCP  or  'D)  +  OCP  or  *D),  [1.3] 

as  4.  X  10  '13  cm3  g-i 

As  interest  in  the  airglow  emission  grew,  other  sources  of  0(^8)  were 
suggested.  The  choices  were  limited  by  the  need  to  provide  4.2  eV  excitation 
energy  to  0(3p).  Excess  vibrational  energy  was  suggested  as  a  source  by 
Krassovsky  ( 1958)  in  the  reactions 
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o;  +  CO  — ^  CO2  +  o('s)  [1.4] 

OH’  +  OH  — ^  H2O  +  0('S)  ,  [1.5] 

where  O2*  and  OH*  are  vibrationally  excited.  This  was  also  the  case  for  the 
reaction 


NO’  +  N  N2  +  O('S)  ,  [1.6] 

proposed  by  Young  and  Clark  (1960).  These  were  refuted  by  Dalgamo  (1963) 
mainly  on  the  grounds  that  the  vibrational  energy  would  mostly  be 
converted  to  translational  energy  in  atom-molecule  collisions  (Bates,  1955). 

Problems  with  the  Chapman  mechanism  as  the  source  of  the  lower 
0(^S)  layer  surfaced  due  to  two  developments.  The  first  was  that  most  of  the 
oxygen  association  was  found  to  be  due  to  the  catalytic  action  of  hydrogen 
(Nicolet,  1964,  1971;  Thomas  and  Bowman,  1972),  and  the  density  of  atomic 
oxygen  in  the  region  was  less  than  supposed.  This  would  make  the  number 
of  associations  with  oxygen  as  the  third  body  considerably  larger.  The 
second  was  of  more  importance.  The  rate  coefficient  for  Reaction  [1.1]  (  yi ) 
was  determined  to  be  less  than  8.0  x  lO'^®  cm®  s'^,  and  quenching  of  0(^S)  by 
O2,  Reaction  [1.2],  was  2.7  x  IQ  l®  cm^  s-l  (Barth  and  Hildebrandt,  1961). 
The  agreement  with  the  value  from  Kvifte  and  Vegard  (1947)  for  the 
quenching  rate  gave  some  confidence  in  the  value  for  yi,  which  is  lower 
than  the  value  required  to  explain  the  nightglow  emission. 
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These  problems  with  the  Chapman  mechanism  led  to  the  suggestion 
of  a  two  step  mechanism  with  excited  O2  (Barth,  1961,  1962,  1964).  The  two 
step  process  starts  with  an  association  reaction, 

0('P)  +  0('P)  +M  — ^  0;  +  M  .  [1.7] 

This  reaction  would  then  be  followed  by  an  excitation  transfer  process, 

O;  +  0('P)  — ^  O2  +  OCS)  .  [1.8] 

The  excited  O2  from  Equation  [1.7]  is  thought  to  be  in  one  of  four  possible 
states:  c^Z„  ,  ,  A^I^,  or  .  These  states  have  been  debated  in  several 

reviews  (  see  Bates,  1978,  1981,  1988;  Torr  and  Torr,  1982)  and  will  be 
discussed  later. 

The  loss  of  O2*  given  in  Equation  [1.8]  is  in  competition  with  several 
other  loss  mechanisms.  The  radiative  loss  is  given  by, 

0;  --g-  >  O2  +  hu  .  [1.9] 

Other  quenching  losses  would  be  given  by 

0;  +  M  — ^  O2  +  M  ^  [110] 

where  M  is  not  0(3P)  and 


O.;  +  0('’P) 


+  0("Por  ’D)  . 


[1.11] 
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This  set  of  reactions  allows  an  effective  production  rate  for  OC^S)  to  be 
determined  as  where  Yb  is  , 


_ Y7ps[M] _ 

(p3  +(3„XO]  +  p.o[M]  +  A,  ’ 


[1.12] 


as  given  by  Bates  (1979). 

The  greenline  photon  emission  rate  is  less  than  the  excitation  rate 
from  either  the  Chapman  or  Barth  mechanisms  due  to  collisional 
deactivation  of  the  OC^S).  Particular  attention  was  given  by  Atkinson  and 
Welge  (1972)  to 


0('S)  +  N.,  0('P  or  ‘D)  +  N2  and 


[1.13] 


0(‘S)  +  O2  — 0(®P  or  'D)  +  O2  . 


[1.14] 


They  found  P13  to  have  an  upper  limit  of  5.  x  cm^  s'l,  in  agreement  with 
Young,  Black,  and  Slanger  (1969).  For  P14  they  determined  a  value  of 


P,4  =  4.9  X  lO'^exp 


r-1700 
RT 


cm  s 


[1.15] 


which  is  consistent  with  the  value  found  by  Zipf  (1967)  at  300  K  of  2.1  x  lO'l^ 
cm3  g-l  This  was  confirmed  by  Slanger,  Wood,  and  Black  (1972)  with  a 
value  of 


p,4  =  4.0  X  10  exp 


-1730 

RT 


cm^  s  ' 


\ 


\ 


[1.16] 
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Slanger  and  Black  (1981)  identified  the  state  of  the  O2  in  the  quenching 
reaction  to  be  02(^Ag). 

Deactivation  of  0(^S)  by  0(3P)  was  proposed  by  Olsen  (1973)  to  occur 
mainly  by 


O('S)  +  OC'’P)  — 0(*D)  +  0(*D)  +  2.55  eV  .  [1.17] 


While  the  earlier  measurements  by  Kvifte  and  Vegard,  1947)  gave  a  value  of 
4  X  10  cm^  s'l  regardless  of  the  state  of  the  products,  Felder  and  Young 
(1972)  found  the  rate  to  be  7.5  x  10  cm^  s‘l  at  300  K.  Measurements  by 
Slanger  and  Black  (1976a)  for  Equation  [1.17]  gave  P17  to  be 


Pi7  =  5.0  X 


10  "exp 


-610^ 
RT  . 


cm  s  , 


[1.18] 


which  agrees  well  with  a  single  measurement  made  by  Lorents  and 
Huestis  (1975). 

On  the  basis  of  measured  profiles  of  0(3p)  and  laboratory  rate 
coefficients,  Slanger  and  Black  (1977)  concluded  that  only  the  Barth 
mechanism  would  reproduce  the  observed  0(3p)  data.  A  study  by  Thomas 
et  al.  (1979)  supported  these  conclusions  and  a  theoretical  study  led  Bates 
(1979)  to  present  arguments  supporting  the  Barth  mechanism  over  the 
Chapman  mechanism.  The  study  of  Thomas  (1981)  used  simultaneous 

0 

nighttime  measurements  of  atomic  oxygen  densities,  5577A  greenline 
emission  and  the  Hertzberg  I  bands  of  O2  to  support  the  Barth  mechanism 

and  show  that  0(3p)  may  not  be  the  predominant  quencher  of  0(3S). 
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Since  the  Barth  mechanism  is  the  generally  accepted  method  of 
production  of  the  0(^S)  in  the  lower  thermosphere,  attention  has  been  given 
to  determining  which  of  the  four  possible  O2*  states  is  produced  in  Equation 
[1.7].  Of  these  states,  Thrush  and  Stott  (1987)  have  evidence  which  shows 
that  02(A,v>6)  can  be  a  precursor  for  OC^S)  but  02(A,v<6)  and  02(A’,v=2,3,4) 
are  not.  This  along  with  the  oxygen-argon  afterglow  studies  of  Slanger  and 
Black  (1976b)  precludes  these  two  states  as  being  possible  precursors  for 
0(lS).  The  measurements  of  Slanger  and  Black  (1976b)  indicate  support  for 
02(c)  as  the  precursor  in  the  afterglows  and  are  consistent  with  it  being  the 
nightglow  precursor. 

The  02(^ng)  state  is  favored  as  the  precursor  by  Wraight  (1982)  and 
Krasnopolsky  (1986).  The  02(c)  state  is  thought  to  be  the  precursor  since  it 
can  produce  0(1S)  for  02(c,v>2).  In  comparison,  the  Venus  atmosphere 
where  the  greenline  emission  is  very  weak,  02(c)  is  in  the  v  =  0  level.  But, 
02(®ng)  can  meet  the  requirements  for  the  precursor  in  the  atmospheres  on 
Venus  and  the  Earth  if  the  quenching  for  the  precursor  by  CO2  is  very 
effective  (Krasnopolsky,  1986).  While  02(^ng)  is  not  significant  as  a 
precursor  in  the  afterglows  of  Slanger  and  Black,  its  importance  in  the 
nightglow  as  a  precursor  of  0(^S)  is  not  known  and  has  no  obvious  difficulty 
except  the  lack  of  positive  evidence  (Bates,  1988). 

The  Photodissociation  of  O2  by  radiation  short  of  1334A  could  provide 
an  additional  source  of  0(^8)  in  the  dayglow.  Early  estimates  for  the 
dayglow  were  given  as  5  Rayleighs  (R)  (Bates  and  Dalgamo,  1954),  with  an 
upper  limit  of  15  R  (Walker,  1965).  Wallace  and  McElroy  (1966)  using  data 
from  Metzger  and  Cook  (1964)  placed  the  resulting  luminosity  below  150 
km.  Hays  and  Sharp  (1973)  found  better  agreement  between  120  km  and  150 
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km  for  twilight  rocket  data  by  using  this  additional  source.  This  source  will 
be  discussed  further  in  the  section  on  the  upper  layer  of  5577A  emission. 

Higher-altitude  component.  The  higher  altitude  layer  of  the  5577A 
emission  peaks  in  the  F-region  of  the  ionosphere.  In  addition  to  being 
much  higher  than  the  lower  peak  it  is  weaker  and  much  broader, 
especially  at  night.  Unlike  the  lower  layer,  which  is  attributed  to  a  single 
source  mechanism,  the  higher  layer  has  several  source  mechanisms  that 
are  effective  over  various  altitude  ranges.  Each  of  these  proposed  reactions 
will  now  be  discussed. 

Photoelectron  impact  on  atomic  oxygen, 

0(=*P)  +  ep-li2_>0('S)  +  ep  ,  [1.19] 

is  an  important  source  of  0(1S)  in  the  thermosphere  (Wallace  and  McElroy, 
1966).  In  fact,  this  is  a  dominant  source  of  0(^8)  in  the  dayglow.  The 
production  rate  for  this  reaction  is  a  function  of  atomic  oxygen  density,  [0]  , 
given  by: 


P,9=[0]Jf,(EJa(EJdE,  .  [1-20] 

The  integral  over  electron  energy  Ee  contains  the  photoelectron  flux,  fe,  and 
a,  the  cross  section  for  electron  impact,  both  of  which  are  functions  of 
electron  energy. 

Feldman  et  al.  (1971)  used  simultaneous  rocket  measurements  of 

0 

photoelectrons  and  the  5577 A  emission  made  by  Doering  et  al.  (1970),  and 
showed  this  reaction  to  be  minor  below  250  km  and  to  account  for  half  of  the 
emission  above  250  km.  Twilight  studies  by  Hays  and  Sharp  (1973)  showed 
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the  electron  impact  source  to  be  important  above  190  km  and  dominant 
above  260  km.  This  was  also  supported  by  Frederick  et  al.  (1976)  using 
multiparameter  Atmosphere  Explorer  Satellite  observations. 

The  cross  section  for  electron  impact  excitation  of  OC^S)  has  been 
measured  by  Shyn  e/ a/.  (1986)  and  Doering  and  Gulcicek  (1989).  The 
magnitudes  of  the  cross  section  obtained  by  Shyn  et  al.  are  higher  than  the 
previous  theoretical  results  (  Smith  et  al.,  1967;  Henry  et  al.,  1969;  Thomas 
and  Nesbet,  1975  and  Vo  Ky  Lan  et  al.,  1972)  by  a  factor  of  about  2,  while 
those  by  Doering  and  Gulcicek  (1989)  are  higher  by  a  factor  of  approximately 
1.5.  Parameterized  values  of  the  cross  section  for  the  excitation  of  0(lS) 
determined  from  the  values  by  Smith  et  al.  (1967)  and  later  supported  by 
Henry  et  al.  (1969)  have  been  given  by  Green  and  Stolarski  (1972).  Recent 
theoretical  calculations  by  Herrington  and  Burke  (unpublished  data)  are  in 
good  agreement  with  those  of  Henry  et  al.  (1969)  and  have  uncertainties  of 
10-20%  (Berrington,  1991).  Figure  1.3  is  a  plot  of  the  theoretical  and 
laboratory  values  of  the  electron  impact  cross  sections  for  atomic  oxygen. 

The  dissociative  recombination  reaction  for  02‘''. 

Q2+e,-^^’-»0(‘S)  +  0  ,  [1.21] 

is  known  to  be  an  important  contributor  to  the  green  line  emission  at  higher 
altitudes  and  the  only  source  in  the  nightglow.  Dissociative  recombination 
was  suggested  as  an  important  process  to  explain  the  observed  nightglow 
spectrum  by  Kaplan  (1931).  The  general  mechanism  for  this  reaction  was 
discussed  by  Bates  and  Massey  (1947),  who  described  the  observed  electron 

o 

recombination  rates.  The  importance  to  5577A  emission  in  the 
thermosphere  was  suggested  by  Nicolet  (1954). 


Impact  Cross  Section  (cm  ) 


20 


1  10  100 

ENERGY  (eV) 


Figure  1.3  Various  cross  sections  for  electron  impact  on  atomic 
oxygen.  Lines  are  theoretical  determinations  while 
symbols  are  experimental. 


This  reaction  has  been  studied  extensively  in  the  past  by  using 
thermospheric  nightglow  measurements  (Hernandez,  1971;  Hays  and 
Sharp,  1973;  Frederick  et  al.,  1976;  Kopp  et  al,  1977;  Sharp  and  Torr,  1979; 
O'Neil  et  al.,  1979;  Zipf,  1979;  Abreu  et  al.,  1983;  Kilieen  and  Hays,  1983;  Yee 
and  Killeen,  1986;  Yee  et  al.,  1989;  Takahashi  et  al.,  1990)  where  it  is  the 
only  production  mechanism  for  5571 A  emission.  Numerous  laboratory 
studies  have  been  carried  out  for  this  mechanism  (ivlehr  and  Binndi,  1969; 
Walls  and  Dunn,  1974;  Zipf,  1980;  Bates  and  Zipf,  1980;  Guberman,  1987; 
Rowe  and  Queffelec,  1989;  Queffelec  ei  aZ.,  1989;  Guberman,  1989  ).  The 
amount  of  OC^S)  produced  by  the  dissociative  recombination  is  dependent  on 
the  vibrational  excitation  of  the  02'''  and  is  not  conclusively  known.  The 
various  studies  and  their  findings  have  been  reviewed  and  summarized  by 
Bates  (1990). 

The  production  rate  for  0(^8)  from  this  reaction  is, 


P2i=p2.k2,[0;][eJ  . 


The  total  reaction  rate  as  a  function  of  electron  temperature. 


[1.22] 


k^i  =  16x  10 


-7 


300 


0.55 


[1.23] 


is  given  by  TorreZ  al.  (1976).  This  expression,  derived  from  measurements 
taken  by  the  Atmosphere  Explorer  Satellite  AE-C,  has  been  shown  to  agree 
well  with  the  measurements  of  Walls  and  Dunn  (1974)  and  other  laboratory 
measurements  as  discussed  by  Walls  and  Dunn  (1974).  Mul  and  McGowan 
(1979)  concluded  that  the  total  rate  coefficient  is  not  substantially  affected  by 
the  vibrational  excitation  of  the  02"^  ion. 
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The  branching  ratio  for  this  reaction,  P21,  is  not  well  known.  Early 
results  of  aeronomical  studies  under  various  conditions  varied  but  were 
typically  less  than  0.10  (Hernandez,  1971;  Hays  and  Sharp,  1973;  Frederick 
et  al.,  1976;  Kopp  et  al.,  1977;  Sharp  and  Torr,  1979;  O'Neil  et  al.,  1979;  Zipf, 
1979).  These  determinations  of  the  production  of  0(1S)  are  not  in  agreement 
with  the  ab  intio  calculations  of  Guberman  (1987,  1989)  (Zipf,  1988;  Yee  et 
al.,  1989;  Bates,  1990).  Paxton  (1983)  used  the  values  of  the  branching  ratios 
known  on  Earth  to  deduce  unknown  branching  ratios  for  other  channels 
from  observations  of  the  upper  atmosphere  of  Venus. 

From  theoretical  considerations,  the  major  channel  for  the 
production  of  0(1S)  is  from  the  v=2  level  of  the  02''’  ion  (Guberman,  1987; 
1989).  This  was  supported  by  Killeen  and  Hays  (1983)  and  Yee  and  Killeen 
(1986)  from  studies  using  the  Fabry-Perot  on  the  Dynamics  Explorer 
satellite  to  study  the  nonthermal  line  profiles  of  OC^S).  The  quantum  yield  of 
0(^8),  P21,  as  a  function  of  the  vibrational  population  of  02^,  can  be 
expressed  in  terms  of  the  ratio  of  electron  density  to  atomic  oxygen  , 


P2. 


oc 


[O] 


300 


,0.7 


[1.24] 


as  discussed  in  Abreu  et  al.  (1983)  and  Yee  et  al.  (1989). 

The  energy  transfer  reaction  of  N2  (A^Lu^)  by  atomic  oxygen, 

N,(A^z;)+0("P)-^^N2+0('S)  ,  [1.25] 


can  produce  a  major  portion  of  the  0(1S)  at  altitudes  near  the  peak  in  the 
thermospheric  emission  layer.  First  proposed  as  an  auroral  source  by 
Vegar  '  and  Tonsberg  (1937)  and  later  rejected  by  Vegard  (1939),  this 
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reaction  was  not  considered  to  be  important  until  the  laboratory 
measurements  by  Meyer  et  al.  (1969,1970).  As  a  source  of  OC^S),  this 
reaction  was  considered  to  be  a  minor  source  by  Frederick  et  al.  (1976).  The 
total  rate  coefficient  used  was  3.0  x  10'^^  cm^  s*^  (Hunten  and  McElroy  , 
1966;  Meyer  et  al.,  1970).  The  value  of  the  branching  ratio  in  Frederick's 
model  was  0.10,  which  was  chosen  to  be  consistent  with  auroral 
observations  at  that  time  (Frederick  et  al.,  1976). 

The  N2  (A)  state  is  produced  both  directly  by  photoelectron  impact  and 
indirectly  by  cascading  from  higher  excited  states  of  ihe  nitrogen  molecule. 
The  effective  rate  of  production  of  0(lS)  from  this  reaction, 

P25  =  p25k26[N,(A^i:)][0],  [1.26] 

can  be  expressed  as  a  product  of  the  total  rate  coefficient  for  the  production 
of  all  states  of  atomic  oxygen  and  the  quantum  yield  or  branching  ratio  for 
the  amount  produced  in  the  0(1S)  state.  The  branching  ratio,  P25,  is  not 
well  known  and  there  is  still  some  disagreement  in  the  total  rate  for  this 
reaction,  k25,  as  well.  Table  1.1  is  a  list  of  some  of  the  rate  coefficients  and 
branching  ratios  for  this  reaction  that  have  been  determined  by  previous 
studies. 

Total  reaction  rates  measured  by  Piper  et  al.  (1981)  for  the  quenching 
of  the  u  =  0,1  states  of  N2(A)  by  0  and  recent  measurements  by  Thomas  and 
Kaufman  (1985)  for  the  u  =  0, 1,2,3  states  of  N2(A)  show  that  as  the 
vibrational  level  increases,  the  rate  constant  for  the  energy  transfer 
reaction  also  increases.  The  two  production  mechanisms  of  N2(A)  tend  to 
fill  different  vibrational  levels  in  the  excited  state.  Direct  excitation  fi-om  the 
ground  state  of  the  N2  molecule  by  electron  impact  will  preferentially  fill  the 


Branching  Ratio 

Study 

Comments 

3.0 

0.10 

Frederick  et  al. 

Branching  ratio  set 
to  match  auroral 
detennination 

u=0  3.5±0.6 
u=l  4.110.5 

u=2  4.610.6 
u=3  5.210.8 

Thomas  and 
Kaufmann 

Laboratory  study 

\)=0  2.810.4 
u=l  3.410.6 

0.7510.13 

Piper  etal. 

Laboratory  study 

u=0  2.8 
u=l  3.3 

\)=2  3.6 

0.75 

De  Souza  et  al. 

Laboratory  study 

0.57 

0.30 

O'  Neil  et  al. 

12.0 

0.30 

0.45 

Sharp  and  Torr 

For  v=0, 1  1 

For  0=0, 1,2  I 

Table  1.1  List  of  various  rate  coefficients  and  branching  ratios 
for  Reaction  [1.25]. 
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higher  vibrational  levels,  while  cascading  from  the  state  and  above 

tend  to  populate  the  lower  levels  (Ahmed,  1969;  Cartwright  et  al., 
1971,1973;Cartwright,  1978).  The  total  rate  coefficient  for  the  reaction  will 
depend  on  the  vibrational  population  distribution  of  the  N2(A)  molecules, 
regardless  of  the  final  state  of  atomic  oxygen  produced. 

Using  laboratory  measurements.  Piper  (1982)  determined  a 
branching  ratio  for  the  production  of  0(^8)  of  0.75±0.13  from  the  u  =  0,1 
states  of  N2(A).  This  value  is  supported  by  laboratory  measurements 
carried  out  by  De  Souza  et  al.  (1985).  Use  of  the  laboratory  determined 
branching  ratios  under  auroral  conditions  would  produce  an  excessive 
amount  of  0(^8),  as  pointed  out  by  McDade  and  Llewellyn  (1984).  The 
quantum  yield  of  this  reaction  of  about  0.30  has  been  deduced  by  O'Neil  et  al. 
(1979)  for  auroral  production  during  the  artificial  auroral  experiment, 
Precede.  Studies  by  Sharp  and  Torr  (1979)  supported  this  value  for  the  v  = 
0,1  levels  of  N2(A),  and  by  including  the  u  =  2  vibrational  level  they  foimd  the 
branching  ratio  increased  to  0.45. 

Photodissociation  of  O2,  whose  possible  importance  was  pointed  out 
by  Bates  and  Dalgarno  (1954), 

0,  +  hv-^""  >0(^8)  + O  .  [1.27] 

is  a  very  minor  source  for  the  production  of  0(^8)  at  altitudes  above  140  km. 
However,  this  reaction  becomes  the  dominant  source  at  120  km  (Hays  and 
Sharp,  1973).  The  cross  sections  for  this  reaction  have  been  given  by 
Lawerence  and  McEwan  (1973).  The  rate  for  this  reaction  as  a  function  of 
O2  density  is. 
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P,7  =  [OJ j f(X)q(X)a(X)dX  ,  [1-28] 


where  ftX.)  is  the  attenuated  solar  flux,  a(X)  is  the  photodissociation  cross 
section,  and  q(A,)  is  the  quantum  yield  of  0(^S)  all  as  a  function  of 
wavelength,  X.  An  expression  for  the  attenuated  solar  flux  is 


f(>.)  =  f«(?^)exp 


[-2J 


Ogi(X,)nj(s)  ds 


[1.29] 


where  s  is  the  distance  along  the  path  from  the  sun  to  the  point  in  the 
atmosphere,  foo  (X)  is  the  unattenuated  solar  flux,  <7ai(^)  is  the  total 
absorption  cross  section,  and  ni(s)  is  the  number  density  of  the  ith  species 
along  the  path.  The  species  included  in  the  solar  flux  attenuation  are 
atomic  oxygen,  O2  and  N2. 

Ion-atom  interchange  between  02'''  and  atomic  nitrogen, 

0;  +  N— >0(*S)  +  N0^  ,  [1.30] 


was  considered  by  Frederick  ef  al,  (1976)  as  the  required  additional  source  of 
0(lS)  at  low  altitudes.  The  rate  of  production  for  this  process  is. 


P3o=P3ok3o[0;][N].  [1.31] 

Atomic  nitrogen  was  found  to  have  a  diurnal  variation  that  matched  the 
variation  of  0(lS)  observed  at  sunset  and  sunrise  (Kopp  et  al.,  1977).  A  rate 
coefficient,  kao,  of  2.5  x  10’^^  cm^  s*^  was  used  by  Frederick  et  al.  (1976)  for 
the  analysis  of  the  5577A  dayglow.  The  uncertainties  in  the  nitrogen 
density  gives  a  range  for  the  rate  coefficient  from  1.2  x  lO'^l  to  5.0  x  lO'^^ 
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cm3  g-l  (Frederick  et  al.,  1976).  This  reaction  may  or  may  not  be  significant 

0 

for  the  production  of  5577A  emission  depending  on  the  solar  flux  and  the 
impact  cross  section  for  electrons  on  atomic  oxygen. 

Another  ion  exchange  reaction,  between  N"*"  and  O2, 

N"  +  O2 - >  NO^  +  0('S)  ,  [1.32] 

was  suggested  by  Strobel  et  al.  (1975).  It  is  unlikely  that  this  reaction  will 
produce  any  appreciable  amount  of  0(^8)  (Tully  et  al.,  1971).  The  N+  and  O2 
exchange  mechanism  has  a  total  rate  coefficient  of  6  x  10"  10  cm3  g-l 
(Howorka  et  al.,  1980)  and  43  ±  5  %  of  the  reactions  form  atomic  oxygen 
(Smith  et  al.,  1978;  Howorka  et  al.,  1980).  However,  less  than  0.1  %  of  the 
atomic  oxygen  is  in  the  form  of  0(1S)  (Langford  et  al.,  1986). 

Other  mechanisms  that  have  been  investigated  or  suggested  as 
possible  sources  for  the  upper  layer  include, 

N('D)  +NO - >  N2  +  0('S)  ,  [1.33] 

which  was  investigated  by  Frederick  et  al.  (1976)  and  then  rejected  since  it 
did  not  fit  the  observations.  For  the  region  below  120  km,  the  source 
responsible  for  the  lower  layer  was  suggested  but  due  to  the  low  densities  it 
is  unlikely  that  either  the  Chapman  or  the  Barth  mechanism  will  produce 
any  appreciable  amount  of  5577k  emission.  The  Chapman  mechanism  is 
especially  doubtful  since  it  also  has  a  small  rate  coefficient. 

The  loss  processes  of  0(1S)  at  the  level  of  the  thermospheric  emission 
are  the  radiative  transitions  to  the  Of^D)  and  0(3p)  states, 
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0('S)— ^^0('D)  +  hv(5577A)  [1.34] 

0(*S)— ^0("P)  +  hv(2972A)  .  [135] 

with  transition  probabilities  of  1.06  s'^  and  0.045  s‘l,  respectively  (Kernahan 
and  Pang,  1975).  Losses  of  0(^8)  due  to  quenching  are  not  important  at  the 
altitudes  of  the  upper  thermospheric  emission.  The  same  quenching 
reactions  (Equations  [1.13],  [1.14],  and  [1.17])  as  occur  in  the  lower  layer 
may  have  some  effect  below  125km. 

Auroral  0(^S).  The  study  of  the  5577k  emission  in  the  airglow 
stemmed  from  observational  investigations  of  the  aurora.  The  emission  in 
the  aurora  has  similar  source  reactions  to  the  upper  airglow  layer  also 
includes  the  quenching  of  the  lower  layer.  This  section  briefly  summarizes 
the  auroral  emission  as  it  compares  to  the  other  airglow  reactions. 

While  electron  impact  on  atomic  oxygen  is  a  major  source  of  0(^8)  in 
the  daytime  airglow,  it  is  usually  of  minor  importance  in  aurora.  In 
auroral  spectra,  both  the  greenline  at  5577A  and  the  red  line  at  6300A  of 
atomic  oxygen  are  strong.  In  laboratory  studies,  however,  the  red  line  is 
absent,  which  led  Vegard  (1932)  to  doubt  the  direct  mechanisms  for 
production  of  0(^8).  Based  on  the  present  knowledge  of  impact  cross 
sections,  Gattinger  et  ai,  (1985)  concluded  that  direct  electron  impact 
(Equation  [1.19])  is  of  minor  importance  in  the  aurora. 

Dissociative  recombination  of  02"^  (Equation  [1.21])  was  suggested  as  a 
small  source  by  Dalgamo  and  Khare  (1967)  when  it  was  shown  that  the  loss 
of  slow  secondary  electrons  to  thermal  electrons  in  the  aurora  could 
compete  strongly  with  Reaction  [1.19].  On  the  basis  of  two  rocket  flights. 
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Donahue  et  al.  (1968)  found  the  secondary  electron  flux  to  be  much  too  small 
to  account  for  the  measured  5577A  emission.  They  examined  the  possible 
source  from  fast  electron  impact  dissociation  of  molecular  oxygen, 

O,  +  ep  - >  O  +  0(’S)  .  [1.36] 

They  also  found  that  this  reaction  would  require  an  unrealistically  large 
cross  section  of  -10*^®  cm^,  while  Zipf  had  determined  the  cross  section  to  be 
less  than  2  x  10  cm^  (Vallance  Jones,  1974).  Donahue  et  al.  (1968) 
concluded  that  dissociative  recombination  (Equation  [1.21])  was  the  major 
source.  Using  rocket-borne  instruments  Parkinson  et  al.  (1970)  found  that 
Reaction  [1.21]  could  explain  the  5577A  emission  at  the  altitude  of  the  FI 
layer,  but  near  the  peak  of  the  emission  there  was  insufficient  O2'*'  to 
account  for  the  emission. 

This  led  to  the  consideration  of  indirect  sources  of  0(^S).  It  was 
shown  by  Olmholt  (1971)  that  any  indirect  source  will  have  a  time  constant 
of  the  order  of  that  for  radiative  decay.  Several  indirect  sources  have  been 
investigated  and  have  been  shown  to  be  minor  sources.  The  exception  to 
this  is  the  quenching  of  N2(A)  by  atomic  oxygen  (Equation  [1.25])  as 
discussed  above. 

Other  proposed  reactions  to  produce  0(^S)  include  the  dissociation  of 
O2  by  EUV  photons  created  by  the  impact  of  electrons  on  N2  (Zipf,  1973). 
However,  laboratory  studies  of  Zipf  and  McLaughlin  (1978)  demonstrated 
that  the  EUV  photons  are  resonantly  trapped  by  N2,  causing  it  to 
predissociate  to  atomic  nitrogen.  Also,  Park  et  al.  (1977)  found  the  auroral 
intensity  between  950A  and  11 00 A  was  inadequate  for  the  0(^S)  production. 
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Yau  and  Shepherd  (1979)  and  Solheim  and  Llewellyn  (1979)  pointed 
out  that  the  electron  impact  cross  section  for  O2  to  the  C^Au,  and  c^Zu 

states  are  large  (Watson  et  al.,  1967;  Konishi  et  al.,  1970;  Vallance  Jones, 
1974)  but  that  the  Herzberg  bands  are  not  enhanced  in  the  aurora.  Yau  and 
Shepherd  (1979)  suggested  that  collisional  excitation  0(lS)  by  these  states 
may  be  important.  Until  the  mechanisms  for  the  quenching  of  0(lS)  are 
understood,  the  importance  of  this  reaction  will  not  be  quantified.  The 
source  may  not  be  needed  to  explain  the  observed  profiles  of  5577 A  emission, 
or  it  may  just  alter  the  efficiency  required  for  Reaction  [1.25], 

As  can  be  seen,  the  complete  answer  for  any  one  of  the  three 

o 

components  of  the  5577 A  emission  is  not  well  known.  As  new  information 
about  any  of  the  reactions  becomes  available,  the  emission  has  to  be 
remodeled.  This  is  compounded  since  a  change  in  the  accepted  chemistry 
in  one  region  can  affect  the  accepted  chemistry  in  the  others.  Any  changes 
in  one  region  must  be  consistent  with  the  chemistry  in  the  others. 

1.2.3  Measurement  Techniques 

Various  instruments  and  platforms  have  been  used  to  study  airglow 
emissions.  Photometers  provide  measurements  of  the  relative  strength  of 
atomic  and  molecular  emissions,  which  can  be  related  to  the  number 
densities  of  the  source  species.  A  spectrometer  can  give  information  on  the 
intensities  of  vibrational-rotational  bands.  The  relative  intensities  of  the 
lines  in  a  band  can  be  related  to  the  rotational  temperature  of  the  emitter, 
which  if  the  emitter  is  a  meta-stable  species  may  be  directly  indicative  of  the 
kinetic  temperature.  Wind  velocities  can  be  determined  from  line  position, 
and  temperatures  may  be  derived  by  the  analysis  line  width  with  an 
interferometer.  All  of  these  instruments  have  been  used  in  studies  from 
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ground  based  observatories,  as  rocket  borne  instruments  and  on  numerous 
satellites.  Measurements  of  the  5577A  "auroral  "  green  line  from  these 
various  instruments  provide  a  large  amount  of  the  information  on  the 
lower  thermosphere.  This  section  will  review  the  history  of  observations  of 
the  5577A  emission  by  these  different  means. 

Observations  of  the  5577k  airglow  emission  from  the  ground 
naturally  extends  from  auroral  observations.  Early  history  of  the  aurora 
and  visible  observations  of  emissions  is  given  in  the  monograph  by  Eather 
(1980).  Estimates  of  the  auroral  altitude  were  given  by  Cavendish  in  1784 
and  spectral  observations  were  carried  out  by  Angstrom  in  1868,  as 
mentioned  earlier.  In  1910  Stomer  made  the  first  studies  of  the  aurora 
using  the  photographic  triangulation  technique  (Vallance  Jones,  1974). 
Visual  reports  of  the  aurora  were  used  to  augment  mapping  of  the  auroral 
zone  as  late  as  1957-1958  (Millman,  1968). 

Ground  based  studies  of  the  airglow  cannot  be  used  to  determine 
altitude  variations  but  can  be  used  in  studies  of  the  long  term  changes  and 
morphology  of  the  greenline.  In  studying  the  green  line  from  the  ground, 
the  strong  emission  layer  at  about  95  km  hampers  observations  of  the  upper 
layer  in  the  F-region.  Photometric  studies  from  the  ground  include  Smith 
and  Steiger  (1968),  Wiens  and  Weill  (1973),  Fukuyama  (1976;  1977a;  and 
1977b),  Misaua  and  Takeuchi  (1978),  Freund  and  Jacka  (1979),  Armstrong 
(1982),  Takahashi  et  al.  (1977,  1984,  1985),  and  Saito  and  Kiyama  (1988). 
Photographic  and  imaging  studies  have  been  carried  out  by  Moreels  and 
Herse  (1977),  Armstrong  (1982)  and  Taylor  al.  (1987).  Freund  and  Jacka 
(1979),  Cogger  et  al.  (1985),  Weins  et  al.  (1988),  and  Lloyd  et  al.  (1990)  have 
also  performed  ground  based  interferometer  studies  of  the  5577k  emission. 
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Rockets  provide  the  opportunity  to  gain  high  resolution  data  of 
emission  profiles  with  height.  They  also  make  it  possible  to  study  the 
airglow  emissions  in  the  daytime  and  at  twilight.  Rockets  have  the 
advantage  over  ground  based  studies  of  being  able  to  provide  in  situ 
measurements  of  the  atmosphere  along  with  the  airglow  information. 

In  spite  of  the  short  measurement  time,  rockets  have  provided  many 
observations  of  the  green  line  emission.  With  photometers,  rockets  can  be 
used  to  determine  variations  of  the  emission  with  altitude  if  the  rocket  is 
traveling  vertically  and  the  emission  is  stable.  Rocket  based  instruments 
can  be  used  to  scan  the  emission  but  the  analysis  of  the  data  becomes  more 
difficult. 

Rockets  have  been  used  to  observe  the  5577A  green  line  emission  in 
airglow  and  aurora.  Measurements  of  the  nightglow  have  been  reported  by 
Berg  et  al.  (1956),  Koomen  et  al  (1956),  Heppner  and  Meridith  (1958), 
Tousey  (1958),  Cooper  et  al.  (1960),  Packer  (1961),  Huruhata  et  al.  (1962), 
Tarasova  (1963),  O'Brien  et  al.  (1965),  Greer  and  Best  (1967),  Baker  and 
Waddoups  (1967,  1968),  Gulledge  et  al.  (1968),  Danderkar  and  Turtle  (1971), 
Offermann  and  Drescher  (1973),  Kulkami  (1967),  Thomas,  L.  et  al.  (1979), 
Thomas,  R.  J.  (1981),  Lopez-Moreno  et  al.  (1982),  Greer  et  al.  (1986)  and  Kita 
et  al.  (1988).  Twilight  and  dayglow  measurements  of  the  01  5577 A  have 
been  presented  by  Silverman  al.  (1964),  Wallace  and  Nidey  (1964),  Wallace 
and  McElroy  (1966),  Lloyd  et  al.  (1968),  Dandekar  (1969),  Feldman,  Doering 
and  Zipf  (1971),  Schaffer  et  al.  (1972)  and  Hays  and  Sharp  (1973). 
Coordinated  satellite  and  rockets  measurements  for  aurora  were  reported 
in  a  series  of  papers  by  Rees  et  al.  (1977),  Sharp  et  al.  (1979),  and  Sharp  and 
Torr  (1979).  Some  characteristics  of  these  rocket  measurements  are 
summarized  in  Table  1.2. 


Long.  Alt  I  Emission  VER 
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Satellite  observations  of  airglow  can  provide  data  on  large  scale  global 
emission  features.  Satellites  have  the  capability  to  provide  observations 
equivalent  to  thousands  of  rocket  launches  and  give  the  spatial  coverage  of 
hundreds  of  ground  observatories.  Due  to  atmospheric  drag  which  would 
shorten  their  lifetime,  satellites  are  not  able  to  provide  the  in-situ 
measurements  of  rockets.  Information  can  be  obtained  about  atmospheric 
conditions  above  the  satellite. 

Satellites  have  provided  measurements  of  the  airglow  emissions 
since  the  late  1960's  (O'Brien,  1967;  Warnecke  et  ai,  1969;  Deehr  and 
Egeland,  1972;)  as  cited  by  Vallance  Jones  0974).  Reed  and  Chandra  (1975) 
presented  measurements  from  the  Orbiting  Geophysical  Observatory-4 
(OGO-4)  satellite.  OGO-6  observations  were  presented  by  Thomas  and 
Donahue  (1972),  Donahue  et  al.  (1973,  1974),  Donahue  (1975)  and  Wasser 
and  Donahue  (1979).  The  International  Satellites  for  Ionospheric  Studies-2 
(ISIS-2)  satellite  used  simple  photometers  to  produce  scans  of  the  aurora 
and  airglow  (Anger  et  ai,  1973;  Shepard  et  al.,  1973;  Cogger  and  Anger, 
1973;  Cogger  and  Murphee,  1980;  Cogger  et  ai,  1981;  Elphinstone,  1986). 
The  Atmosphere  Explorer  (AE)  satellite  program  studied  aurora  and 
airglow  emissions  for  a  period  from  near  solar  minimum  to  solar 
maximum.  This  satellite  program  is  discussed  in  detail  in  Chapter  2. 

The  two  satellite  Dynamics  Explorer  (DE)  program  used  a  Fabry- 
Perot  interferometer  (Hays  et  al.,  1981)  to  measure  airglow  and  auroral 
emission  (Killeen  and  Hays,  1983).  The  High  Resolution  Doppler  Imager 
(HRDI)  on  the  recently  launched  Upper  Atmosphere  Research  Satellite 
(UARS)  carries  a  combination  instrument  that  contains  a  photometer  and  a 
triple  etalon  Fabry-Perot  that  can  measure  5577A  emissions  (described  by 
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Satellite 

Launch  Date 

Instrument 

Start  of  Data 

End  of  Data 

Orbit 

Inclination 

OGO-4 

28JUL67 

Photometer 

30AUG67 

10JAN68 

86 

OGO-6 

5JUL69 

Photometer 

9JUN69 

24JIJL70 

-82 

ISIS-2 

1APR71 

Photometer 

31DEC71 

AE-C 

16DEC73 

Photometer 

16DEC73 

11DEC78 

-58 

AE-D 

60CT75 

Photometer 

60CT75 

-90 

AE-E 

20NOV75 

Photometer 

21NOV75 

07JUL81 

-20 

DE-2 

3AUG81 

4AUG81 

18FEB83 

90 

Intercosmos 

Bulgaria 

1300 

7AUG81 

Photometer 

81.9 

UARS 

12SEP91 

liiiii 

I^HI 

HHI 

-58 

Table  1.3  Satellite  measxirements  of  5577A  emission. 
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Bucholtz,  1991).  Table  1.3  summarizes  some  characteristics  of  these 
programs. 

1.3  Statement  of  Thesis  Problem 

The  major  purpose  of  this  thesis  is  to  provide  the  details  of  the 
photochemical  model  developed  for  the  upper  layer  of  the  5577 A  airglow 
emission.  In  presenting  the  background  and  history  of  the  emission,  the 
many  potential  source  and  loss  mechanisms  have  been  discussed.  The 
relative  importance  of  some  of  these  reactions  has  not  been  conclusively 
determined.  This  research  provides  information  on  which  reactions  are 
required  to  reproduce  the  observations  of  the  5577 A  dayglow  made  by  the 
Atmosphere  Explorer  satellites. 

Two  reactions  have  been  shown  to  be  major  sources  of  the  emission 
in  the  upper  layer.  The  photoelectron  impact  source  [Reaction  1.19]  is 
important  at  the  higher  altitudes  of  the  layer,  and  energy  transfer  from 
Equation  [1.25]  becomes  more  important  at  the  lower  part  of  the  layer.  Both 
of  these  reactions  have  parameters  that  have  not  been  conclusively 
established. 

The  impact  cross  section  for  Reaction  [1.19]  has  as  much  as  a  factor 
of  two  difference  in  magnitude  between  theory  and  laboratory 
measurements.  The  magnitude  of  the  cross  section  that  best  reproduces 
the  observed  data  will  be  determined.  Also,  Reaction  [1.25]  has  both  a  total 
rate  coefficient  and  the  branching  ratio  for  the  production  of  0(^S)  that  are 
not  well  known.  Determinations  of  the  branching  ratio  in  auroral  studies 
give  relatively  low  values  of  0.10  to  0.30,  while  laboratory  determinations 
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place  the  value  at  0.75.  The  model  is  used  to  determine  those  parameters 

that  best  reproduce  the  observed  values  of  the  emission. 

0 

The  model  is  then  used  to  study  the  5577A  emission  in  the  twilight. 
The  model  has  been  run  for  the  conditions  in  the  terminator  and  the  results 
compared  to  the  observations  of  the  Atmosphere  Explorer  satellites.  The 
output  of  the  model  is  integrated  to  provide  models  of  the  integrated 
brightness  which  can  then  be  compared  to  those  observed  by  the  satellite. 


CHAPTER  n 


ATMOSPHERE  EXPLORER  SATELLITE  OBSERVATIONS 


2.1  Overview 

The  Atmosphere  Explorer  satellites  (AE-C,  -D,  and  -E)  were  used  to 
provide  measurements  both  of  the  5577A  airglow  and  of  atmospheric 
parameters  to  be  used  for  the  modelling  of  the  emissions.  This  chapter 
begins  with  a  brief  history  and  overview  of  the  satellite  program.  Next,  the 
Visible  Airglow  Experiment  (VAE)  is  described  and  data  handling  for  the 
instrument  is  presented.  A  general  discussion  of  the  inversion  methods  to 
derive  Volume  Emission  Rate  (VER)  profiles  from  the  brightness 
observations  is  also  included.  The  chapter  concludes  with  a  discussion  of 
the  AE  instruments  used  to  make  in  situ  measurements  of  neutral  and 
ionized  atmospheric  composition. 

The  Atmosphere  Explorer  satellite  program  has  a  rich  history.  The 
satellite  program  has  been  summarized  by  Burgess  and  Torr  (1987).  This 
review  places  the  emphasis  on  the  later  missions,  which  were  the  first 
examples  of  problem-dedicated  missions.  In  this  case,  a  team  approach 
was  used  to  study  the  system  of  the  incoming  solar  flux,  the  neutral 
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atmosphere  and  the  resulting  ionization.  The  Atmosphere  Explorer  team 
consisted  of  theorists  and  investigators  making  the  correlated 
measurements  required  to  understand  this  complex  region  of  the 
atmosphere. 

In  the  mid  1960's,  NASA  launched  Explorer  17  and  Explorer  32,  also 
called  AEl  and  AE2.  The  lower  thermosphere  was  largely  unknown  and 
the  need  to  study  this  region  specifically  led  to  the  Atmosphere  Explorer 
Program.  The  program  consisted  of  three  satellites  designated  AES,  AE4, 
and  AES  (or  AE-C,  -D,  and  -E). 

The  purpose  of  the  Atmosphere  Explorer  program  was  to  provide 
simultaneous  measurements  of  the  parameters  needed  to  study  the 
physical  processes  that  govern  the  lower  thermosphere  and  ionosphere. 
Three  satellites  were  required  to  provide  the  measurements,  each  with 
different  orbital  characteristics.  The  spacecraft  were  to  occupy  elliptical 
orbits  during  the  first  year  with  perigee  as  low  as  135  km  and  apogee  about 
4000  km,  and  subsequently  the  orbits  were  to  be  circularized.  Each  satellite 
had  a  propulsion  system  that  would  maintain  the  orbit  both  during  the 
elliptical  phase  and  during  the  circular  phase  where  the  altitude  could  be 
adjusted.  Table  2.1  is  a  listing  of  the  history  of  the  orbital  altitudes  for  the 
satellites.  Even  during  the  circular  phase,  the  altitude  of  the  satellites 
could  vary  by  a  few  kilometers  during  a  single  orbit. 

The  orbital  inclinations  of  the  three  satelites  were  chosen  to  provide 
different  spatial  and  temporal  coverages.  The  inclination  of  AE-C  was  to  be 
at  about  68,  degrees  which  allowed  measurements  of  low  and  midlatitudes 
for  a  range  of  local  times.  AE-D  was  a  polar  orbiting  satellite  with  an 
inclination  near  90  degrees.  Where  AE-C  could  study  the  thermosphere 
with  good  local  time  coverage,  AE-D's  high  inclination  allowed  a  complete 
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Table  2.1  Orbital  history  of  the  Atmosphere  Explorer  satellites. 
Perigee  altitudes  are  approximate.  Elliptical  orbital 
eccentricities  are  between  0.05  and  0.25,  Apogees  from  400 
to  2200  km. 
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latitudinal  coverage.  AE-E  was  placed  in  an  orbit  of  about  20  degrees, 
which  allowed  continuous  measurements  at  low  latitudes  to  be  used  in 
correlative  studies  with  the  incoherent  backscatter  station  at  Arecibo, 
Puerto  Rico.  The  effects  of  perigee  drift  for  each  satellite  as  a  function  of 
local  time  and  latitude  is  illustrated  in  Figure  2.1. 

The  satellites  each  carried  a  variety  of  instruments  for  remote 
sensing  and  in-situ  measurements  of  the  atmosphere.  These  instruments 
included  various  spectrometers,  photometers  and  ion  probes.  Listed  in 
Table  2.2  is  the  various  instruments  used  on  the  different  satellites.  These 
instruments  are  reviewed  in  a  special  issue  of  Radio  Science  (Volume  8, 
Number  4,  1973).  The  satellites  were  initially  launched  with  a  minimiun 
design  lifetime  of  one  year.  AE-C  and  AE-E  exceeded  this  by  many  years, 
but  AE-D  failed  totally  after  a  few  months  and  never  reached  the  circular 
phase  of  the  mission.  Table  2,3  is  a  summary  of  the  orbital  characteristics 
for  each  satellite. 

The  satellites  were  spin  stabilized  and  could  be  used  in  either  a 
spinning  or  despun  mode.  In  the  spinning  mode,  the  satellite  rotated  with 
a  period  of  15  seconds,  with  the  spin  vector  either  parallel  or  antiparallel  to 
the  orbit  normal.  This  mode  enabled  the  instruments  to  sweep  through  all 
angles  of  attack  within  each  spin  period.  In  the  despun  mode,  the  satellite 
would  rotate  once  per  orbit  in  order  to  maintain  the  same  attitude  for  the 
instruments  relative  to  the  local  horizon. 

A  major  achievement  for  the  Atmosphere  Explorer  program  was  the 
handling  of  the  data  stream  from  the  spacecraft.  The  Data  was  received 
from  the  spacecraft  and  processed  in  near  real  time  to  allow  investigators  to 
modify  the  mission  plans  quickly  in  order  to  take  advantage  of  any 
observations  that  were  found  to  address  particular  aeronomic  problems. 
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Fiifure  2.1  Perigee  dnfls  for  the  Atmosphere  Explorer  satellites. 
Altitudes  ranged  from  150  to  4000  km.  (From  Dalaamo  et 
al.,  1973)  ® 


1  Instrument 

II 

EUVS 

Solar  EUV  Filter 
Photometer 

ESUM 

UV  Nitric  Oxide 

UVNO 

Visible  Airglow 
Photometer 

VAE 

Open  Source  Neutral 
Mass  Spectrometer 

OSS 

Closed  Source 
Neutral  Mass 
Spectrometer 

NACE 

1 

Neutral  Atmospheric 
Temperature 
Experiment 

NATE 

Atmospheric  Density 
Accelerometer 

MESA 

Planar  Ion  Trap 

RPA 

L 

Cylindrical 
Electrostatic  Probe 

CEP 

L 

Magnetic  Ion  Mass 
Spectrometer 

MIMS 

Positive  Ion  Mass 
Spectrometer 

BIMS 

Low  Energy  Electron 
Experiment 

LEE 

1 

Photoelectron 

Spectrometer 

PES 

Capacitance 

Manometer 

Temperature  Alarm  | 

IB 

1 

Satellite 

Detector 

1  Parameters 

C,D.E 

Channel  Electron 
Multipliers 

C.D.E 

Spiraltron  Electron 
Multitpiers  and 
EUV  diodes 

llllliHHI 

C.D 

Photomultiplier 

Tubes 

IIB3SBS3I 

C.D.E 

Photomultiplier 

Tubes 

6300^77.42783371, 
5200,7319  to  733QA 

C.D.E 

Electron  Multiplier 

lto46AMU 

mm 

Electron  Multiplier 

lto46AMU 

C,D,E 

Electron  Multiplier 

Tg.N2.Vw 

C.D.E 

Accelerometer 

Neutral  Density 

C.D,E 

Electrometer 

C,D.E 

Electrometer 

Te.Ne.  Ni,  Mi 

C,D 

Electron  Multiplier 

lto64AMU 

1 

C.E 

Electrometer 

0.5  to  72  AMU 

mu 

Spiraltron  Electron 
Multipliers 

0.2to25keV 

C,D.E 

Johnson  Electron 
Multiplier 

Photoelectron 

Spectra 

C,D,E 

Electrometer 

Pressure 

C,D,E 

Diaphragm 

Pressure 

mu 

m^mmi 

Table  2J2  Instrument  package  lists  for  the  Atmosphere  Explorer 
satellites. 
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Orbital  1 
Inclination  I 

8 

19.9  1 

Circular 

Phase 

DEC74 

1 

N6V76 

Re-entry 

date 

12DEC78 

12MAR76 

10JUN81 

Launch 

Date 

13DEC73 

60CT75 

20NOV75 

5 

U 

Q 

U 

CO 

CO 
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The  data  from  the  Atmosphere  Explorer  program  is  available  in  different 
formats.  The  complete  database  known  as  the  Unified  Abstract  database 
contains  15  second  averages  of  the  data  from  all  of  the  various  instrument 
and  other  important  geophysical  parameters.  This  database  is  available  on 
tapes  provided  by  the  National  Space  Science  Data  Center  (NSSDC)  or 
through  the  Space  Physics  Analysis  Network  (SPAN).  Higher  resolution 
databases  are  available  for  each  of  the  instruments  on  magnetic  tape  and/or 
in  some  cases  on  microfilm  plots.  These  databases  and  their  applicability  to 
the  airglow  model  will  be  discussed  later. 

The  Atmosphere  Explorer  satellites  were  very  successful.  Most  of  the 
instruments  provided  data  for  much  longer  than  the  expected  lifetime  of  the 
satellites,  and  due  to  the  built-in  redundancy  failures  could  be  replaced  by 
other  instruments.  The  knowledge  of  the  structure,  chemistry  and 
dynamics  of  the  upper  atmosphere  was  greatly  advanced  by  this  program 
and  provided  the  opportunity  for  correlative  observations  by  hundreds  of 
rockets  flights  and  several  other  satellite  programs.  The  data  from  the 
Atmosphere  Explorer  satellites  is  still  used  to  study  the  upper  atmosphere 
and  is  still  being  published.  The  AE  measurements  have  also  been 
employed  by  semi-empirical  models  like  the  Mass  Spectrometer  and 
Incoherent  Scatter  (MSIS)  model. 

2.2  The  Visible  Airglow  Experiment 


2.2.1  Instrument  Description 

The  Visible  Airglow  Experiment  (VAE)  was  flown  on  all  three  of  the 
Atmosphere  Explorer  satellites  and  has  been  described  by  Hays  et  al.  (1973) 
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and  by  Burgess  and  Torr  (1987).  The  instrument  consisted  of  a  simple  two 
channel  filter  wheel  photometer  similar  to  that  flown  on  the  ISIS-B 
program  (Burgess  and  Torr,  1987).  It  was  designed  to  monitor  dayglow, 
nightglow  and  auroral  emissions  in  the  upper  atmosphere.  The  results  of 
previous  studies  using  the  VAE  instrument  have  been  summarized  by 
Hays  et  al.  (1988). 

In  the  VAE  instrument,  shown  in  Figure  2.2,  the  two  channels  were 
perpendicular  to  each  other  and  oriented  in  the  satellite  such  that  in  the 
despun  mode  channel  1  was  pointed  back  along  the  satellite  track  and 
channel  two  looked  either  up  (radially  out  from  the  Earth)  or  down 
(radially  in)  depending  on  whether  the  satellite  was  in  the  normal  or 
inverted  position,  respectively.  When  the  satellite  was  in  a  spinning  mode, 
the  two  channels  turned  in  the  plane  of  the  orbit,  with  channel  one  always 
leading  channel  two  by  90  degrees.  When  the  satellite  was  in  the  normal 
position  and  spinning  it  would  have  a  "skidding"  motion,  while  the  inverted 
position  would  give  a  rolling  motion.  The  orientation  of  the  instrument  and 
the  motion  is  shown  in  Figure  2.3. 

The  two  photometer  channels  of  the  VAE  instrument  were 
essentially  the  same  with  a  combination  of  single-objective  lens  and  field 
stop  used  to  define  two  different  fields  of  view  for  the  detectors.  The  channel 
1  photometer  had  a  half  angle  cone  of  0.75®  and  was  used  to  measure 
dayglow  and  nightglow  horizons  and  to  observe  auroral  features  with 
strong  spatial  gradients.  Channel  2  was  a  high-sensitivity  instrument  with 
a  3°  half  angle  cone  that  could  be  used  to  measure  nightglow,  dayglow  above 
the  satellite  and  other  weak  emissions  but  without  information  on  spatial 
gradients.  The  integration  periods  for  the  two  channels  were  matched  to 
the  spin  rate  of  the  satellites  so  that  the  field  of  view  was  roughly  equal  to 
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Figure  2.2  Schematic  diagram  of  the  Visible  Airglow 
Experiment  instrument  used  on  the  Atmosphere  Explorer 
satellites. 
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Channel  2 


Satellite  in  Normal  Mode 

y 


Channel  2 


Satellite  in  Inverted  Mode 


Fi|nu^  2.3  Illustration  of  VAE  photometer  orientation  and  spin 
directions  for  normal  and  inverted  modes. 
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the  angular  interval  of  the  measurement.  Channel  1  was  integrated  for 
0.031  second  to  give  a  sensitivity  of  about  20  Rayleighs/count/integration. 
For  Channel  2  the  resolution  was  about  0.3  Rayleighs/count/integration 
with  an  integration  time  of  0.125  seconds. 

The  two  channels  of  the  VAE  instrument  shared  a  common  filter 
wheel.  The  filter  wheel  contained  six  interference  filters  with  spectral 

o 

resolution  of  15  to  30A,  a  dark  count  position  and  a  phosphor  calibration 
source  position.  The  different  satellites  carried  filters  for  different 
emissions  as  well  as  some  common  emission  filters.  The  emissions  and 
the  source  on  each  satellite  are  listed  in  Table  2.4.  The  filter  wheels  could 
be  operated  in  either  a  fixed  mode,  or  in  a  stepping  mode  that  would  change 
the  filter  of  the  channels  every  4,  8  or  16  seconds.  Because  of  the  large 
spectral  range  of  the  filters,  photon  contamination  from  sources  other  than 
those  of  Table  2.4  could  be  detected  by  the  photometers. 

The  photometers  in  the  VAE  instrument  were  protected  from  high 
dynode  currents,  which  is  the  primary  cause  of  failure,  by  two  methods. 
The  first  was  a  1/100  attenuator  that  acts  to  extend  the  dynamic  range  by  2 
orders  of  magnitude  when  the  count  rate  exceeds  an  upper  threshold.  The 
second  system  was  a  "squint"  mechanism.  This  system  used  a  secondary 
counting  system  with  a  1.0  msec  integration  period  that  adds  a  bias  voltage 
to  the  cathode  when  the  count  rate  exceeds  an  upper  threshold.  The  back 
bias  would  decay  slowly  to  normal  and  the  instrument  would  "squint" 
repeatedly  until  the  count  rate  was  below  a  critical  value.  This  mechanism 
was  used  to  prevent  damage  arising  from  the  passage  of  the  sun  through 
the  instrument  field  of  view  or  from  the  bright  dayside  of  the  Earth.  With 
these  systems,  the  VAE  instrument  could  measure  a  dark  feature  within 
120.0  msec  after  directly  viewing  the  sun. 
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Emission  (A) 

Source 

C 

6300 

Ol(’P-‘D) 

5577 

Ol('D-‘S) 

7319-7330 

OIl(*D-*P) 

5200 

NI("S-*D) 

3371 

Njl(B“n,.v  =  0-C"n„,v  =  0) 

4278 

D 

6300 

Ol(®P-'D) 

5577 

0I(‘D-‘S) 

7319-7330 

OIl{*D-*P) 

5200 

NI(*S-*D) 

4278 

Nji(x*z;.v  =  i-B*i:,v=o) 

4861 

HI(n  =  2-4) 

E 

6300 

0I(’P-‘D) 

5577 

Ol(‘D-'S) 

7319-7330 

0II(*D-*P) 

5200 

NI('S-*D) 

2802 

MgIl(*S-*P) 

6563 

HI(n  =  2-3) 

Table  2.4  Emissions  measured  by  the  Visible  Airglow 
Experiment  on  the  Atmosphere  Explorer  satellites  with 
the  emission  sources. 
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As  noted  before,  there  can  be  sources  of  photons  that  can  pass  through  the 
interference  filters  other  than  those  listed  in  Table  2.4.  In  addition  to  the 
photons  from  other  emissions  that  fall  within  the  bandpass  of  the  filters, 
there  may  be  other  non-airglow  optical  contamination.  The  problem  of 
vehicle  glow,  from  spacecraft  atmosphere  interaction,  can  be  severe  near 
perigee,  especially  during  elliptical  phases,  and  is  most  pronounced  toward 
the  red  end  of  the  spectrum.  This  source  of  contaminationhas  been 
discussed  for  the  Atmosphere  Explorer  satellites  by  Yee  and  Abreu  (1983). 
The  contribution  from  galactic  and  zodical  background  or  from  a  planet  in 
the  field  of  view  can  be  another  source  of  contamination.  The  third  source 
of  spurious  signal  is  from  scattering  by  the  optical  baffle  systems  of  the 
telescopes.  This  is  a  pai-ticular  problem  when  observing  near  the  bright 
limb  of  the  Earth.  Curves  showing  valid  regions  of  dayglow  limb 
observations  were  given  by  Fesen  (1981).  The  influence  of  these 
contamination  sources  on  this  work  will  be  discussed  later. 

2.2.2  Analysis  of  Visible  Airglow  Data 

The  data  gathered  by  the  Visible  Airglow  Experiment  is  not  a  direct 
measure  of  the  volume  emission  rate  at  a  point  in  the  atmosphere.  The 
recorded  quantity  is  the  integrated  surface  brightness  along  the  path  of  the 
observation  in  Rayleighs.  A  Rayleigh  is  a  unit  of  surface  brightness  equal 
to  1x10®  photons  /  (cm2  [column]  sec).  For  the  satellite  observations,  the 
column  is  along  the  line  of  sight  of  the  instrument.  To  turn  the  measured 
surface  brightness  into  a  volume  emission  rate  (VER)  the  data  is  inverted 
by  different  techniques,  depending  on  the  type  of  orbit.  Inversion  methods 
recover  a  set  of  physical  quantities  from  a  set  of  integral  measurements. 
The  term  inversion  refers  to  the  inverse  transform  which  is  sought  to  undo 
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the  effects  of  the  measurement  process.  In  this  section,  some  of  the 
techniques  that  have  been  used  to  recover  volume  emission  rates  from  the 
measurements  made  by  the  Atmosphere  Explorer  satellites  will  be 
discussed.  The  list  will  not  be  exhaustive  but  will  serve  to  present  some  of 
the  various  classes  of  the  techniques. 

Since  the  Atmosphere  Explorer  satellites  employed  several  different 
types  of  orbits  (  spinning  or  de-spun  and  elliptical  or  circular),  there  are 
several  "inversion"  processes  that  may  be  used  to  recover  volume  emission 
rates  from  the  brightnesses.  As  ordered  here  these  techniques  fall  into 
three  basic  categories.  The  first  group  is  not  actual  inversion  techniques 
but  uses  different  mathematical  methods  to  determine  the  volume  emission 
rate  profiles  from  brightnesses.  The  other  two  groups  are  analytical 
inversion  methods  and  linear  inversion  methods.  Examples  of  methods  in 
each  group  these  will  be  discussed  in  terms  of  the  type  of  orbit  and  geometry 
to  recover  the  desired  data. 

One  of  the  simplest  means  to  recover  the  volume  emission  rate  from 
elliptical,  de-spun  orbits  is  to  use  data  from  channel  two  looking  radially 
out  from  the  Earth  in  a  technique  similar  to  that  used  for  rocket  data.  This 
technique  uses  a  series  of  brightnesses  as  the  satellite  crosses  various 
layers  of  the  atmosphere  (see  Figure  2.4).  The  brightness  profile  generated 
this  way  can  be  differentiated  to  produce  volume  emission  rate  profiles. 
This  technique  requires  the  emission  to  be  homogeneous  and  spherically 
distributed.  Another  restriction  is  that  the  background  does  not  change 
appreciably  between  two  successive  measurements.  These  conditions  can 
usually  be  met  by  restricting  the  data  to  midlatitudes  and  away  from  the 
terminators.  This  technique  may  be  used  on  spinning  elliptical  orbits  by 


Figure  2.4  Geometry  of  differential  "inversion"  used  in  despun 
normal  modes. 
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using  measurements  when  the  photometer  is  pointing  radially  out  within 
some  small  angle. 

A  technique  that  uses  circular  spinning  orbits  for  the  brightness 
measurements  but  is  not  an  inversion  in  the  true  sense  is  the  chord 
subtraction  technique  (Abreu  et  aL,  1983;  Yee  and  Abreu,  1987).  The  basic 
geometry  of  this  technique  is  shown  in  Figure  2.5.  When  the  satellite  is 
spinning,  measurements  are  made  with  the  photometer  looking  in  the 
same  direction  for  different  positions  in  the  orbit.  These  measurements  are 
subtracted  and  divided  by  the  path  length  to  determine  the  VER  along  the 
chord.  The  same  calculation  is  carried  out  for  another  chord  and  the  two 
are  averaged  to  give  the  VER  at  the  point  desired.  This  technique  corrects 
for  the  background  since  the  difference  in  the  measurements  is  just  the 
emission  along  the  chord  bounded  the  satellite  orbit  if  the  background 
changes  slowly  compared  to  the  time  between  the  observations. 

Another  technique  which  is  similar  to  the  chord  subtraction  method 
is  the  "onion  peeling"  method  (Fesen,  1981;  Fesen  and  Hays,  1982).  It  is  a 
two  dimensional  inversion  that  finds  a  non-iterative  solution  to  the  limb 
scan  data.  Figure  2.6  shows  the  inversion  geometry  and  parameters.  The 
method  uses  the  nature  of  the  limb  scan  and  orbit  in  the  algebraic  domain 
by  starting  at  the  top  and  working  down  through  the  atmosphere.  By 
starting  at  the  satellite,  the  contribution  to  the  emission  from  the  first  layer 
below  the  satellite  is  found.  This  is  then  used  to  determine  the  contribution 
from  the  second  layer,  and  knowledge  of  both  contributions  is  then  applied 
tot  the  next  layer.  This  procedure  is  continued  down  through  the 
atmosphere,  but  to  prevent  oscillations  from  occurring,  a  smoothness 
constraint  is  applied. 
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The  majority  of  the  inversion  techniques  are  based  on  brightness 
measurements  from  circular  spinning  orbits.  The  general  geometry  for 
this  type  of  observation  is  shown  in  Figure  2.7.  The  measured  brightness  is 
the  integral  of  the  volume  emission  rates  along  the  line  of  sight.  The  basic 
equation  for  this  is 


t2.1] 

where  B  is  the  measured  brightness  along  a  path  with  a  tangent  height  rt 
and  q  is  the  volume  emission  rate  along  the  pathlength  s.  Using  the 
geometry  shown  in  Figure  2.8  for  a  single  path,  equation  [2.1]  can  be  re¬ 
written  as 


B(r.)  =  dr.  [2.2] 

'■(r'-r?)! 

In  this  case,  rg  is  the  satellite  height  and  the  volume  emission  rate  is 
expressed  as  a  function  of  height.  Equation  [2.2]  is  the  basic  equation  used 
in  most  of  the  inversion  processes. 

The  Abel  inversion  (Roble  and  Hays,  1972)  was  the  first  inversion 
applied  to  VAE  data  from  spinning  orbits.  The  Abel  inversion  is  loosely  a 
tomographic  inversion  since  it  can  be  considered  as  a  one  dimensional  case 
of  the  Cormack  inversion  where  only  variation  in  altitude  is  allowed.  The 
geometry  for  the  Abel  inversion  is  shown  in  Figure  2.9.  Equation  [2.2]  is 
essentially  one  of  a  pair  of  equations  that  make  up  the  Abel  Transform. 
This  equation  has  an  analytical  solution  which  may  be  expressed  in  first  or 
second  derivative  form.  The  first  derivative  form  of  the  inverse  is. 


Figure  2.8  Geometry  of  a  single  surface  brightness 
measurement  from  satellite  at  height  rg  and  photometer 
zenith  angle  z. 
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B(r,) 

dr, 

dr^  , 


[2.3] 


from  r  to  the  satellite  altitude.  The  second  derivative  form  of  the  inverse 
equation  is, 


dr. 


dr. 


B(r,) 


dr^  , 


[2.4] 


over  the  same  limits  of  integration. 

The  Cormack  inversion  is  a  two  dimensional  inversion  that  has 
many  applications  in  different  fields  (Cormack,  1963;  1964;  197u).  It  was 
developed  for  use  in  the  medical  fields  and  is  also  known  as  tomography. 
In  1963,  A.  M.  Cormack  solved  the  problem  of  obtaining  two  dimensional 
density  function  from  multiple  x-ray  measurements  by  estimating  the 
function  from  knowledge  of  line  integrals.  The  problem  had  been  solved 
differently  much  earlier  by  Radon  (1917)  but  was  forgotten  since  it  would  be 
an  academic  interest  at  the  time  (Berry  and  Gibbs,  1970).  The  Abel 
inversion  is  a  special  case  of  the  Cormack  inversion.  In  the  Abel  case  the 
horizontal  gradients  are  assumed  to  be  small.  For  the  Cormack  inversion 
the  spatial  dependance  is  allowed  to  vary.  Figure  2.10  shows  the  geometry 
for  the  Cormack  inversion. 

The  Cormack  equation,  which  is  similar  to  that  for  the  Abel 
inversion,  is 


le  Cormack  inversii 
all  g(r,6)  is  8pecifie( 
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F.(rJ  =  2j 


G.(r)T. 


I'i'j 

.  r  j 


(r^- 


dr  , 


[2.5] 


T„  =  Chebyshev  Polynomial 
F„  =  Fourier  Series  Expansion  of  Brightness 
G„  =  Fourier  Series  Expansion  of  Volume  Emission  Rate 
The  inverse  of  the  Cormack  equation  is, 


G„(r)  = 


[2.6] 


Discussion  of  the  application  of  this  technique  to  satellite  data  has  been 
given  by  Solomon  et  al.  (1984,  1985);  Solomon  (1987);  Abreu  et  al.  (1989).  The 
Cormack  inversion  is  particularly  useful  in  the  cases  where  there  is  a  high 
degree  of  spatial  variability  of  the  emission  as  in  aurora. 

Equation  [2.2]  can  itself  be  re-written  into  a  matrix  representation. 


G  =  A  F  . 


{2:n 


In  this  case,  the  observed  brightness  is  the  vector  G  and  the  corresponding 
volume  emission  rates  are  in  the  vector  F.  The  matrix  A  contains  the 
weighting  function  or  kernel  which  relates  the  two  vectors. 


A  = 


[2.8] 


In  the  case  of  limb  scans,  the  kernel  gives  the  contribution  of  each  layer  of 
the  atmosphere  to  the  integrated  brightness  due  to  the  geometry  of  the 
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observation.  Figure  2.11  is  an  example  of  the  weighting  functions  at  three 
tangent  heights  for  the  geometric  kernel. 

The  task  of  the  inversion  is  to  find  a  matrix  B  such  that 

F  =  G  ,  [2.9] 

where  B  can  be  identified  as  the  inverse  of  A,  if  A  is  a  square  matrix  or  can 
be  made  square.  Only  in  the  simplest  cases  can  the  matrix  be  inverted 
directly.  The  weighting  functions  must  be  sharp  and  not  have  significant 
overlap  to  be  useful  in  the  inversion  process.  In  typical  remote  sensing 
applications,  the  measurement  G  contains  some  errors  and  the  matrix  A 
is  nearly  singular.  The  matrix  becomes  singular  when  some  of  the 
weighting  functions  are  so  broad  that  there  is  significant  overlap  between 
tangent  heights.  If  this  occurs  the  rows  of  the  matrix  may  not  be  linearly 
independant.  Finding  F  in  this  case  may  lead  to  large  oscillations.  The 
mathematical  explanation  of  this  is  that  if  the  eigenvalues  of  A  are  small 
(the  matrix  is  nearly  singular)  then  the  eigenvalues  of  B  will  be  large;  this 
will  lead  to  instability  in  the  calculation  of  F  (Twomey,  1977;  Menke^  1989). 

By  using  various  linear  methods,  the  problem  of  solving  Equation 
[2.7]  in  less  than  the  ideal  cases  can  be  reduced.  In  remote  sensing 
applications  the  number  of  measurements  is  larger  than  the  number  of 
values  to  be  recovered.  This  will  cause  the  problem  to  be  over  determined 
but  it  can  be  made  solvable  by  linear  methods  by  forming  the  least  square 
solution  which  is. 


F  =  (A'A)’'a’  G  , 


[2.10] 
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Figure  2.11  Geometric  weighting  kernels  for  the  tangent  height 
measurements  used  by  the  Abel  inversion  method. 
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where  A*  is  the  transpose  of  the  matrix. 

Since  Equation  [2.10]  allows  for  the  solution  of  non-square  matrix 
equations  it  will  allow  for  the  solution  to  be  constrained  by  adding  more 
equations  to  the  system.  A  smoothness  constraint  may  be  applied  in  this 
manner.  The  solution  to  Equation  [2.7]  is  now, 

F  =  {A*A  +  G  .  [2.11] 

In  this  solution,  y  is  the  weight  applied  to  the  constraint  and  H  is  the 
matrix  constructed  by  the  differences  that  will  give  the  desired  level  of 
smoothness. 

Another  form  of  the  solution  for  Equation  [2.7]  is  one  that  allows  for 
the  system  to  be  constrained  to  a  priori  information  such  as  a  model  or 
climatology.  When  allowing  for  the  measurements  and  the  constraint  to 
have  errors  the  solution  to  the  system  becomes, 

F  =  (A*E-‘A  +  y/)''(A*E''  G  -t-  y//)  .  [2.12] 

In  this  solution,  y  is  the  weight  applied  to  the  constraint  and  E  is  the 
covariance  matrix  of  the  measurements.  /  is  the  covariance  matrix  of  the 
a  priori  information  and  H  is  the  a  priori  information  vector. 

A  final  linear  method  that  can  be  applied  to  the  measurements  of  the 
Atmosphere  Explorer  satellites  is  from  optimal  estimation  theory  or 
sequential  estimation  (Rodgers,  1976;  Houghton  ef  aZ. ,  1984;  Rodgers,  1990). 
Using  optimal  estimation,  the  solution  to  Equation  [2.7]  becomes, 

t  =  S.A’'(ylS.A^  +  S,)'‘(G-AF,) 


[2.13] 
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S  =  S,  -  S,A'^(AS,A'^  +  SJ'aS,  ,  [2.14] 

where  F  is  the  returned  volume  emission  rate  and  S  is  the  covariance  of 
the  returned  volume  emission  rate.  The  a  priori  estimate  for  the  volume 
emission  rate  and  its  covariance  is  given  by  Fq  and  respectively.  The 
error  covariance  for  the  measured  values  is  .  Since  these  equations  can 
be  solved  sequentially,  it  is  not  necessary  to  perform  a  matrix  inversion  or 
solve  simultaneous  linear  equations  to  retrieve  a  volume  emission  rate. 
Besides  the  savings  in  computer  time,  the  concept  of  sequential  estimation 
allows  the  use  of  continuity  along  an  orbit  to  improve  the  accuracy  of  the 
retrievals.The  application  of  this  method  has  been  discussed  by  Rodgers 
(1976)  and  Houghton  et  al.  (1984).  This  method  has  been  used  to  analyze 
Atmosphere  Explorer  data  by  Abreu  et  al.  (1991). 

There  are  numerous  techniques  to  recover  the  volume  emission  rate 
from  the  measured  surface  brightnesses.  The  choice  of  inversion  method  is 
primarily  determined  by  the  type  of  orbit  used  to  measure  the  data.  The 
methods  described  here  have  been  applied  to  data  taken  by  the  Visible 
Airglow  Experiment  data  on  the  Atmosphere  Explorer  satellites.  In  the 
next  section  of  this  chapter,  instruments  used  for  other  measurements 
from  the  Atmosphere  Explorer  satellites  needed  for  the  modelling  of  the 
0(^8)  emission  will  be  described. 

2.3  Other  Atmospheric  Parameters 

The  various  instruments  on  the  Atmosphere  Explorer  satellites  were 
able  to  provide  in-situ  measurements  of  several  important  atmospheric 
parameters.  The  neutral  and  ion  densities  were  measured  along  with  their 
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associated  temperatures.  In  most  cases  each  of  these  parameters  were 
determined  by  more  than  one  instrument.  In  this  section,  these 
parameters  are  broken  down  into  neutral  measurements  and  ion 
measurements.  Each  of  these  are  discussed  separately,  as  are  the 
individual  instruments  and  the  inter-relationship  between  the 
instruments. 

2.3.1  Neutral  Measurements 

The  neutral  atmospheric  parameters  that  are  important  to  the 
production  of  5577A  emission  are  the  molecular  nitrogen  density,  atomic 
oxygen  density,  and  neutral  temperature.  These  were  measured  by  a 
combination  of  three  instruments  on  the  satellites,  the  Open-Source 
Neutral-Mass  Spectrometer  (OSS),  the  Neutral-Atmosphere  Composition 
Experiment  (NACE),  and  the  Neutral-Atmosphere  Temperature 
Experiment  (NATE).  Each  of  these  instruments  and  their  operation  will  be 
briefly  discussed  below. 

Open-Source  Neutral-Mass  Spectrometer.  On  the  Atmosphere 
Explorer  satellites  the  Open-Source  Neutral-Mass  Spectrometer  provided 
measurements  of  the  molecular  nitrogen  density  and  atomic  oxygen 
density.  Open-source  spectrometers  can  best  detect  reactive  species  such  as 
oxygen  and  hydrogen  but  have  difficulty  detecting  particles  that  are  not 
approaching  along  the  instrument  axis.  This  instrument  has  been 
described  by  Neir  et  al.  (1973)  and  Burgess  and  Torr  (1987).  It  is  a  double- 
focusing  Mattauch-Herzog  magnetic-deflection  instrument  that  is  similar 
to  those  flown  on  earlier  sounding  rockets.  This  type  of  instrument  was 
chosen  for  compactness  and  low  power  consumption. 
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The  major  components  of  the  spectrometer  are  the  ion  source,  the 
electric  and  magnetic  analyzers  and  the  ion  collection  system.  The  ambient 
gas  enters  through  a  series  of  grids  that  isolate  the  electric  and  magnetic 
fields  of  the  ion  chamber  from  the  ambient  gas.  The  last  of  these  grids  acts 
to  prevent  any  ions  formed  in  the  chamber  from  exiting  through  the  grids. 
An  electron  beam  in  the  ion  chamber  produces  positive  ions  which  are 
focussed,  accelerated  and  collimated  before  entering  the  electric  and 
magnetic  analyzing  region.  The  ions  pass  through  the  electric  and  then 
the  magnetic  analyzers  where  energy  focusing  and  mass  separation  takes 
place.  The  two  detectors  of  the  OSS  would  give  mass  ranges  in  the  ratio  of 
1:8.  The  same  multiplier  detectors  were  used  for  each  region.  The  high 
mass  range  multiplier  would  be  turned  off  a  low  altitudes  and  a  grid  and 
electrometer  would  be  used  to  determine  the  high  mass  density.  This 
collector  system  allowed  the  density  measurements  to  be  made  over  a  large 
dynamic  range  caused  by  the  low  perigee  orbits. 

Neutral-Atmosphere  Composition  Experiment.  The  Neutral- 
Atmosphere  Composition  Experiment  is  a  closed  source  instrument  that 
complements  the  OSS  instrument  described  above.  For  5577 A  emission,  the 
same  constituents  measured  by  the  OSS  are  measured  by  the  NACE.  It 
produces  precise  measurements  of  non-reactive  species  and  total  oxygen 
(atomic-t-molecular)  due  to  surface  recombination  of  the  atomic  oxygen. 
The  NACE  can  provide  measurements  over  the  entire  forward  hemisphere 
independent  of  the  side  energies.  The  NACE  instrument  was  described  by 
Pelz  et  al.  (1973)  and  Burgess  and  Ton*  (1987).  The  system  was  based  on  the 
previously  flown  instruments  of  the  Orbiting  Geophysical  Satellite  (OGO-6) 
and  the  San  Marco-3  satellites.  The  spectrometer  had  a  mass  range  from 
one  to  46  amu  with  a  better  than  1  amu  resolution. 
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The  NACE  instrument  consisted  of  an  inlet  system,  an  ion  chamber, 
a  quadrupole  analyzer,  an  electron  multiplier  and  a  vacuum  pump  system. 
The  vacuum  pump  was  not  used  during  flight  but  was  used  to  maintain  a 
vacuum  in  the  instrument  prior  to  launch.  During  flight,  the  antechamber 
and  ion  source  are  connected  to  the  ambient  atmosphere  by  a  knife-edge 
orifice.  After  the  ambient  gas  is  thermalized  in  the  antechamber,  it  passes 
into  the  ion  chamber.  After  ionization,  the  gas  is  focused  and  accelerated 
into  the  quadrupole  analyzer.  The  ions  are  decelerated  and  separated  by 
their  charge  to  mass  ratio.  The  ions  are  then  accelerated  into  the  electron 
multiplier.  The  count  rate  is  proportional  to  the  chamber  density  of  the 
selected  gas.  An  additional  feature  of  the  NACE  instrument  was  a 
contaminant  cover  for  the  orifice  to  prevent  exhaust  products  from  the 
hydrazine  motors  from  entering  the  instrument  during  maneuvering. 

Neutral-Atmosphere  Temperature  Experiment.  The  Neutral- 
Atmosphere  Temperature  Experiment  was  used  to  provide  the  inputs  to 
allow  calculation  of  the  temperature  of  the  neutral  atmosphere  at  the 
location  of  the  satellite  and  the  density  of  molecular  nitrogen  outside  the 
spacecraft.  The  instrument  was  based  on  those  developed  for  the  San 
Marcos-3  and  AEROS  programs.  The  instrument  and  the  calculation  of  the 
temperature  and  density  are  discussed  by  Spencer  et  al.  (1973)  and  Burgess 
and  Torr  (1987). 

The  NATE  instrument  was  very  similar  to  the  NACE.  The  NATE 
featured  different  modes  of  operation  that  allowed  the  determination  of  the 
velocity  distribution  of  the  ambient  gas  either  by  using  the  spinning  motion 
of  the  satellite  or  in  the  despun  mode  by  using  a  baffle  that  could  be  stepped 
across  the  orifice.  The  method  of  density  and  temperature  determination  is 
given  by  Spencer  et  al.  (1973).  While  the  main  atmospheric  constituent  used 
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by  the  instrument  was  N2,  oxygen  and  hydrogen  were  studied  as  possible 
constituents  to  be  used  at  high  altitudes  to  determine  temperature. 

2.3.2  Ion  Measurements 

The  ionic  atmospheric  parameters  that  are  important  to  the 
production  of  5577A  emission  are  total  ion  or  electron  density,  electron  and 
ion  temperature  and  singly  ionized  molecular  oxygen.  A  combination  of 
four  instruments  on  the  satellites  provided  overlapping  measurements  of 
these  parameters,  the  Bennett  Ion-Mass  Spectrometer  (BIMS),  the 
Magnetic  Ion-Mass  Spectrometer  (BIMS),  the  Retarding-Potential  Analyzer 
(RPA),  and  the  Cylindrical  Electrostatic  Probes  (CEP).  Each  of  these 
instruments  are  briefly  described  below.  Also,  the  inter-calibration  between 
these  instruments  will  be  discussed. 

Bennett  Ion-Mass  Spectrometer.  The  Bennett  Ion-Mass 
Spectrometer  as  flown  on  the  Atmosphere  Explorer  satellites  was  derived 
from  those  used  on  the  Orbiting  Geophysical  Observatories  and  the 
Atmosphere  Explorer-B  satellites.  The  instrument  made  measurements  of 
thermal  positive  ions  between  1  and  72  amu  in  combinations  of  three  ranges 
with  densities  between  5  and  5  million  ions  per  cubic  centimeter  densities. 
The  ion  densities  measured  by  the  instrument  were  calibrated  by 
correlation  with  the  RPA  and  CEP  instruments  on  the  satellite.  The 
instrument  was  described  by  Brinton  et  al.  (1973)  and  Burgess  and  Torr 
(1987). 

The  Bennett  Ion-Mass  Spectrometer  was  located  on  the  forward- 
looking  part  of  the  satellite.  The  ions  are  drawn  in  from  the  ambient  gas  by 
a  negative  electric  field  at  the  orifice.  The  ions  are  then  accelerated 
through  the  spectrometer  gaining  energy  from  the  radio  frequency  fields. 
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passing  through  a  retarding  potential  to  the  collector.  The  flow  through  the 
instrument  was  enhanced  by  two  features,  the  use  of  a  multigrid  ion 
current  collector  instead  of  a  solid  plate  and  the  venting  at  the  rear  of  the 
instrument.  The  venting  was  especially  important  at  low  altitudes  by 
reducing  the  ram  buildup  in  the  instrument. 

Magnetic  Ion-Mass  Spectrometer.  The  Magnetic  Ion-Mass 
Spectrometer  was  designed  to  measure  the  abundance  of  the  ambient 
positive  ions  at  the  satellite  location.  The  spectrometer  flown  on  the 
Atmosphere  Explorer  satellites  was  similar  to  the  one  flown  on  the 
International  Satellite  for  Ionospheric  Studies  (ISIS-2)  spacecraft.  It  was  a 
magnetic  deflection-type  spectrometer  with  a  mass  range  of  1  to  90  amu. 
Like  the  BIMS,  the  MIMS  was  calibrated  in-flight  to  the  RPA  and  CEP 
instruments.  This  instrument  has  been  described  by  Hoffman  et  al.  (1973) 
and  Burgess  and  Torr  (1987). 

The  MIMS  instrument  was  a  small  magnetic  deflection  system 
consisting  of  an  entrance  aperture,  a  magnetic  analyzer  and  a  detector 
system.  The  entrance  was  positioned  normal  to  the  spin  axis  and  looking 
forward  in  the  despun  mode.  The  entrance  slit  formed  the  inlet  port  to  the 
instrument  and  used  ram  pressure  with  a  large  vent  to  minimize  internal 
pressure.  This  venting  was  adequate  to  allow  operations  down  to  120  km. 
The  instrument  had  three  detectors  that  could  collect  ions  in  the  ratio  of 
1:4:16.  The  detectors  could  be  used  to  study  specific  mass  ions  or  could 
study  the  spectrum  of  ions  by  sweeping  the  voltage. 

Retarding  Potential  Analyzer,  The  Retarding  Potential  Analyzer 
was  included  on  all  of  the  Atmosphere  Explorer  satellites.  It  is  also  known 
as  a  planar  ion  trap  and  was  an  improved  version  of  an  instrument  used  on 
the  Orbiting  Geophysical  Observatory  and  in  high  altitude  sounding 
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rockets.  Primarily  intended  to  measure  ion  temperature  and  ion 
concentrations  on  the  Atmosphere  Explorer  satellites,  it  also  was  used  to 
measure  ion  drifts  and  energy  spectra  of  the  thermal  and  supra-thermal 
electrons.  For  the  5577A  emission,  the  RPA  provided  electron  densities  and 
ion  temperatures.  The  other  ion  instruments  on  the  satellites  were  cross 
calibrated  with  this  instrument.  The  RPA  has  been  described  by  Hanson  et 
al.  (1973)  and  Burgess  and  Torr  (1987). 

The  instrument  consisted  of  four  sensors,  each  of  which  was  a  series 
of  grids  along  a  small  cylinder  and  a  solid  detector.  Sensor  heads  one  and 
four  were  located  forward-facing  when  the  satellite  was  despun  and  sensor 
heads  two  and  three  were  at  110  and  130  degrees,  respectively.  The  grids  in 
the  sensors  could  be  controlled  to  screen  the  different  particles  entering  the 
detector.  All  of  the  sensors  were  essentially  the  same  with  just  a  few  minor 
differences.  Collector  four  was  called  the  ion  drift  meter  and  was 
segmented  to  measure  off  axis  flows.  On  sensors  two  and  three,  an  extra 
grid  was  used  to  protect  the  sensors  from  ion  bombardment  when 
measuring  electrons.  In  addition,  the  heads  could  be  slightly  positively 
biased  to  measure  thermal  electrons. 

Cylindrical  Electrostatic  Probes.  The  Cylindrical  Electrostatic  Probes 
were  a  continuation  of  the  use  of  Langmuir  probes  used  since  1947  on  high- 
altitude  sounding  rockets  and  numerous  satellites,  including  Tiros, 
Atmosphere  Explorers  A  and  B,  and  the  International  Satellite  for 
Ionospheric  Studies.  The  instrument  was  used  to  measure  the 
temperature  and  concentrations  of  electrons.  It  consisted  of  two 
independent  probes,  one  mounted  parallel  to  the  velocity  vector  and  the 
other  perpendicular  to  it.  This  instrument  was  described  in  detail  by  Brace 
et  al,  (1973)  and  Burgess  and  Torr  (1987). 
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The  method  of  using  the  probe  to  measure  temperatures  and 
densities  was  described  by  Brace  et  al,  (1973).  The  probes  were  modified  for 
the  different  environment  the  satellite  would  encounter  at  low  altitudes. 
The  length  of  the  probe  was  reduced  by  a  factor  of  three  compared  to 
previous  probes  to  minimize  the  induced  voltages  from  the  motion  of  the 
instrument  through  the  geomagnetic  field.  The  diameter  of  the  probes 
were  increased  by  the  same  amount  to  add  a  heating  element  to  drive  off 
any  exhaust  contaminants  from  the  hydrazine  motors.  A  new  collector 
material  of  highly  oriented  tungsten  crystals  was  used  to  improve  the 
surface  potential  and  reduce  energy  smearing. 

Instrument  inter-calibration.  The  calibration  of  the  various  ion 
measuring  instruments  is  important  for  using  the  measurements  from  the 
Atmosphere  Explorer  satellites.  Since  the  instruments  provided  some 
redundancy  in  the  measurements,  if  a  parameter  from  a  given  instniment 
was  not  complete  another  instrument  could  be  substituted  in  its  place  if  the 
inter-calibration  was  known.  The  normalization  between  the  instruments 
was  based  on  using  the  Retarding  Potential  Analyzer  as  ground  truth  and 
then  determining  the  factors  for  the  other  instrument  to  make  their 
measurements  match  those  of  the  RPA. 

Prior  to  this  study,  the  normalization  factors  were  determined  as  a 
single  parameter  fit  for  the  instruments  on  satellite  AE-E  (Yee,  Private 
Communication).  These  factors  were  determined  using  a  ten  day  period  in 
the  satellite  lifetime.  The  normalization  factors  for  AE-C  and  AE-D 
instruments  were  assumed  to  be  the  same  as  for  AE-E.  Since  AE-E  did  not 
have  a  MIMS  instrument,  the  normalization  for  the  BIMS  from  AE-E  was 
used  on  the  other  satellites.  These  factors  were  multiplied  by  the 
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appropriate  satellite  measurement  to  get  a  corrected  value  to  match  the 
RPA  and  are  shown  in  Table  2.5. 

For  this  study,  these  normalization  factors  were  recalculated 
separately  for  each  satellite.  This  eliminates  the  assumption  that  each  of 
the  instruments  on  the  three  satellites  had  the  same  inter-calibration.  The 
RPA  instrument  on  each  satellite  was  assumed  to  be  the  correct 
measurement.  Measurements  from  each  of  the  other  instruments  were 
then  compared  to  the  RPA  as  shown  in  Figure  2.12,  and  a  straigl/  line  was 
fit  to  the  logarithm  of  the  values.  Table  2.6  shows  the  new  normalization 
factors  determined  for  the  various  satellites  and  instruments.  The  two 
parameters  are  Ubed  in  equation  [2.15], 


log,„(Mc)  =  S(log„(M„)]  -  I,  [2.15] 

where  Me  is  the  corrected  measurement  and  Mu  is  the  uncorrected 
measureii  .it.  These  correction  factors  were  calculated  for  the  elliptical 
period  of  each  satellite.  The  normalization  factors  <  \  ange  over  the  life  of 
the  satellite  as  shown  by  Figure  2.13  for  the  instruments  on  satellite  AE-E. 
If  data  is  to  be  used  from  other  periods,  the  normalization  factors  will  have 
to  be  updated  for  the  a^ipropriate  period. 
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Figure  2.12  Normalization  plots  for  satellite  C  for  a  single  day 
orbits. 
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Fi^re  2.13  Equivalent  RPA  densities  over  life  time  of  AE-E. 
Assumed  measurement  is  1x10®  cm'3  by  the  indicated 
instrument.  Each  point  uses  normalization  for  a  15  day 
period. 


CHAPTER  m 


THEORETICAL  MODEL  OF  5577A  OI  EMISSION 


3.1  Introduction 

In  this  chapter  the  development  of  a  theoretical  model  for  the 
production  of  0(lS)  will  be  discussed.  The  model  will  be  constructed  to  use 
data  from  the  Atmosphere  Explorer  satellites  if  available  to  calculate 
volume  emission  rate  profiles.  The  simulated  volume  emission  rate  (VER) 
profiles  will  be  used  to  study  the  various  sources  of  0(^S)  in  order  to 
determine  their  importance  to  the  dayglow  emission.  The  model  will  also 
be  used  to  study  the  individual  reactions  themselves. 

The  photochemical  model  requires  some  inputs  that  are  not  available 
from  the  satellites.  These  inputs  are  produced  by  other  models  and  include 
solar  flux,  neutral,  and  ionic  parameters.  The  models  used  have  been 
developed  for  and  used  in  a  number  of  investigations  and  can  be  used  to  fill 
in  gaps  in  the  satellite  data,  as  well  as,  provide  vertical  profiles  of  densities 
and  temperatures.  Figure  3.1  gives  a  general  flow  for  the  various  models 
used  in  the  production  of  OC^S). 
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Fiirui^  3.1  Simple  flow  diagram  showing  the  models  used  with 
the  satellite  data  to  calculate  the  production  of  0(1S). 
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Each  of  the  individual  models  will  be  discussed  separately  in  this 
chapter.  The  general  model  will  be  covered  and  examples  of  the  outputs 
shown.  The  use  of  the  individual  model  sources  to  supplement  or  replace 
the  satellite  data  will  be  covered  in  Section  3.3.  As  seen  in  Section  1.2,  there 

e 

are  numerous  reactions  that  may  be  applicable  to  the  5577A  dayglow 
emission.  The  selection  of  the  reactions  for  production  and  losses  in  the 

o 

photochemical  model  of  5577A  emission  will  also  be  discussed  in  Section 
3.2. 

3^  Photochemical  Modelling  of  5577A0I  Emission 

The  various  source  and  loss  reactions  that  have  been  suggested  for 
the  production  of  the  OC^S)  in  the  upper  atmosphere  were  discussed  in 
Chapter  I.  The  reactions  were  discussed  from  a  historical  standpoint.  To 
produce  the  photochemical  model,  the  reactions  were  examined  to 
determine  those  that  would  be  applicable  to  the  thermospheric  5577A 
emission.  The  reactions  concluded  of  possible  importance  are  listed  in 
Table  3.1.  These  source  and  loss  reactions  will  be  briefly  reviewed  then  the 
use  of  the  reactions  in  the  photochemical  model  will  be  discussed. 

The  major  source  of  upper  thermospheric  0(^8)  in  the  dayglow  is 
direct  excitation  of  0(3p)  by  photoelectron  impact  (Wallace  and  McElroy, 
1966), 


OC^P)  +  ep^^O('S)  +  ep  . 


[3.1] 


The  production  rate  for  this  reaction  is  a  function  of  atomic  oxygen  density, 
[0],  given  by 
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layer  0(^S). 
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P,=[0]j4(E^)o(EJdE,.  [3.2] 

The  integral  over  electron  energy  Eg  contains  the  photoelectron  flux,  fg,  and 
o,  the  cross  section  for  electron  impact,  both  of  which  are  functions  of 
electron  energy. 

As  discussed  in  Chapter  1,  the  cross  section  for  excitation  of  0(^8)  by 
electron  impact  has  been  studied  by  using  laboratory  measurements  (Shyn 
et  al.  ,1986;  Doering  and  Gulcicek,  1989)  and  from  theoretical  calculations 
(Smith  et  al.,  1967;  Henry  et  al.,  1969;  Thomas  and  Nesbet,  1975  and  Vo  Ky 
Lan  et  al.,  1972;  Berrington,  1991).  The  theoretical  cross  sections  are  lower 
than  the  measurements  of  Shyn  et  al.  (1986)  by  a  factor  of  approximately  2 
and  those  by  Doering  and  Gulcicek  (1989)  by  a  factor  of  about  1.5. 

Since  this  is  the  dominant  source  of  0(^8)  in  the  dayglow,  it  is 
important  to  know  the  magnitude  of  the  cross-section  for  this  reaction. 
Over  the  energy  range  of  1  to  100  eV,  a  correction  factor  for  the  magnitude  of 
the  parameterized  cross  sections  of  Green  and  Stolarski  (1972)  will  be 
determined  by  doing  a  least  squares  fit  to  the  VAE  observations.  This 
derived  cross  section  will  then  be  compared  to  those  determined  by  the 
previous  studies. 

The  dissociative  recombination  reaction  for  02'^, 

0;  +  e,— ^0(‘S)  +  0  ,  [3.3] 

is  an  important  source  of  green  line  emission  at  night  and  at  higher 
altitudes  in  the  dayglow  but  is  a  minor  source  near  the  peak  of  the 
emission.  The  production  rate  for  0(^S)  from  this  reaction  is. 
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P3=p3Y3[0;][eJ. 


[3.4] 


The  total  reaction  rate  as  a  function  of  electron  temperature, 


Ya  =  16x  10 


-7 


( 300^0  55 

vTe  y 


[3.5] 


is  taken  from  Torre/  al.  (1976).  The  amount  of  0(^S)  produced  by  the 
dissociative  recombination  is  dependant  on  the  vibrational  excitation  of  the 
02^  and  is  not  conclusively  known.  Mul  and  McGowan  (1979)  concluded  the 

total  rate  coefficient  is  not  substantially  affected  by  the  vibrational  excitation 
of  the  02'*'  ion.  The  quantum  3deld  of  0(1S),  Pa,  is  a  function  of  the 
vibrational  population  of  62^  and  can  be  expressed  in  terms  of  the  ratio  of 
electron  density  to  atomic  oxygen. 


/  \0.7 

R  -liilfMl 
[0]I,tJ 


[3.6] 


as  discussed  in  Abreu  et  al.  (1983)  and  Yee  et  al.  (1989). 

The  energy  transfer  reaction  of  N2  (A^Z^^)  by  atomic  oxygen, 

(A-’  i;)+0(=’P)^^N,  +  0('S),  [3.7] 

can  produce  a  major  portion  of  the  0(^S)  at  altitudes  near  the  peak  in  the 
thermospheric  emission  layer.  The  effective  rate  of  production  of  0(^8) 
from  this  reaction, 


P,  =  P,y,[N,(A^i:)][0], 


[3.8] 


89 


can  be  expressed  as  a  product  of  the  total  rate  coefficient  for  the  production 
of  all  states  of  atomic  oxygen  and  the  quantum  yield  or  branching  ratio  for 
the  amount  produced  in  the  0(^S)  state.  As  was  discussed  in  Chapter  1,  the 
branching  ratio,  P7,  is  not  well  known  and  there  is  still  some  disagreement 
in  the  total  rate  for  this  reaction,  Y7,  as  well.  Table  3.2  is  a  list  of  some  of 
the  rate  coefficients  and  branching  ratios  for  this  reaction  that  have  been 
determined  by  previous  studies. 

As  can  be  seen  from  the  discussion  in  Chapter  1  and  from  Table 
XXX,  this  reaction  has  been  repeatedly  accepted  and  rejected  as  an 
important  source  of  0(^S)  in  aurora.  The  importance  to  the  dayglow  5577 A 
emission  has  not  been  conclusively  determined.  Frederick  et  al.  (1976) 
assumed  the  values  for  the  total  rate  coefficient  and  branching  ratios  based 
on  auroral  studies.  Using  these  values,  a  small  contribution  to  the  total 
0(^S)  production  was  found  and  led  to  the  search  for  other  source  reactions 
in  the  dayglow. 

This  reaction  will  be  studied  in  two  parts.  First,  the  total  rate 
coefficient  for  the  quenching  of  N2(A)  will  be  determined  by  a  study  of  the 
3371 A  emission  of  N2  in  Chapter  4.  The  value  of  the  total  rate  coefficient 
will  be  compared  with  previous  values  determined.  In  Chapter  5,  this 
value  for  the  total  rate  coefficient  will  be  used  to  find  the  branching  ratio  for 
the  production  of  0(  ^S)  due  to  this  reaction.  An  effective  branching  ratio  for 
the  production  of  OC^S)  from  energy  transfer  from  N2(A)  under  dayglow 

conditions  will  be  determined  and  the  results  compared  to  recent  laboratory 
and  aeronomical  values. 


Total  Rate 

(X 10*^^) 

Branching 

Ratio 

Study 

Comments 

1  Laboratory  Values  | 

\)=0  3.5±0.6 
u=l  4.1±0.5 
\)=2  4.6±0.6 
u=3  5.2±0.8 

Thomas  and 
Kaufmann 

u=0  2.810.4 

u=l  3.410.6 

0.7510.13 

Piper  et  al. 

u=0  2.8 
u=l  3.3 
u=2  3.6 

_ 

0.75 

De  Souza 
et  al. 

Branching  ratio 
determined  from 
numerical  solution  of 
differential  equations 

Aeronomical  Values  | 

3.0 

0.10 

Frederick 
et  al. 

Branching  ratio  set  to 
match  auroral 
determination 

2.0 

0.29 

O'  Neil  et  al. 

Precede  artificial  aurora 
experiment  using 
Koneshea  et  al.  (1978) 
atmospheric  model 

12.0 

0.30 

0.45 

Sharp  and 
Torr 

For  -0=0,1 

For  0=0, 1,2 

Table  3.2  Summary  of  previous  values  of  the  total  rate  coefficient 
and  branching  ratios  used  for  Reaction  [3.7]. 
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Photodissociation  of  O2  whose  possible  importance  was  pointed  out  by 
Bates  and  Dalgarno  (1954), 

O^  +  hv— ^0('S)  +  0,  [3.9] 

is  a  very  minor  source  for  the  production  of  0(1S)  at  altitudes  above  140  km. 
However,  this  reaction  becomes  the  dominant  source  at  120  km.  The  cross 
sections  for  this  reaction  have  been  given  by  Lawrence  and  McEwan  (1973). 
The  rate  for  this  reaction  as  a  fimction  of  O2  density  is, 

Pg  =  [O2]  j  ra)qa)aa)dx,  [3.10] 

where  flk)  is  the  attenuated  solar  flux,  o(X)  is  the  photodissociation  cross 
section,  and  q(>.)  is  the  quantum  yield  of  0(1S)  all  as  a  function  of 
wavelength,  X.  An  expression  for  the  attenuated  solar  flux  is 

f(?i)  =  fjA)expj^-^Ja,.(\)n,(s)dsj,  [3.11] 

where  s  is  the  distance  along  the  path  from  the  sun  to  the  point  in  the 
atmosphere,  foo  iX)  is  the  unattenuated  solar  flux,  ciai(^)  is  the  total 
absorption  cross  section,  and  ni(s)  is  the  number  density  of  the  ith  species 

along  the  path.  The  species  included  in  the  solar  flux  attenuation  are 
atomic  oxygen,  O2  and  N2.  At  or  above  perigee  for  the  AE  satellites,  the 

attenuation  of  the  flux  for  this  reaction  is  small. 

Ion-atom  interchange  between  02"^  and  atomic  nitrogen, 

0;-(-N  0('S)-hN0\ 


[3.12] 
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was  considered  by  Frederick  et  al.  (1976)  as  the  required  additional  source  of 
0(^8)  at  low  altitudes.  The  rate  of  production  for  this  process  is, 


P,2={ii2Y,2[0^][N].  [3.13] 

A  rate  coefficient,  Y12,  of  2.5  x  lO'^^  cm^  s‘^  was  used  in  Frederick’s 
model  for  the  analysis  of  the  5577A  dayglow.  The  uncertainties  in  the 
nitrogen  density  gives  a  range  of  the  rate  coefficient  from  1.2  x  10' to  5.0  x 
lO'll  cm3  g-l  (Frederick  et  al.,  1976).  This  reaction  may  or  may  not  be 
significant  for  the  production  of  5577 A  emission  depending  on  the  solar  flux 
and  the  impact  cross  section  for  electrons  on  atomic  oxygen. 

The  loss  processes  of  0(^8)  at  the  level  of  the  thermospheric  emission 
are  the  radiative  transitions  to  the  0(lD)  and  0(3p)  states, 

0('8)— ^^0(‘D)  +  hv(5577A),  [3.14] 

0('8)— ^^0('P)  +  hv(2972A),  [3.15] 

with  transition  probabilities  of  1.06  s‘l  and  0.045  s’^,  respectively 
(Kernahan  and  Pang,  1975).  Losses  of  0(^8)  due  to  quenching  are  not 
important  at  the  altitudes  of  the  thermospheric  emission.  Quenching  of 
0(^8)  by  0(3p)  has  a  reaction  rate  of  2  x  lO'^'^  cm3  g-l  (Krauss  and 
Neumann,  1975)  and  can  be  ignored  in  modeling.  The  quenching  by  N2 
rate  of  about  5  x  10'^^  cm3  §-1  ^g  noted  in  Chapter  1  and  that  by  O2  is  a  very 

small  factor  above  about  125  km. 

A  photochemical  model  was  developed  using  the  reactions  discussed 
in  the  introduction  as  the  sources  and  losses  initially.  8ome  of  the  reactions 
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discussed  were  later  removed  from  the  model  since  they  could  not  produce 
enough  0(^S)  to  be  of  any  consequence  in  the  dayglow.  The  model  developed 
includes  the  following  reactions.  The  sources  of  0(^8)  are  photoelectron 
impact  on  atomic  oxygen  (1),  quenching  of  N2(A)  by  atomic  oxygen  (7), 
dissociative  recombination  of  02'''  (3),  and  photo  dissociation  of  O2  (9).  The 
only  losses  are  the  radiative  transitions  from  0(^8)  to  either  0(lD)  or  0(3p). 
Using  these  reactions  the  production  of  0(l8)  is 


Where  A5577  and  A2072  are  the  Einstein  transition  probabilities  for  the 
radiative  loss  processes  and  Pi,  P3,  P7,  and  P9  are  the  production  rates 
previously  discussed.  Ion-atom  interchange  process  in  Reaction  [3.12]  is 
included  in  the  model,  Equation  [3.16]  becomes, 

1^5577  ■^2972/ 

where  P12  is  the  production  of  0(18)  from  the  reaction.  The  importance  of 

0 

this  reaction  to  the  5577A  dayglow  emission  will  be  discussed  in  8ection  3.4. 

The  photochemical  model  calculates  the  volume  emission  rate  at 
each  given  point  in  space.  A  schematic  diagram  of  the  model  is  presented 
in  Figure  3.2  and  shows  the  general  flow  of  the  photochemical  model.  The 


Figure  3.2  Schematic  diagram  of  the  photochemical  model 
showing  the  various  possible  sources  of  the  atmospheric 
constituents  used  to  calculate  the  population  of  0(^S). 
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schematic  (Figure  3.2)  indicates  the  sources  for  the  constituents  required 
for  the  chemical  reactions.  The  model  employs  in-situ  atmospheric 
composition  measurements  from  the  AE  satellites  when  they  are  available. 
The  instruments  used  to  collect  the  data  were  discussed  in  Chapter  II. 
Data  were  obtained  for  most  of  the  chemical  species  from  more  than  one 
instrument  to  insure  the  most  complete  measured  profiles  possible. 
Atomic  oxygen  and  N2  densities  were  measured  by  the  Open-Source 

Spectrometer  (OSS)  (Nier  et  al.,  1973)  and  the  Neutral-Atmosphere 
Composition  Experiment  (NACE)  (Pelz  et  al.,  1973).  The  Neutral- 
Atmosphere  Temperature  instrument  (NATE)  (Spencer  et  al.,  1973) 
provided  temperature  measurements.  Ion  composition  data  were  collected 
by  the  Bennett  Ion-Mass  Spectrometer  (BIMS)  (Brinton  et  al.,  1973)  and  the 
Magnetic  Ion-Mass  Spectrometer  (MIMS)  (Hoffman  et  al.,  1973).  Electron 
temperatures  are  from  the  Cylindrical  Electrostatic  Probes  (CEP)  (Brace  et 
al.,  1973)  and  electron  densities  were  obtained  from  the  Retarding  Potential 
Analyzer  (RPA)  (Hanson  et  al.,  1973).  To  avoid  possible  inconsistencies 
between  the  data  from  the  various  sources,  the  ion  composition 
measurements  were  normalized  to  those  made  by  the  RPA.  This  was  done 
by  performing  a  least  squares  fitting  of  the  data  between  the  RPA  and  the 
ion  mass  spectrometers  individually  over  the  elliptical  phase  of  the 
satellite's  lifetime. 

When  an  in-situ  measured  parameter  is  not  available,  the 
photochemical  model  employs  the  Mass  Spectrometer  and  Incoherent 
Scatter  model  (MSIS-86)  (Hedin,  1987)  to  provide  neutral  densities  and 
temperatures  and  the  International  Reference  Ionosphere  (IRI)  (Belitza, 
1986;  1990a)  for  plasma  densities  and  temperatures.  When  these  model 
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parameters  are  used,  they  are  normalized  to  the  available  in-situ  satellite 
measurements. 

Since  both  of  the  important  reactions  being  studied  are  dependant  on 
the  solar  flux,  the  model  employs  a  solar  flux  spectrum  based  on  the 
Hinteregger  Reference  Spectra  (Hinteregger  et  al.,  1981)  and  is  scaled  by  the 
level  of  solar  activity.  The  solar  flux  below  250A  is  corrected  by  a  factor  of 
two  as  discussed  by  Richards  and  Torr  (1984).  The  photoelectron  spectrum 
is  produced  using  the  photoelectron  model  based  on  the  two  stream 
approximation  (Nagy  and  Banks,  1970).  The  photoelectron  model  is  used  to 
calculate  the  production  of  0(^S)  from  direct  excitation  of  atomic  oxygen 
and  the  N2(A)  population  due  to  both  direct  excitation  and  cascading  from 

the  higher  excited  states  of  the  nitrogen  molecule. 

For  a  given  model  run,  measured  or  modeled  values  for  the 
atmospheric  parameters  are  obtained  at  15  second  intervals  along  the  orbit, 
as  the  satellite  moves  towards  and  then  away  from  perigee.  At  each  point 
along  a  satellite  orbit  where  the  volume  emission  rate  is  determined,  a 
complete  profile  in  5  km  steps  between  80  and  500  km  is  generated  by  the 
MSIS-86  and  IRI  models  which  are  scaled  by  the  in-situ  measurements 
made  at  the  point  the  volume  emission  rate  is  being  determined.  These 
vertical  profiles  are  used  to  produce  the  photoelectron  fluxes  and  direct 
excitations  used  by  the  model  to  determine  the  volume  emission  rate  at  the 
point  of  measurement.  This  procedure  compensates  for  the  variation  in 
photoelectron  production  as  the  satellite  moves  along  the  path  to  different 
latitudes,  longitudes  and  solar  zenith  angles.  Once  the  profiles  are 
constructed,  the  model  determines  the  volume  emission  rate  at  5577A  for 
that  point  using  the  measured  in-situ  constituent  data.  This  procedure  is 
repeated  at  each  point  where  in-situ  measurements  are  made.  Since  each 
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point  is  calculated  independently,  a  profile  of  the  model  volume  emission 
rate  can  be  constructed  from  a  series  of  constituent  measurements  along 
the  satellite  track. 

As  mentioned  before,  the  model  is  constructed  to  be  run  if  some  of  the 
satellite  data  are  not  available.  The  photochemical  model  will  use  various 
atmospheric  model  to  fill  the  data  gaps  from  the  satellite.  These  models 
can  also  supply  all  of  the  needed  inputs  to  the  photochemical  model  to  allow 
it  to  run  independant  of  the  satellite.  These  various  models  will  be 
discussed  individually  in  the  next  section.  Examples  of  the  model  output 
will  be  given  in  section  3.4.  These  examples  will  give  the  results  both  with 
and  without  satellite  data.  This  model  will  be  used  to  examine  the  major 
reactions  responsible  for  producing  the  5577A  dayglow  emission.  In 
Chapter  VI,  the  model  will  be  used  to  examine  the  emission  profile  during 
twilight  conditions. 

3.3  Models  for  Atmospheric  Parameters 

The  photochemical  model  of  the  thermospheric  5577A  emission 
requires  inputs  that  either  are  not  available  from  the  satellite  observations 
or  are  missing  at  a  given  point  in  the  orbit  This  is  also  required  if  the 
model  is  to  be  used  without  satellite  input.  These  parameters  are  supplied 
by  separate  models  which  have  been  incorporated  into  the  photochemical 
model.  Where  possible  these  inputs  are  scaled  to  fit  the  observed 
parameters.  The  models  used,  discussed  in  this  section  separately,  include 
a  solar  flux  model,  the  Mass  Spectrometer  and  Incoherent  Scatter  (MSIS) 
model,  the  International  Reference  Ionosphere  (IRI)  model,  and  a 
photoelectron/electron  transport  model. 
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3.3.1  The  Solar  Flux  Model 

During  the  daytime  at  midlatitudes,  the  upper  atmospheric  process 
are  dominated  by  those  that  are  either  directly  or  indirectly  by  the  sun.  The 
solar  radiation  is  absorbed  by  the  gasses  of  the  atmosphere  which  are  then 
either  excited,  ionized  or  dissociated.  The  result  of  the  energy  deposited  in 
the  upper  atmosphere  is  heating  of  the  upper  atmosphere.  The  part  of  the 
solar  spectrum  that  is  most  important  to  the  thermosphere  is  the  x-ray  and 
ultraviolet  regions  below  about  2000A.  The  bulk  of  this  radiation  is  absorbed 
above  100  km. 

The  solar  radiation  shorter  than  lOOOA  has  enough  energy  to  ionize 
the  predominant  species  in  the  atmosphere.  While  the  region  between 
lOOOA  and  2000A  is  almost  completely  absorbed  in  the  dissociation  of  O2. 
This  short  wavelength  radiation  is  a  small  part  of  the  total  radiation  flux 
from  the  sum,  less  than  1  part  in  10^  (Ratclifle,  1972;  Wallace  and  Hobbs, 
1977).  But,  the  mass  of  the  absorbing  part  of  the  atmosphere  is  a  small 
fraction  of  the  total  atmospheric  mass,  so  the  input  of  energy  per  unit  mass 
is  extremely  large.  This  energy  input  becomes  even  larger  during  solar 
disturbances  when  wavelengths  less  than  about  lOOOA  can  increase  by 
orders  of  magnitude.  The  changes  in  solar  flux  due  to  solar  activity  are 
shown  in  Figure  3.3. 

o 

To  model  the  5577A  dayglow  emission  it  is  necessary  to  have  the  solar 
flux  as  a  direct  or  indirect  input  to  several  of  the  reactions.  The  model  used 
here  gives  the  solar  EUV  and  FUV  flux  at  the  top  of  the  atmosphere  in  the 
range  of  1  to  1750A  based  on  a  specification  of  the  level  of  solar  activity.  This 
solar  flux  is  used  to  calculate  the  direct  excitation  of  OC^S)  as  well  as  various 
ions  and  the  photoelectrons  that  can  be  used  to  produce  0(lS).  The  solar 
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Figure  3.3  Spectral  distribution  of  the  solar  irradiance  and  its 
variation  with  solar  activity.  This  representation 
emphasizes  the  x-ray  and  EUV  contribution  (After  Smith 
and  Gottleib,  1974). 
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flux  is  input  to  other  models,  which  will  be  discussed  later,  to  produce  the 
desired  quantities  for  modelling  the  5577 A  emission. 

The  model  has  several  methods  of  estimating  the  solar  flux 
depending  on  the  level  of  solar  activity,  each  of  these  methods  will  be  briefly 
described  here.  Details  of  these  methods  are  given  in  the  listing  in 
Appendix  A  along  with  the  various  data  sets  used  by  the  different  methods. 
After  the  different  methods  are  described,  examples  of  the  models  will  be 
shown  for  different  activity  cases. 

In  the  first  case,  the  flux  is  scaled  using  parametrization  methods 
based  on  the  daily  10.7  cm  radio  flux  and  the  81  day  centered  average  10.7 
cm  flux.  For  the  ionizing  EUV  flux,  Hinteregger’s  contrast  ratio  method 
(Hinteregger  et  al.,  1981)  is  used,  based  on  the  Torr  and  Torr  (1985)  bin 
Structure  for  reference  spectrum  SC#21REFW.  The  1026A  (Hydrogen 
Lyman  Beta)  and  335A  (Iron  XVI)  enhancement  ratios  are  calculated  firom 
Hinteregger's  formula.  The  flux  spectrum  is  then  calculated  from  these 
emissions  using  Hinteregger's  method.  For  Hydrogen  Lyman  alpha  line  at 
12 16 A,  the  correlation  relationship  derived  from  SME  data  by  Tobiska  and 
Barth  (1990)  is  used.  For  the  Schumann-Runge  continuum,  the  Torr  et  al. 

o 

(1980)  50A  bin  structure  is  used  but  the  coefficients  have  been  adjusted  to 
reflect  SME  and  LASP  rocket  data  (Rottman,  1981,  1988;  Mount  and 
Rottman,  1983,  1985). 

Another  method  of  determining  the  solar  flux  is  by  linear 
interpolation  between  high  and  low  activity  spectra.  This  method  is  based 
on  daily  10.7  cm  radio  flux  alone.  It  is  assumed  that  the  low  activity 
spectrum  corresponds  to  a  10.7  cm  flux  of  68  and  the  high  activity  spectrum 
to  243.  The  Hinteregger  SC#21REFW  and  F79050  spectra  as  binned  by  Torr 
&  Torr  (1985)  are  used  for  ionizing  EUV.  For  the  1050A-1350A  region,  the 
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SC#21REFW  spectrum  for  low  solar  activity  is  scaled  by  the  contrast  ratios 
to  obtain  the  high  activity  spectrum.  For  Lyman  alpha  and  the  Schumann- 
Runge  continuum,  linear  interpolation  is  essentially  the  same 
parametrization  as  above. 

The  next  two  methods  do  not  cover  the  spectral  range  of  the  previous 
methods.  To  keep  the  range  the  same  these  methods  use  one  of  the  above 
methods  and  then  replaces  the  values  for  their  spectral  range.  The  first  of 
these  methods  is  the  model  by  Tobiska  and  Barth  (1990)  which  uses  the 
linear  interpolation  method  to  complete  the  spectrum  of  solar  fluxes  and 
model  the  spectrum  between  32A  to  1050A.  The  procedure  involves  using 
10.7  cm  flux  and  the  Hydrogen  Lyman  alpha  flux,  considered  as 
representative  of  the  coronal/transition  region  and  the  chromospheric 
emissions,  respectively.  The  chromospheric  fluxes  are  calculated  from  a 
linear  correlation  relationship  with  the  Hydrogen  Lyman  alpha  flux. 
While  the  coronal  fluxes  are  calculated  from  a  combination  of  10.7  cm  flux 
and  a  "modeled"  10.7  cm  flux  based  on  Hydrogen  Lyman  alpha  flux.  The 
coronal  and  chromospheric  scaled  fluxes  are  then  summed  to  obtain  the 
value  of  the  solar  flux. 

The  last  method  uses  the  Hinteregger  contrast  ratio  method  for  the 
spectrum  and  the  Woods  and  Rottman  rocket  spectrum  for  10  November, 

o 

1988  for  the  region  from  300-1050A.  None  of  the  above  models  extends 
shortwards  of  18A.  An  amalgam  of  sources  are  used  to  estimate  the  flux  in 
this  region,  e.g.,  DeJager,  1964;  Smith  &  Gottlieb,  1974;  Manson,  1977; 
Kreplin  et  ai,  1977;  Horan  and  Kreplin,  1981;  Wagner,  1988.  For  the  first 
two  methods,  the  EUV  fluxes  between  250A  and  50A  are  normalized 
upwards  (after  Richards  and  Torr,  1984).  The  EUV  normalization 
coefficient  is  2.0  for  a  10.7  cm  flux  of  68  and  reduces  linearly  to  1.0  for  a  10.7 
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cm  flux  of  243.  Examples  of  the  solar  flux  are  shown  in  Figure  3.4  for 
different  solar  conditions  and  each  of  the  methods  above. 

3.3.2  The  Neutral  Atmosphere  Model 

As  can  be  seen  from  Figures  3.1  and  3.2,  the  modelling  of  0(^8 ) 
requires  a  neutral  atmosphere  model.  The  model  serves  two  roles,  the  first 
is  to  provide  density  and  temperature  profiles  for  use  in  the  photoelectron 
model  discussed  below.  These  profiles  are  scaled  to  match  the  satellite 
measurements  where  possible.  The  model  is  also  used  to  fill  in  for  the 
satellite  when  some  of  the  parameters  normally  measured  are  missing. 
The  neutral  model  can  also  be  used  to  calculate  the  profiles  of  emission 
above  and  below  the  satellite  or  in  cases  when  the  satellite  is  not  used. 

The  neutral  atmosphere  model  used  is  the  Mass  Spectrometer  and 
Incoherent  Scatter  (MSIS)  model  of  Hedin  (1983,  1987).  The  MSIS-86  model 
uses  data  collected  by  satellites,  rocket  measurements  and  ground  based 
incoherent  scatter  radars  to  produce  a  composite  model  of  the  densities  for 
selected  species  and  temperatures.  Table  3.3  is  a  summary  of  the  data 
sources  used  in  the  production  of  the  MSIS-86  model.  The  MSIS-86  model  is 
an  extension  of  the  MSIS-83  model  into  the  polar  regions  using  the 
additional  data  from  the  Dynamics  Explorer  satellite  and  the  addition  of 
atomic  nitrogen  as  an  output  species. 

In  the  MSIS  model,  data  are  fit  to  equations  for  the  temperature  and 
density  by  least  squares  methods.  The  fitting  equations  are  spherical 
harmonic  expansions  in  terms  of  local  time  and  geographic  latitude.  The 
coefficients  for  the  equations  are  dependant  on  solar  and  magnetic  activity 
as  well  as  season  (Hedin  et  ai,  1974).  The  temperature  profile  is  created  by 
joining  a  Bates  (1959)  temperature  profile  for  the  upper  thermosphere  to  an 
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Figure  3.4  Model  solar  flux  from  the  various  methods  for  (a) 
low  solar  activity  and  (b)  high  solar  activity.  [HR- 
Hinteregger  Contrast  Ratio;  LIN-  Hinteregger  Linear 
Interpolation;  LYA-  Lyman  a  scaling  of  Tobiska  and 
Barth;  WR-  Woods  and  Rottman  Rocket  Spectrum] 
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Table  3.3  List  of  data  sources  used  to  produce  the  Mass 
Spectrometer  and  Incoherent  Scatter  model  (After 
Hedin,  1983, 1987). 
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inverse  polynomial  profile  for  the  lower  thermosphere  (Hedin,  1983).  These 
profiles  are  joined  at  116.5  km  by  matching  the  temperature  and  the 
temperature  gradient. 

Using  the  temperature  profiles  and  and  average  density  that  is  a 
function  of  location  and  solar  and  magnetic  activity,  the  densities  are 
determined.  Using  functional  forms  of  the  temperature  profiles,  the 
hj'drostatic  equation  is  integrated  exactly  to  give  an  expression  for  the 
density  in  tenns  of  the  average  density.  To  account  for  the  homopause  near 
100  km,  the  models  give  density  profiles  that  are  sums  of  mixed  and 
diffusive  densities.  The  sum  is  multiplied  by  one  or  more  factors  to  allow 
for  changes  due  to  chemistry  and  dynamics  for  each  species.  Details  of  the 
methods  of  calculation  for  these  parameters  and  the  expansion  coefficient 
data  sets  are  given  in  Hedin  (1983,  1987)  and  are  based  on  the  methods  of  the 
OGO-6  model  of  Hedin  et  al.  (1974). 

The  MSIS  model  produces  densities  for  N2,  O2,  He,  O,  H,  Ar,  and  N 
at  a  given  altitude.  It  also  produces  two  temperatures,  the  exospheric 
temperature  and  the  neutral  temperature  at  the  altitude  desired. 
Comparison  of  the  model  to  subsets  of  the  data  used  to  produce  the  model 
show  reasonably  good  agreement  (Hedin,  1983’  1987).  While  the  model 
densities  can  vary  from  the  measurements  by  as  much  as  20%  and  the 
temperatures  by  as  much  as  40  K,  an  average  error  of  about  10%  in  the 
density  and  about  15  K  in  temperature  can  be  assumed.  Figure  3.5  gives 
examples  of  the  various  density  profiles  for  low  and  high  solar  activity. 

3.3.3  The  Ionospheric  Model 

Several  of  the  possible  source  reactions  for  0(^S)  require  knowledge  of 
temperatures  and  densities  of  the  ionosphere.  These  include  thermal 
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electron  density  and  temperature  and  ionic  species  densities.  As  discussed 
before  this  information  is  available  from  the  Atmosphere  Explorer  satellites 
but  a  model  is  needed  if  the  information  is  not  available.  The  model 
profiles  can  also  be  also  be  used  to  calculate  the  0(^8)  production  above  and 
below  the  satellite  track. 

To  provide  this  information,  the  International  Reference  Ionosphere 
(IRI)  model  is  used.  The  IRI  model  describes  monthly  averages  of  the 
electron  density,  electron  temperature  and  ion  composition  in  an  altitude 
range  from  50  to  1000  km.  The  model  is  designed  to  operate  under 
magnetically  quiet  conditions  in  non-auroral  regions.  The  IRI  model  has 
been  described  by  Rawer  et  al.  (1981).  Since  its  development  the  model  has 
been  extended  and  improved  by  the  addition  of  ground  and  space  data  as 
described  in  a  series  of  issues  of  Advances  in  Space  Research  (see  Volume 
2,  No.  10,  1982;  Volume  4,  No.  1,  1984;  Volume  5,  No.  7  and  No.  10,  1985; 
Volume  7,  No.  6,  1987;  Volume  8,  No.  4,  1988;  Volume  10,  No.  8  and  No.  11, 
1990).  The  newest  version  of  the  model  has  been  described  by  Belitza 
(1990a).  Table  3.4  gives  a  summary  of  the  various  data  sources  used  in  the 
development  of  the  IRI  model. 

The  IRI  model  can  be  separated  into  three  parts,  the  electron  density 
profile,  the  temperature  profiles  and  the  ion  density  profiles.  Each  of  these 
will  be  discussed  below.  The  model  uses  a  variety  of  analytical  functions 
and  expressions  to  represent  the  variations  of  the  densities  and 
temperatures.  In  most  cases,  these  are  in  the  form  of  spherical  harmonic 
function.  A  group  of  functions  known  as  the  Epstein  family  have  been  used 
to  represent  the  additional  variations  and  other  special  latitudinal  and 
diurnal  features  as  described  by  Rawer  (1987,  1988).  The  approach  of 
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Data  Sources 

Measurement 

lonosondes 

Peak  Plasma  Frequencies 

Incoherent  Scatter  Radars 
Jicamarca 

Arecibo 

St.  Santin 

Millstone  Hill 

Malvern 

E  Valley 

F2  Peak  Height 

F2  Bottomside,  Topside 
Electron  and  Ion  Temperature 

Allouette  1  &  2 

Topside  Sounder  Measurements 

Atmosphere  Explorer  C 

Electron  Temperatures 

Electron  temperatures 

AEROS 

Electron  and  Ion  Temperatures 

Ion  Composition 

Various  Rocket  Measurements 

D  &  E  Region  Electron  Density 
and  Ion  Composition 

Beacon  Satellites 

Ionospheric  Electron  Content 
Testing 

Ground  Based  Absorption 

D  &  E  i^gion  Variability 

Table  3.4  Data  sources  used  in  the  development  of  the 
International  Reference  Ionosphere  model  (After  Belitza, 
1990). 
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Booker  (1977)  is  used  to  construct  the  altitude  profiles.  Both  the  Epstein 
functions  and  the  Booker  application  are  discussed  in  Appendix  B. 

The  electron  density  profiles  have  1000  km  as  an  upper  limit  and 
either  60  or  80  km  as  a  lower  limit  for  da5dime  or  night  respectively.  It  is 
built  from  six  sub-regions.  These  sub-regions  are  joined  using  the  Booker 
approach.  The  middle  ionosphere  is  determined  by  various  layer 
parameters  such  as  peak  densities  and  heights,  thicknesses,  and  valley 
parameters  (Bilitza  and  Rawer,  1990). 

The  temperature  profiles  extend  from  120  km  to  1000  km  in  the  IRI 
model.  At  the  lower  boundary,  thermal  equilibrium  is  assumed  between 
the  neutral  species,  ions  and  electrons.  As  altitude  increases, 
temperatures  are  maintained  with  the  electron  temperature  greater  than 
the  ion  temperature  which  is  greater  than  the  neutral  temperature. 
Neutral  temperatures  are  derived  from  the  CIRA-72  model.  Both  the  ion 
and  electron  temperature  profiles  are  constructed  of  sub-intervals  which 
are  defined  using  Booker  functions.  The  methods  of  determining  the 
electron  temperatures  are  discussed  by  Bilitza  (1984,  1990)  and  the  ion 
temperatures  in  Bilitza  (1990). 

IRI  determines  the  ion  density  profiles  as  percentages  relative  to  the 
total  ion  density  or  the  electron  density  profile.  The  model  uses  the  Booker 
approach  to  determine  the  O'*’  and  02'''  profiles  and  fills  the  rest  of  the 
ionosphere  with  NO+  at  lower  altitudes  and  uses  H+  and  He"*"  at  higher 
altitudes.  At  higher  altitudes,  it  is  assumed  that  He'*'  is  10%  and  H+  is  90% 
of  the  remaining  ions.  The  ion  density  parameters  have  variation  with 
solar  zenith  angle  during  the  daytime  but  because  of  the  database 
limitations  has  no  variations  if  the  solar  zenith  angle  is  greater  than  90 
degrees.  Details  of  the  ion  composition  calculations  is  given  by  Bilitza 
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(1990a,  1990b).  Variations  of  temperatures  are  not  explicitly  included  in  the 
model  but  may  occur  through  the  CIRA  neutral  temperatures. 

The  IRI  model  is  designed  to  give  ionospheric  parameters  in  non- 
auroral,  magnetically  quiet  times.  IRI  has  been  shown  to  give  reasonable 
results  for  total  electron  content  (TEC)  for  mid-latitude  stations  but 
underestimates  daytime  low  latitude  TEC  (McNamara,  1984).  IRI  was 
tested  along  with  other  various  ionospheric  models  by  Brown  et  al.  (1984). 
It  was  shown  that  IRI  gives  reasonable  results  for  solar  minimum 
conditions  but  consistently  over  estimates  TEC  during  solar  maximum. 
Temperatures  from  IRI  are  in  good  agreement  with  observations  with  most 
discrepancies  attributed  to  the  insufficient  description  of  changes  with 
season  and  solar  activity  (Buonsanto,  1989;  Bilitza  and  Hoegy,  1990).  The 
following  Figures  3.6  and  3.7  show  profiles  of  densities  and  temperatures 
for  electron  and  ion  species  from  the  IRI  model  for  various  conditions. 

3.3.4  The  Photoelectron  Model 

During  the  daytime,  the  upper  atmosphere  is  dominated  by  solar 
driven  processes.  The  solar  flux  interacts  with  the  upper  atmosphere 
having  three  possible  consequences;  photodissociation  of  molecules, 
photoexcitation  of  atoms  and  molecules,  and  photoionization  of  atoms  and 
molecules.  The  photelectrons  produced  by  the  ionization  processes  have 
sufficient  energy  to  in  turn  dissociate,  excite  or  ionize  the  neutral 
atmosphere  further.  This  cascading  process  of  photoelectrons  is  shown 
schematically  in  Figure  3.8. 

The  photoelectron  model  used  is  essentially  the  one  described  by 
Solomon  (1987),  Solomon  et  al.  (1988),  and  Solomon  and  Abreu  (1989).  It 
uses  the  two  stream  electron  transport  model  of  Nagy  and  Banks  (1970)  to 
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Figure  3.6  Results  for  the  IRI  model  (a)  electron  temperature 
and  number  density  and  (b)  Ion  densities  for  low  solar 
activity. 
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Figure  3.7  Results  for  the  IRI  model  (a)  electron  temperature 
and  number  density  and  (b)  ion  densities  for  high  solar 


113 


Figure  3.8  Simple  flow  diagram  illustrating  the  cascading  of 
photoelectrons  from  photoionization  of  the  neutral 
atmospheric  constituents. 
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calculate  the  photoelectron  fluxes.  The  solar  flux  model  described  above  is 
used  as  an  input  at  the  top  of  the  atmosphere  and  is  attenuated  to  a  given 
altitude.  The  MSIS-86  model  can  be  used  to  generate  the  density  profiles  for 
the  major  neutral  species  for  attenuation  of  the  solar  flux  and  for 
ionization,  dissociation,  and  excitation  calculations.  The  photoionization 
cross  sections  are  from  Samson  and  Pareek  (1985)  and  Kirby  (1979).  The 
electron  impact  cross  sections  are  those  used  by  Solomon  (1988)  and  were 
calculated  using  the  analytic  fit  parameters  of  Green  and  Stolarski  (1972) 
and  Jackman  et  al.  (1977).  Appendix  C  is  a  detailed  discussion  of  the  model 
and  has  tables  showing  the  cross  section  parameters  used. 

The  electron  flux  is  calculated  for  the  upward  and  downward 
hemispherical  fluve''  along  a  magnetic  field  line  as  a  function  of  energy 
and  altitude,  j  effects  of  gravity  and  electric  and  magnetic  fields  are 
neglected.  The  model  does  not  allow  for  variations  in  the  pitch  angle,  so  the 
value  of  0.577  was  assumed  for  the  average  cosine  of  the  pitch  angle  to 
characterize  the  transport  as  recommended  by  Banks  et  al.  (1974).  Discrete 
energy  losses  from  inelastic  collisions  and  elastic  collisions  with  the 
neutral  atmosphere  are  considered.  The  backscatter  ratios  for  the  inelastic 
collisions  were  assumed  to  be  equal  to  those  for  elastic  collisions. 

To  calculate  the  electron  fluxes  at  a  given  point  in  the  atmosphere, 
the  electrons  are  binned  according  to  energy  in  variable  sized  energy  bins. 
Starting  with  the  highest  energy  electrons,  the  secondary  and  degraded 
primary  electrons  are  computed  and  added  to  the  electron  flux  at  each 
altitude  and  each  bin  of  lower  energy.  The  excitation  and  ionization  due  to 
this  energy  of  electrons  is  also  found.  Then  the  next  highest  energy  bin  is 
calculated  and  the  results  added  to  those  below  it  in  energy.  This  process  is 
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repeated  until  the  lowest  energy  bin  is  reached.  Figure  3.9  gives  examples 
of  the  modeled  fluxes. 

3.4  Model  Simulations  of  5577 A  Dayglow  Emission 

In  the  previous  sections  of  this  chapter  the  photochemical  model  of 
the  5577A  dayglow  emission  was  described.  Also,  details  of  several 
subsiderary  atmospheric  models  used  in  the  photochemical  model.  In  this 
section  some  examples  of  the  output  from  the  model  will  be  presented. 
These  examples  will  use  data  from  the  Atmosphere  Explorer  satellites. 
The  model  output  will  also  be  shown  using  input  profiles  from  the 
atmospheric  models  in  place  of  the  in-situ  measurements  of  the  satellites 
under  the  same  geophysical  conditions.  Figure  3.10  shows  the  relative 
importance  of  the  selected  reactions  in  a  typical  model  run.  Reaction  [3.12] 
is  included  in  this  figure. 

The  atmospheric  constituent  profiles  for  the  in-situ  measurements 
are  shown  in  Figures  3.11  and  3.12  for  the  up  and  down  legs  from  AE-C 
orbit  0584  and  AE-D  orbit  0912  respectively.  Orbit  AE-C  0584  varied  in  solar 
zenith  angle  from  75  to  41  degrees  during  the  down  and  up  leg  of  the  orbit 
respectively  with  an  F10.7  of  80.0  and  an  Ap  of  10.0.  For  AE-D  0912,  the 
solar  zenith  angle  changed  from  54  degrees  to  85  degrees  and  the  flO.7  and 
Ap  was  72.0  and  7.0  respectively.  The  resulting  5517 A  emission  from  the 
various  reactions  are  shown  in  Figures  3.13  and  3.14  for  the  photochemical 
model  output  from  the  in-situ  data. 

As  discussed  before,  the  production  from  the  two  major  reactions  in 
the  dayglow  has  not  been  conclusivly  established.  The  cross  sections  for 
electron  impact  on  atomic  oxygen  varies  by  as  much  as  a  factor  of  two 
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Figure  3.10  Plot  of  the  model  output  based  on  constituent  profiles 
provided  by  the  subsidiary  models  showing  the  emission 
due  to  the  various  production  reactions  of  0(^S) . 
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400 


119 


(ui5()  apn^i^iv 


Figure  3.12  In-Situ  measurements  of  atmospheric  constituents 
provided  by  the  various  instruments  for  orbit  AE-D  0912  up 
and  down  leg  of  the  orbit. 
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Figfure  3.14  Plot  of  the  model  output  for  orbit  AE-D  0912  based  on 
cor^tituent  profiles  provided  by  the  in-situ  measurements 
in  Figure  3.12. 
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between  the  theoretical  and  experimental  values.  The  quenching  of  N2(A) 
by  atomic  oxygen  has  both  the  total  rate  coefficient  and  the  0(lS)  branching 
ratio  in  question.  In  Chapter  4,  the  total  reaction  rate  for  the  quenching 
reaction  (Reaction  [3.7])  will  be  discussed  in  more  detail.  Chapter  5  will 
deal  with  the  cross  section  of  electron  impact  (Reaction  [3.1])  and  the 
branching  ratio  of  Reaction  [3.7]  using  the  total  rate  coefficient  determined 
in  Chapter  4. 


CHAPTER  rV 


RATE  FOR  ENERGY  EXCHANGE  FROM  N2(A)  BY  0(¥) 


4.1  Introduction 

As  shown  by  this  study,  one  of  the  important  mechanisms  for  the 
production  of  OC^S)  near  the  peak  of  the  dayglow  is  the  energy  transfer 
reaction, 


N2(A'i:)+0("P)— ^N^+OCS).  [4.1] 

There  is  disagreement  about  the  value  of  the  total  rate  coefficient  for  this 
reaction  and  for  the  branching  ratio  of  the  0(^8)  channel.  Details  of  the 
reaction  and  the  rate  coefficient  and  branching  ratio  were  given  in 
Chapters  1  and  3.  In  this  chapter,  a  separate  study  of  this  reaction  is  done 
using  a  model  that  gives  the  vibrational  distribution  of  the  excited  triplet 
states  of  N2.  This  model  is  used  in  conjunction  with  a  band  model  of  N2,  the 
photoelectron  model  discussed  in  Chapter  2  and  observations  from  the 
Atmosphere  Explorer  satellites.  The  value  of  the  total  rate  coefficient  is 
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determined  for  the  reaction  in  Equation  [4.1]  and  will  be  used  in  the  rest  of 
this  work. 

The  N2(A3Zu''")  state  is  the  lowest  excited  triplet  state  of  N2  with  a 
threshold  energy  of  6.5  eV.  Figure  4.1  shows  the  N2A  state  and  other 
excited  triplet  states  which  are  important  to  the  N2A  state.  This  state  is 
populated  by  two  mechanisms,  direct  excitation  from  the  ground  state  of  N2 
and  by  cascading  from  the  higher  N2  triplet  states,  involving  the  B^flg, 
W^Au,  B’^Su  .  C^riu.  These  population  processes  tend  to  fill  different 
vibrational  levels  of  the  N2A  state  (Ahmed,  1969;  Cartwright  et  al.,  1971; 
Cartwright,  1978).  The  direct  mechanism  fills  the  higher  vibrational  levels 
of  the  N2A  state  while  the  cascading  processes  fill  the  lower  levels. 

The  N2A  is  a  metastable  state  with  a  lifetime  of  1.36  to  2.6  S 
(Shemansky  and  Carleton,  1969;  Shemansky,  1969).  The  emission  due  to 
this  state  going  to  the  ground  state  is  known  as  the  Vegard-Kaplan  (VK) 
bands  and  is  weak  in  the  dayglow.  The  loss  of  N2A  is  due  to  quenching  by 
O,  O2,  N2  (Cartwright,  1978).  Above  about  120  km,  quenching  by  0(3p) 
would  be  an  important  removal  mechanism  for  N2A.  Figure  4.2  shows 
altitude  profiles  of  the  direct  and  cascading  sources  for  N2A  and  the  loss  by 
radiation  and  quenching  with  0. 

The  determination  of  the  rate  coefficient  for  this  mechanism  requires 
knowledge  of  the  distribution  of  the  vibrational  population  of  the  N2A  state. 
The  rate  coefficient  for  the  energy  transfer  reaction  has  been  determined  by 
Piper  (1982)  for  the  v  =  0,1  vibrational  levels  and  by  Thomas  and  Kaufmann 
(1985)  for  v  =  0, 1,2,3.  Both  studies  show  the  rate  coefficient  tends  to  increase 
as  the  population  of  higher  vibrational  levels  increases. 

In  addition  to  the  radiation  and  quenching  losses  of  the  N2A  state, 
the  vibrational  levels  can  have  processes  that  will  affect  the  population 
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Figure  4.1  Potential  energy  curves  for  the  lowest  triplet  states  of 

N2. 
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Figure  4.2  Sources  and  sinks  of  N2(A)  in  the  thermosphere. 
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distribution  of  the  N2A  state.  For  vibrational  levels  higher  than  7,  an 
inverse  transition  to  the  N2(B3ng)  has  been  suggested  by  Campbell  and 
Thrush  (1967)  and  Gilmore  (1969).  This  mechanism  can  account  for  the 
lack  of  emissions  from  the  N2A  states  higher  vibrational  levels  even  though 
these  levels  are  preferably  populated  by  the  direct  excitation  mechanism  as 
discussed  above.  This  mechanism  would  tend  to  remove  the  population 
from  the  higher  vibrational  levels  of  the  N2A  state  and  return  it  into  the 
lower  vibrational  levels  by  cascading  from  the  B^rig  state.  This  transition 
can  be  either  a  cascading  process  or  may  be  induced  by  collisions  as  in 
Equation  [4.2], 


N2  ( A"  i; ),  u  =  7  +  N2  N2  (B^^Rg ),  u  =  0  +  N2  +  AE .  [4.2] 

This  reaction  is  spin  allowed  and  should  be  a  fast  process  that  will  have 
AE=-0.0156  eV  (Dreyer  and  Perner,  1973).  The  modelling  of  the  inverse 
transition  mechanism  has  been  done  by  Cartwright  et  al.  (1971)  and 
Cartwright  (1978). 

Another  method  which  can  redistribute  the  vibrational  population  of 
N2A  is  a  relaxation  process  that  could  follow  the  reaction  in  Equation  [4.3], 

N2(A),Ui  >2  +  N2(X),u  =  0f^N2(A),Ui-2  +  N2(X),u=l,  [4.3] 

as  described  by  Dreyer  and  Pemer  (1973).  Setser  et  al.  (1970)  suggested  a 
two  quanta  process  to  account  for  the  observations  of  emissions  from  only 
the  V  =  0  and  1  levels  at  elevated  pressures  of  N2  which  would  correspond 
with  a  rapid  relaxation  of  the  bands  from  higher  vibrational  levels.  The 
fundamental  vibration  of  the  ground  state  N2(Xllg+)  is  2330.7  cm-1  which 
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is  nearly  equal  to  an  energy  difference  of  Av  =  2  for  IN 2 A  and  the  resonance 
becomes  a  better  match  as  the  vibrational  level  increases  (Dreyer  and 
Pemer,  1973).  This  process  may  become  important  at  lower  altitudes. 

4^  Modelling  of  N2  Triplet  States 

The  model  for  the  vibrational  population  of  the  N2  triplet  excited  state 
is  adapted  from  the  model  used  by  Cartwright  et  al.  (1971,1973)  and 
Cartwright  (1978)  for  the  auroral  studies  of  the  molecular  nitrogen 
emissions.  In  this  section,  the  model  will  be  briefly  described  and  a 
discussion  of  the  application  of  the  model  to  determine  the  total  rate 
coefficient  for  reaction  [4.1].  Examples  of  the  output  from  the  model  will  be 
shown  in  general  and  for  the  N2A  state  in  general. 

The  vibrational  population  model  takes  into  account  the  excited  states 
lying  within  13  eV  of  the  ground  state.  The  cross  sections  for  excitation  to 
energies  of  higher  states  become  very  small.  Using  this  energy  level  as  a 
cutoff  gives  seven  excited  triplet  states  to  be  considered  (A^Lu'*’,  B^flg,  W^Au, 
C^Ifu,  E^Zg"*",  and  D^Zu"^)  which  are  shown  in  Figure  4.1.  The 
population  mechanisms  for  the  excited  states  are  direct  excitation  from  the 
ground  electronic  state  (X^Zg"*’,  v=0)  and  by  various  intrasystem  cascading 
processes.  Collisional  transfer  excitation  can  be  neglected  above  100  km  but 
may  become  important  at  lower  altitudes  (Hays  and  Oskam,  1973; 
Cartwright,  1978). 

The  vibrational  population  model  has  been  described  by  Cartwright  et 
al.  (1971,1973)  and  Cartwright  (1978).  In  this  model  the  populations  of  the 
triplet  states  are  governed  by  statistical  equilibrium  as  the  entire  system 
cannot  be  characterized  by  a  kinetic  temperature.  By  assuming  the  only 
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direct  excitation  to  the  triplet  states  is  from  the  ground  state  v=0,  the 
statistical  equilibrium  equation  for  the  relative  number  density  (ni“/no^)  of 
the  v'  vibrational  level  of  the  a  excited  state  can  be  written  as, 


k;,+siASr.^ 


[4.4] 


where  no^  is  the  number  density  of  the  ground  state  N2  (Cartwright,  1978). 
The  Greek  and  Arabic  symbols  denote  the  excited  electronic  states  and 
vibrational  levels  respectively.  The  first  term  on  the  left  is  the  direct 
excitation  to  the  v'  level  of  the  a  excited  state  from  the  ground  state  with 
kv'o“  being  the  excitation  rate  from  electron  impact  in  s-l.The  second  term 
on  the  left  represents  the  cascading  transitions  from  the  i^^  vibrational  level 
of  the  state  to  to  the  v’  vibrational  level  of  the  a  state  with  AijP“  as  the 
Einstein  transition  probability  between  levels  and  states.  The  terms  on  the 
right  side  of  Equation  [4.4]  represent  the  losses  of  the  population  of  the  v' 
level  of  the  a  state.  These  are  the  spontaneous  transitions  to  the  r^h  level  of 
the  6^^  state  and  the  collisional  quenching  of  the  state,  by  species  such  as 
0,  n,  O2,  NO,  and  N2,  Qv 

The  population  of  each  vibrational  level  for  certain  excited  states  are 
coupled  to  those  in  some  other  electronic  state  by  dipole  allowed  intrasystem 
cascading  terms  of  Equation  [4.4].  A  set  of  93  coupled  equations  is  produced 
when  each  level  is  expressed  as  a  separate  equation  given  by  Equation  [4.4]. 
This  set  of  equations  may  be  solved  iteratively  for  the  populations  relative  to 
the  ground  state  lowest  level.  Iteration  is  used  to  solve  the  equations  since 
the  states  are  initially  populated  by  direct  excitation  and  then  allowed  to 
evlove  into  steady  state  populations  by  intrasystem  cascading.  To  solve  the 
set  of  equations,  knowledge  of  the  photoelectron  flux  which  is  assumed  to 
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change  slowly  compared  to  the  time  it  takes  the  system  to  reach  steady  state 
populations,  the  excitation  cross  sections,  the  Einstein  transition 
probabilities,  and  the  quenching  rates  are  required.  Each  of  these  items 
will  be  briefly  discussed  below. 

The  calculation  of  direct  excitation  rate  from  the  ground  state  can  be 
expressed  as. 


k“o  =  q“oJf(eKo{eM®' 


[4.5] 


where  qv'o“  is  the  Franck-Condon  factor  between  the  ground  state  v=0  and 
the  v’  level  of  the  state,  fte)  is  the  photoelectron  flux  as  a  fxinction  of 
energy,  Ov’o®(e)  is  the  rotationally  averaged  excitation  electron  impact  cross 
section  for  the  v'  level  of  the  ati'  state. 

The  Franck-Condon  factors  between  the  electronic  states  are  those 
used  by  Cartwright  (1970;  1978);  Cartwright  et  a/.(1971,  1973).  By  using  the 
v=0  of  the  ground  state  as  the  only  source  level  for  the  direct  excitations,  the 
use  of  the  Franck-Condon  approximation  is  accurate  to  better  than  10% 
(Miller  and  Krauss,  1967;  Cartwright,  1970). 

The  photoelectron  spectrum  used  by  the  vibrational  population  model 
is  generated  by  the  photoelectron  model  discussed  in  Chapter  3.  The 
upward  and  downward  fluxes  are  summed  to  give  the  total  photoelectron 
flux  from  0  to  245  eV  and  are  binned  in  continuously  variable  sized  bins 
ranging  from  0.50  to  11.93  eV.  The  range  of  photoelectrons  more  than 
adequately  covers  the  energies  for  the  impact  cross  sections  for  excitation  of 
the  triplet  states  which  peak  between  10  and  20  eV  and  fall  off  quickly 
beyond  30  eV. 
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The  electron  impact  excitation  cross  sections  for  the  excited  states  are 
those  used  by  Cartwright  (1978).  These  cross  sectiop^;,  u  gained  by 
Cartwright  et  al.  (1977),  do  not  contain  contributions  from  cascading 
processes.  This  allows  the  cascading  contribution  to  be  spp''*rated  from  the 
direct  process  for  individual  study  (Cartwright,  1978).  For  this  work  the 
impact  cross  section  for  the  D3Su+  state  is  assumed  to  be  the  same  as  the 
cross  section  for  the  E3lg+,  both  being  small  relative  to  the  lower  state  cross 
sections  since  they  are  the  lowest  Rydberg  states  for  N2. 

The  transition  probabilities  used  in  this  work  are  the  same  as  those 
used  by  Cartwright  (1978).  Thirteen  sets  of  probabilities  are  used  including 
three  for  the  inverse  intrasystem  transitions.  The  ground  state  is  only 
connected  to  the  triplet  excited  states  through  the  A3Su+  level  which  give 
rise  to  the  Vegard-Kaplan  bands  and  are  of  interest  here.  The  transitions 
used  for  the  triplet  states  are  shown  in  Figure  4.3.  In  Figure  4.3,  the  three 
inverse  transitions  are  shown  as  dashed  lines. 

The  quenching  of  the  excited  N2  triplet  states  is  well  known  only  for 
the  and  the  B^flg  states.  Estimates  can  be  made  of  the  effects  of 

quenching  on  the  populations  of  the  other  excited  states  based  on  orbital 
similarities  (Cartwright,  1978).  For  this  work,  quenching  is  only  important 
for  N2(A3Iu+)  since  it  is  a  metastable  state  but  it  is  also  important  for  the 
W^Au  state.  Of  the  species  available  for  quenching,  only  0  is  important  for 
quenching  since  the  densities  of  the  other  species  is  low  above  130  km, 
especially  for  N  and  NO.  Along  with  decreasing  densities,  O2  has  a 
quenching  rate  an  order  of  magnitude  smaller  than  that  of  O.  N2  may  be 
important  at  lower  altitudes  for  vibrational  relaxation  of  the  N2A  state  as 
discussed  before  but  for  low  vibrational  levels  the  quenching  rate  is  very 
small. 


Figure  4.3  Transit 
vibrational  pop 
shown  as  dashe< 
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Figure  4.4  shows  the  general  flow  of  the  vibrational  population  model 
in  the  dayglow.  The  figure  indicates  the  sources  of  the  input  data  at  the 
various  steps.  Typical  resultant  populations  of  the  three  lowest  triplet 
states  is  shown  in  Figures  4. 5-4. 7.  In  Figure  4.5,  the  N2A  state  shows  the 
population  is  primarily  in  the  lowest  vibrational  levels.  Also  in  this  figure, 
the  inverse  transition  to  the  N2B  state  is  seen  at  vibrational  levels  above  6. 
The  shapes  of  the  N2A  and  N2B  states  (shown  in  Figure  4.6)  are  determined 
by  both  direct  production  and  cascading  processes  while  the  N2C  state  in 
Figure  4.7  is  only  due  to  direct  excitation. 

4,3  DeterminationofTotal  Rate  Coefficient 


4.3.1  Background 

The  determination  of  the  total  rate  coefficient  for  Equation  [4.1]  is 

o 

approached  by  a  study  of  the  3371 A  emission  from  the  second  positive 
transition  of  N2(C3nu,  n=0)  to  N2(B3ng,  u=0).  The  3371A  emission  was 
measured  by  the  Visible  Airglow  Experiment  on  the  Atmosphere  Explorer 
satellites.  The  measurements  of  the  2'^^  positive  transitions  were 
contaminated  by  the  Vegard  Kaplan  bands,  primarily  due  to  the  (0-9) 
transition.  Figure  4.8  is  a  low  resolution  spectrum  of  the  N2  dayglow 
emission  from  3000  to  4000A.  The  overlap  of  the  bands  at  3371 A  is  shown. 

The  sources  of  the  population  of  the  N2(A)  state  were  discussed  above. 
The  N2(C)  state  in  the  dayglow  has  photoelectron  excitation  as  the  only 
source  with  minimal  cascading  from  higher  states.  Since  the  second 
positive  transitions  are  dipole  allowed,  quenching  is  not  important  to  the 
N2(C)  state.  The  use  of  337 lA  emission  is  ideal  to  study  N2  densities  and  as 
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Figure  4.4  Diagram  showing  the  general  flow  of  the  vibrational 
population  model  with  the  various  inputs  and  sources  of 
data. 
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Figure  4.6  Resultant  population  of  the  N2(A)  state  for  the 
vibrational  levels  less  than  12  from  the  model. 
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Figure  4.6  Resultant  population  of  the  N2(B)  state  for  the 
vibrational  levels  less  than  12  from  the  model. 
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Fi^re  4.7  Resultant  population  of  the  N2(C)  state  for  the 
vibrational  levels  less  than  12  from  the  model. 
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Figure  4.8  Low  resolution  spectrum  from  3000  to  4000  A  with  the 
Second  Positive  and  Vegard-Kaplan  bands  (After  Conway 
and  Christensen,  1984). 
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an  indicator  of  the  photoelectron  fluxes  (Kopp  et  al.,  1977;  Hernandez  et  ai, 
1983;  Conway, 1983;  Conway  and  Christensen,  1985).  Both  Kopp  et  al.  (1977) 
and  Hernandez  et  al.  (1983)  concluded  that  the  observations  by  the 
Atmosphere  Explorer  satellites  consistently  overestimated  the  modelled 
emission  at  3371 A  (Conway,  1983). 

By  using  the  model  for  the  vibrational  population  for  N2  described 
above  and  a  band  model  for  the  nitrogen  molecule,  simulated  brightnesses 
for  the  positive  and  Vegard-Kaplan  (0-9)  emissions  can  be  made  and  the 
results  compared  with  the  measured  brightnesses  of  the  AE  satellites. 
Since  atomic  oxygen  will  be  the  only  important  quencher  for  the  N2(A)  state, 
by  adjusting  the  quenching  rate  the  modelled  emissions  can  be  made  to 
agree  with  those  measured  by  the  instrument.  This  will  give  the  total  rate 
coefficient  in  the  dayglow  for  Equation  [4.1], 

4.3.2  Method  of  Determination 

The  measurements  of  3371 A  emission  by  the  VAE  instrument  also 
contained  emission  from  the  (0-9)  Vegard-Kaplan  band.  Using  these 
measurements,  the  quenching  of  the  N2(A),u=0  state  may  be  determined  by 
separating  the  Vegard-Kaplan  band  contribution  from  the  3371A 
brightness.  For  this  purpose,  orbits  were  selected  in  which  the  satellite  was 
despun  with  channel  two  of  the  VAE  instrument  oriented  radially  outward 
from  the  earth.  Data  were  limited  to  orbits  which  remained  south  of 
auroral  latitudes  to  eliminate  sources  for  the  N2(C)  state  other  than 
photoelectron  impact. 

o 

The  3371  A  brightness  data  were  taken  using  channel  2  (wide  angle) 
of  the  Visible  Airglow  Experiment  (VAE)  on  Atmosphere  Explorer  satellite 
(AE-C)  (Hays  et  al.,  1973)  which  was  described  in  Chapter  2.  The  center  of 
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the  transmission  was  at  3371  A  for  a  filter  temperature  of  19.3  and 
shifted  to  longer  wavelengths  by  0.1577  A  per  oC.  The  transmission  of  the 
filter,  which  peaks  at  46.4%,  is  shown  in  Figure  4.9  along  with  each  of  the 
two  bands  normalized  to  give  a  total  intensity  of  1.0.  Atomic  oxygen  and 
molecular  nitrogen  densities  were  acquired  from  the  Open  Source  Mass 
Spectrometer  (OSS)  (Nier  et  al.,  1973)  aboard  AE-C. 

Using  a  band  model  for  molecular  nitrogen,  the  second  positive  and 
Vegard-Kaplan  bands  can  be  produced  for  a  given  rotational  temperature. 
These  S3mthesized  bands  are  shown  in  Figure  4.10  with  the  total  intensity  of 
both  of  the  bands  normalized  to  1.0  for  a  rotational  temperature  of  500  K. 
Increasing  the  rotational  temperature  broadens  the  band  structure  and 
will  alter  the  contribution  of  the  bands  to  a  measurement  of  the  337 lA 
emission.  Since  the  transmission  function  of  the  filter  shifts  with  its 
temperature,  the  fractional  contribution  of  the  bands  to  a  given 
measurement  can  be  written  as, 

0(T„.T„,)  =  s|| r{>-,.T„)f,pa(>.„T„)dX,  + 1 T(Xj,T„,)fv„(Xj,T„)dX,|,  [4.6] 

where  S  is  the  sensitivity  of  the  instrument,  T  is  the  filter  transmission 
function  of  the  filter  at  the  filter  temperature,  and  f  is  the  line  intensity  for 
the  band.  Tf,i  and  Trot  are  the  filter  and  rotational  temperatures 
respectively.  Figures  4.11  and  4,12  show  the  contribution  factors  of  the  N2 
second  positive  and  Vegard-Kaplan  bands  which  are  observed  as  a  fimction 
of  filter  and  rotational  temperatures  from  the  two  terms  in  Equation  [4.6]. 

The  vibrational  population  model  for  the  N2  triplet  states  described 
above  provides  the  steady-state  populations  for  the  N2(A),u=0  and  N2(C),u=0 
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Figure  4.9  The  VAE  3371  A  filter  transmission  for  AE-C  with 
the  Second  Positive  and  Vegard-Kaplan  bands  shown 
normalized  to  1. 


Figure  4.10  The  results  of  the  band  model  for  molecular 
nitrogen  for  the  two  bands  normalized  to  1  at  a  rotational 
temperature  of  500  K.  (a)  V-K  (0-9)  (b)  2PG  (0-0). 
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excited  states.  These  populations  can  be  used  to  determine  the  emission 
rate  as  a  function  of  altitude  for  the  particular  transitions  of  interest  (VK  0- 
9  and  2PG  0-0).  The  emission  rates  can  be  combined  with  the  terms  for  the 
contribution  for  each  band  in  Equation  [4.6]  and  integrated  along  a  line  of 
sight  path  to  calculate  the  observed  brightness  as  measured  by  the  VAE 
instrument, 


)  =  S| || e.,p^T(X.  )f.pc  I e JfvK  ,)dX,  |dl .  [4.7] 

In  Equation  [4.7],  <I>  is  the  integrated  surface  brightness,  e  is  the  emission 
rate  for  the  bands  from  the  vibrational  population  model,  and  the  integral 
would  be  from  the  satellite  to  infinity  along  the  line  of  sight  of  the 
instrument.  The  other  variables  are  defined  as  in  Equation  [4.6]. 

The  rate  for  quenching  of  the  N2(A)  state  by  atomic  oxygen  can  be 
foimd  by  using  Equation  [4.7]  to  model  the  3371A  observations  from  the  VAE 
instrument.  The  quenching  will  afTect  only  the  contribution  from  the 
Vegard-Kaplan  emissions  to  the  integrated  brightness.  By  adjusting  the 
modelled  brightness  to  best  reproduce  the  observations  the  quenching  rate 
may  be  determined. 

The  value  of  the  quenching  rate  may  be  determined  by  finding  the 
value  that  will  give  the  minimum  between  the  data  and  modelled 
emission.  The  data  for  the  x^  calculation  was  taken  from  four  orbits  from 
satellite  AE-C,  with  the  altitudes  between  150  and  165  km.  This  altitude 
range  was  selected  to  avoid  data  from  higher  altitudes  where  the  quenching 
is  not  important  and  the  lower  altitudes  where  vheicle  glow  contamination 
may  be  included  in  the  data.  This  contamination  problem  will  be  briefly 
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covered  later.  The  results  of  the  determination  will  be  discussed  in  the  next 
section  and  comparisons  with  previous  values  will  be  given.  The  rate 
coefficient  will  be  used  in  Chapter  5  to  determine  the  branching  ratio  for  the 
production  of  0(^S)  by  Equation  [4.1]. 

4.3.3  Results 

Orbits  from  the  AE-C  satellite  were  used  to  determine  the  quenching 
rate  of  the  N2(A)  state.  Using  the  vibrational  population  model  for  these 
orbits,  Figures  4.13  and  4.14  show  the  altitude  profiles  of  the  emissions  rate 
for  the  Vegard-Kaplan  (0-9)  and  the  Second  Positive  (0-0)  bands.  The 
quenching  rates  used  for  the  N2(A)  u=0  level  shown  in  the  figures  are  from 
Beiting  and  Feldman  (1979)  and  Piper  and  Caledonia  (1981).  The  effects  of 
the  quenching  rate  on  the  N2(A)  state  can  be  seen  by  the  change  in  emission 
in  the  two  figures. 

The  I'esults  of  the  n^inimization  is  shown  in  Figure  4.15  with  a 
total  rate  coefficient  of  2.75x10'^^  cm^  s'^.  The  observed  337lA  emission  is 
simulated  using  these  emission  profiles  in  Figures  4.13  and  4.14.  These 
have  been  compared  to  the  observations  from  AE-C  for  orbits  457  and  1911  as 
shown  in  Figures  4.16  and  4.17.  The  smaller  quenching  rate  which  agrees 
with  the  values  of  Piper  and  Caledonia  (1981)  gives  the  best  fit  to  the 
observations. 

The  model  tends  to  under  predict  the  emission  at  low  altitudes  near 
perigee  for  orbit  1911.  The  discrepancy  may  be  due  to  vehicle  glow  which 
has  been  studied  for  the  Atmosphere  Explorer  satellites  (Yee  et  al.,  1984) 
and  recently  for  the  Space  shuttle  (Viereck  et  al.,  1991).  The  residual  of  the 
brightness  for  orbit  1911  is  a  linear  function  of  [N2]^[0]  (private 
communication  Yee)  as  shown  by  Figure  4.18.  The  vehicle  glow  reactions 
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Figure  4.13  Plot  of  emission  rates  from  the  model  for  orbit  AE-C 
0457  based  on  perigee  conditions  for  the  two  bands.  The 
effects  of  the  change  in  quenching  for  the  VK  (0-9)  is 


350 


148 


{[uy\)  epmuiv 


Figure  4.14  Plot  of  emission  rates  from  the  model  for  orbit  AE-C 
1911  based  on  perigee  conditions  for  the  two  bands.  The 
effects  of  the  change  in  quenching  for  the  VK  (0-9)  is 


Fisfure  4.16  The  results  of  the  determination  of  the  total  rate 
coefficient  for  Reaction  [4.1]. 
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Figure  4.16  Plot  of  integrated  brightness  measurements  from 
the  model  and  as  measured  for  orbit  AE-C  0457.  The 
effects  of  the  change  in  quenching  for  the  VK  (0-9)  on  the 
integrated  brightness  is  shown. 
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Figure  4.17  Plot  of  integrated  brightness  measurements  from 
the  model  and  as  measured  for  orbit  AE-C  1911.  The 
effects  of  the  change  in  quenching  for  the  VK  (0-9)  on  the 
integrated  brightness  is  shown. 
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Fi^re  4.18  Comparison  of  the  residual  airglow  brightness 
versus  the  [N2]^[0]. 
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have  not  been  confirmed  but  this  dependance  would  support  the  suggested 
reactions  for  vehicle  glow  suggested  by  Green  (1984)  and  Green  et  al.  (1986) 
and  is  similar  to  the  results  of  Meyerott  and  Swenson  (1990). 

The  rate  of  quenching  for  Equation  [4.1]  determined  by  this  study 
supported  the  value  of  2.8xl0'll  cm'3  s'^  as  determined  by  Piper  and 
Caledonia  (1981)  and  DeSouza  et  al.  (1985)  for  the  u=0  vibrational  level.  This 
value  will  be  used  in  the  following  chapters  to  determine  the  cross  section 
for  electron  impact  on  atomic  oxygen  and  the  branching  ratio  of  Equation 
[4.1]  for  the  production  of  0(^8). 


CHAPTER  V 


ANALYSIS  OF  5577A  OI DAYGLOW  EMISSION 

5.1  Introduction 

As  was  discussed  in  Chapters  1  and  3,  the  major  reactions  for 
production  of  0(lS)  in  the  dayglow  have  parameters  which  are  not  known 
well.  The  total  rate  coefficient  for  the  quenching  of  N2(A)  by  0(3p)  (Reaction 
[4.1])  was  determined  in  Chapter  4  from  observations  of  the  3371  A  emission 
by  the  Visible  Airglow  Experiment  on  AE-C.  The  electron  impact  cross 
section  of  atomic  oxygen  (Reaction  [3.1])  and  the  branching  ratio  for  the 
quenching  reaction  (Reaction  [4.1])  have  large  variations  between  various 
studies  of  these  reactions.  In  this  chapter,  these  two  parameters  will  be 

o 

determined  by  a  study  of  the  5577 A  dayglow  emission. 

The  model  developed  in  Chapter  3  will  be  used  to  determine  the 
magnitude  of  the  electron  impact  cross  sections  for  atomic  oxygen  and  the 
branching  ratio  for  the  quenching  of  N2(A).  The  total  rate  coefficient  for  the 
quenching  reaction  from  Chapter  4  will  be  used.  Observations  of  the  0(lS) 
emission  by  the  Atmosphere  Explorer  satellites  are  used  with  the 
photochemical  model  to  determine  the  two  parameters  by  least  squares 
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techniques.  The  effect  of  the  vibrational  population  of  the  N2(A)  state  on  the 
production  of  0(lS)  will  be  discussed.  Since  the  chemistry  for  the 
production  of  0(1S)  in  the  dayglow  must  extend  into  the  twilight,  in  the  final 
part  of  this  chapter  the  values  for  the  various  reactions  determined  for  the 

o 

dayglow  will  be  used  to  model  the  5577A  emission  in  the  twilight. 

5^  Determination  of  Reaction  Parameters 


5.2.1  5577A  Dayglow  Model 

The  photochemical  model  was  developed  in  Chapter  3  using  some  of 
the  reactions  discussed  in  Chapter  1  as  the  sources  and  losses.  Some  of  the 
reactions  discussed  in  Chapter  1  were  not  used  in  the  model  since  their 
production  of  OC^S)  was  not  of  any  consequence  in  the  dayglow.  The  model 
developed  for  the  dayglow  includes  the  following  reactions.  The  sources  of 
OC^S)  are  photoelectron  impact  on  atomic  oxygen  (Reaction  [3.1]), 
dissociative  recombination  of  02'*'  (Reaction  [3.3]),  quenching  of  N2(A)  by 
atomic  oxygen  (Reaction  [3.7]),  photo  dissociation  of  O2  (Reaction  [3.9]),  and 
ion-atom  exchange  (Reaction  [3.12]).  The  only  losses  are  the  radiative 
transitions  from  0(lS)  to  either  0(^D)  or  0(3p).  Using  these  reactions  the 
production  of  0(lS)  is 


[ocs)] 


{P.+P3+P7+P9+P12} 

[  Akc77  -I-  A  2079} 


[5.1] 


Where  A5577  and  A2972  are  the  Einstein  transition  probabilities  for  the 
radiative  loss  processes  and  Pi,  P3,  P7,  Pg^  and  P12  are  the  production  rates 
previously  discussed  in  Chapter  3. 
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The  photochemical  model  calculates  the  volume  emission  rate  at 
each  given  point  in  space.  A  schematic  diagram  of  the  model  as  used  in 
Equation  [5.1]  is  presented  in  Figure  5.1  for  the  source  reactions  of  0(^S) 
used  in  this  Chapter.  The  schematic  (Figure  5.1)  indicates  the  constituents 
required  for  the  chemical  reactions  and  the  possible  by  products  that  result 
from  the  reactions.  The  model  employs  in-situ  atmospheric  composition 
measurements  from  the  AE  satellites  where  possible. 

Data  were  obtained  for  most  of  the  chemical  species  from  more  than 
one  instrument  to  insure  the  most  complete  measured  profiles  possible. 
These  instruments  were  described  in  Chapter  2.  Atomic  oxygen  and  N2 
densities  were  measured  by  the  Open-Source  Spectrometer  (OSS)  (Nieref 
a/.,  1973)  and  the  Neutral- Atmosphere  Composition  Experiment  (NACE) 
(Pelz  et  ai,  1973).  The  Neutral-Atmosphere  Temperature  instrument 
(NATE)  (Spencer  gf  a/.,  1973)  provided  temperature  measurements.  Ion 
composition  data  were  collected  by  the  Bennett  Ion-Mass  Spectrometer 
(BIMS)  (Brinton  et  al.,  1973)  and  the  Magnetic  Ion-Mass  Spectrometer 
(MIMS)  (Hoffman  et  al.,  1973).  Electron  temperatures  a^e  from  the 
Cylindrical  Electrostatic  Probes  (CEP)  (Brace  et  al.,  1973)  and  electron 
densities  were  obtained  from  the  Retarding  Potential  Analyzer  (RPA) 
(Hanson  et  al.,  1973). 

When  an  in-situ  measured  parameter  is  not  available,  the 
photochemical  model  employs  the  Mass  Spectrometer  and  Incoherent 
Scatter  model  (MSIS-86)  (Hedin,  1987)  to  provide  neutral  densities  and 
temperatures  and  the  International  Reference  Ionosphere  (IRI)  (Belitza, 
1986)  for  plasma  densities  and  temperatures.  These  models  were  discussed 
in  Chapter  3.  When  these  model  parameters  are  used,  they  are  scaled  to 
the  available  in-situ  satellite  measurements. 
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Since  both  of  the  important  reactions  being  studied  are  dependant  on  the 
solar  flux,  the  model  employs  a  solar  flux  spectrum  based  on  the 
Hinteregger  Reference  Spectra  (Hinteregger  et  ai,  1981)  and  is  scaled  by  the 
level  of  solar  activity.  The  solar  flux  below  250A  is  corrected  by  a  factor  of 
two  as  discussed  by  Richards  and  Torr  (1984).  The  photoelectron  spectrum 
is  produced  by  a  model  based  on  the  two  stro:  oproximation  (Nagy  and 

Banks,  1970).  The  photoelectron  model  is  used  tc  c  '.culate  the  production  of 
0(1S)  from  direct  excitation  of  atomic  oxygen  and  the  N2(A)  population  due 

to  both  direct  excitation  and  cascading  from  the  higher  excited  states  of  the 
nitrogen  molecule.  The  solar  flux  model  was  described  in  Chapter  3  along 
with  the  photoelectron  model  which  v/as  described  by  Solomon  cfaZ.  (1988) 
and  by  Solomon  and  Abreu  (1989). 

For  a  given  model  run,  measured  or  modeled  values  for  the 
atmospheric  parameters  are  obtained  at  15  second  intervals  along  the  orbit, 
as  the  satellite  moves  towards  and  then  away  from  perigee.  At  each  point 
along  a  satellite  orbit  where  the  volume  emission  rate  is  determined,  a 
complete  profile  in  5  km  steps  between  80  and  500  km  is  generated  by  the 
MSIS-86  and  IRI  models  which  are  scaled  by  the  in-situ  measurements 
made  at  the  point  the  volume  emission  rate  is  being  determined.  These 
vertical  profiles  are  used  to  produce  the  photoelectron  fluxes  and  direct 
excitations  used  by  the  model  to  determine  the  volume  emission  rate  at  the 
point  of  measurement.  This  procedure  compensates  for  the  variation  in 
photoelectron  production  as  the  satellite  moves  along  the  path  to  different 
latitudes,  longitudes  and  solar  zenith  angles.  Once  the  profiles  are 
constructed,  the  model  determines  the  volume  emission  rate  at  5577A  for 
that  point  using  the  measured  in-situ  constituent  data.  This  procedure  will 
be  repeated  at  each  point  along  the  satellite  path  where  in-situ 
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measurements  are  made.  Since  each  point  is  calculated  independently,  a 
profile  of  the  model  volume  emission  rate  can  be  constructed  from  a  series 
of  constituent  measurements  along  the  satellite  track. 

5.2.2  Observations  of  5577A  Dayglow 

The  observations  of  the  5577A  emission  were  taken  by  the  Visible 
Airglow  Experiment  (VAE)  photometer  on  the  Atmosphere  Explorer  (AE) 
satellites  (Hays  et  al.,  1973)  which  were  described  in  Chapter  2.  The 
measurements  of  the  5577A  airglow  employed  in  this  study  were  obtained 
by  the  AE-C  and  -D  satellites,  which  occupied  elliptical  orbits  during  the 
first  year  of  their  lifetimes.  Selected  orbits  from  this  period  provide  the 
emission  data  which  were  used  in  the  development  and  validation  of  the 
photochemical  model. 

Data  were  selected  for  elliptical  orbits,  when  the  spacecraft  was  de- 
spun.  For  the  de-spun  normal  mode,  the  high-sensitivity  channel  (channel 
2)  of  the  VAE  instrument  pointed  radially  out  from  the  Earth  as  the  satellite 
moved  along  its  track.  This  orientation  allowed  the  measurements  to  be 
made  without  significant  contamination  from  either  vehicle  glow  or  by 
Rayleigh  scattered  light  from  below  the  satellite  (Yee  and  Abreu,  1983).  In 
an  elliptical  orbit  it  is  possible  to  obtain  a  profile  of  the  emission  as  the 
satellite  intersects  layers  of  the  atmosphere  at  different  altitudes.  Data 
were  selected  for  mid-latitude  and  equatorial,  daytime  conditions  (latitude 
between  ±60  degrees  and  solar  zenith  angle  <  85  degrees). 

Volume  emission  rates  were  calculated  from  the  brightness 
measurements  in  steps  of  4  km  from  the  satellite  perigee  altitude  up  to  400 
km.  It  was  assumed  that  the  emission  layer  is  horizontally  homogeneous 
and  there  are  no  changes  in  galactic  background  intensity  between  two 
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successive  measurements,  by  the  instrument.  The  measurements  of  the 
column  brightness  made  in  this  fashion  can  be  inverted  by  a  simple 
derivative  technique  schematically  illustrated  in  Figure  5.2  which  was 
discussed  in  Chapter  2. 

An  example  of  the  measurements  of  the  5577A  emission  and  the 
inverted  volume  emission  rate  profile  from  the  derivative  "inversion" 
technique  are  presented  in  Figure  5.3.  The  figure  shows  the  measured 
surface  brightness  for  AE-C  orbit  0455  as  the  satellite  moves  away  from 
perigee.  As  the  satellite  moves  away  from  perigee  the  latitude  changes 
from  30®  N  to  15®  N  and  the  solar  zenith  angle  changes  from  69®  to  about  66®. 

5.2.3  Method  of  Determination 

Photochemical  model  simulations  have  been  carried  out  using  orbits 
from  the  AE  satellites  to  determine  the  magnitude  of  the  impact  cross- 
section  of  electrons  on  atomic  oxygen  (Reaction  [3.1])  and  the  branching 
ratio  for  0(^8)  by  quenching  of  N2(A)  by  atomic  oxygen  using  the  total  rate 
from  Chapter  4  (Reaction  [3.7])  that  will  best  reproduce  the  measured  5577A 
emission  profiles.  Table  5.1  is  a  list  of  the  orbits  used  along  with  some 
orbital  characteristics.  All  of  the  orbits  were  used  to  find  the  parameters 
for  Reactions  [3.1]  and  [3.7]  that  will  best  represent  the  emission  profile  for 
0(^8)  using  data  between  perigee  and  250  km.  Above  250  km,  dissociative 
recombination  (Reaction  [3.3])  becomes  as  important  as  Reaction  [3.1]  for 
the  production  of  0(^8).  The  best  fit  is  determined  by  finding  the 
parameters  that  will  give  the  minimum  deviation  for  the  orbits  listed. 
The  equation  for  is 
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Figure  5.2  The  geometry  used  for  the  brightness  measurements 
of  the  5577A  emission. 
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inverted  volume  emission  rates  for  AE-C  0455U. 
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Satellite 

Orbit 

Date 

Alt 

LST 

Leg 

C 

455 

74026 

155.5 

15.17 

34.0 

U/ 

C 

584 

74038 

157.8 

12.90 

41.8 

U/D 

C 

606 

74040 

158.2 

12.54 

42.2 

U/D 

C 

628 

74042 

158.5 

12.13 

44.3 

D 

C 

650 

74044 

159.0 

11.75 

45.4 

D 

C 

656 

74044 

159.0 

11.64 

56.0 

U/D 

C 

3354 

74269 

135.3 

7.36 

-21.1 

U 

C 

3410 

74273 

136.1 

6.36 

-28 

U 

D 

672 

75337 

142.5 

9.39 

57.1 

D 

D 

814 

75349 

142.8 

8.53 

26.5 

D 

D 

816 

75349 

143.0 

8.52 

22.6 

D 

D 

818 

75349 

143.0 

8.51 

37.3 

D 

D 

846 

75351 

141.8 

8.43 

20. 

D 

D 

912 

75357 

141.6 

7.94 

0.4 

D 

D 

918 

75357 

141.6 

7.90 

8.0 

D 

D 

942 

75359 

4 

D  1 

D 

984 

75363 

143.0 

7.50 

-14.2 

VfD 

D 

1021 

76001 

144.6 

7.27 

-13.5 

U/D 

D 

1045 

76003 

146.0 

7.13 

-23.5 

U/D 

D 

1069 

76005 

147.4 

.99 

-32.9 

U/D 

D 

1094 

76007 

149.3 

6.84 

-27.2 

U/D 

D 

1190 

76015 

142.0 

6.28 

-45.0 

U 

D 

1192 

76015 

142.2 

6.27 

-35.9 

U 

Table  6.1  List  of  orbits  used  from  Satellites  AE-C  £md  AE-D  in 
this  study  along  with  some  of  the  parameters. 


where  n  is  the  number  of  data  points,  ag  is  the  error  in  the  VER  from  the 
measurements  of  the  satellite,  Om  is  the  errors  attributed  to  the  model.  Vs  is 
the  VER  from  the  satellite,  and  V^  is  the  VER  determined  by  the  model. 
Correction  factors  for  the  magnitude  of  the  electron  impact  cross  section 
and  values  of  the  effective  branching  ratio  for  the  production  of  OC^S)  by 
N2(A)  quenching  are  determined  that  will  best  represent  the  data.  The 
model  VER  profile  was  interpolated  to  be  on  the  same  altitude  grid  as  the 
inverted  data,  then  the  total  deviation  for  all  of  the  selected  orbits  were 
calculated. 

The  total  rate  coefficient  for  Reaction  [3.7]  determined  in  Chapter  4  is 
for  the  n=0  vibrational  level.  This  value  for  the  total  rate  coefficient  is  in 
close  agreement  with  the  values  for  u=0  of  Piper  (1982)  and  of  Thomas  and 
Kaufman  (1985).  The  branching  ratio  is  varied  from  0.0  to  1.0  and  the  value 
to  give  the  best  fit  to  the  observed  data  is  determined.  The  effective  rate  of 
production  of  0(^S)  will  then  be  the  product  of  the  total  rate  coefficient  and 
the  best  fit  branching  ratio  for  the  reaction  from  the  simulations. 

The  magnitude  of  the  electron  impact  cross  section  for  electrons  on 
atomic  oxygen  was  determined  by  using  Equation  5.2.  The  magnitude  of 
the  parameterized  cross  section  of  Green  and  Stolarski  (1972)  is  varied  by 
factors  from  0.5  to  3.5  in  steps  of  0.25.  This  allows  the  cross  section 
magnitudes  to  vary  from  half  the  values  of  Henry  et  al.  (1969)  to  almost 
twice  those  of  Shyn  et  al.  ( 1986).  The  value  of  the  multiplying  factor  that 
gives  the  best  fit  is  determined  by  this  method. 
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The  magnitude  of  the  cross  section  and  the  yield  of  0(^8)  from 
Reaction  [3.7]  were  varied  simultaneously  to  find  the  best  fit  parameters  by 
a  least  squares  fitting  to  the  inverted  data.  The  result  is  a  plot  of  the 
magnitude  of  the  a  function  of  the  multiple  factor  for  the  impact  cross 
section  and  the  effective  branching  ratio  for  a  given  assumed  total  rate 
coefficient.  An  example  of  results  for  determining  of  the  model  parameters 
are  shown  in  Figure  5.4.  The  next  section  will  discuss  the  results  of  the 
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parameters  determined  here  as  they  affect  the  dayglow  5577A  emission. 

5^  Analysis  of  Results  for  5577ADay^ow 


5.3.1  Effects  of  the  Vibrational  Population  of  N2(A3Zu'^) 

The  total  rate  coefficient  for  the  quenching  reaction  is  dependant  on 
the  distribution  of  the  population  of  the  vibrational  levels  of  the  N2(A)  state. 
As  noted  before,  the  vibrational  levels  are  preferentially  populated  by  direct 
and  cascading  processes.  The  results  from  laboratory  studies  and 
aeronomical  determinations  of  the  total  rate  coefficient  cannot  be  directly 
compared  unless  the  vibrational  population  is  taken  into  account.  The 
effective  branching  ratio  for  the  reaction  can  be  determined  if  a  value  of  the 
total  rate  coefficient  is  assumed  and  the  distribution  of  the  vibrational 
population  is  known. 

The  production  of  0(  ^S)  by  the  quenching  of  N2(A)  from  Reaction  [3.7] 
can  be  expressed  as, 


P[0('S)1  =  ZP>4N=(A)1[0(’p)]  =  P,„k„[N,(A)L[o(>P)] 


[5.3] 


NORMALIZED  CHI  SQUARE  FOR  AE-C  AND  AE-D 


lee 


Figure  6.4  Example  of  the  normalized  for  the  orbits  listed  in 
Table  5. 


167 


where  P  is  the  branching  ratio  and  k  is  the  rate  coefficient  for  the  reaction. 
Allowing  for  the  change  in  rate  coefficient  by  vibrational  level  as  shown  by 
Piper  et  al.  (1981),  De  Souza  et  al.  (1985),  and  Thomas  and  Kaufman  (1985), 
solving  Equation  [5.3]  for  the  effective  branching  ratio  gives, 

=  iP.f.r.  ,  [5.41 


where  fy  is  the  fractional  population  of  the  vibrational  level  and  r-o  is  the 
ratio  of  the  rate  coefficient  for  the  level  and  the  total  rate  coefficient. 

The  vibrational  population  of  the  N2(A)  state  was  calculated  from  the 
model  by  Cartwright  et  al.  (1971,  1973)  and  Cartwright  (1978)  which  was 
adapted  to  used  the  photoelectron  flux  spectrum  from  a  two  stream  model 
for  photoelectron  fluxes  (Nagy  and  Banks,  1970).  The  vibrational  population 
model  used  statistical  equilibrium  with  both  direct  and  cascading  processes 
to  calculate  the  population  of  the  lowest  triplet  states  of  N2.  The 
photoelectron  model  was  discussed  in  Chapter  3  and  the  vibrational 
population  model  was  described  in  detail  in  Chapter  4. 

The  results  for  the  distribution  of  the  vibrational  population  of  the 
N2(A)  state  at  150  km  is  shown  in  Figure  5.5.  The  figure  shows  the 
population  is  highest  for  the  two  lowest  vibrational  levels.  The  rapid 
decrease  in  population  for  vibrational  levels  above  x>=l  due  to  the  inverse 
transition  to  the  N2(B)  state.  This  population  distribution  is  important  in 
determining  the  production  of  0(lS)  due  to  Reaction  [3.7].  Even  though  the 
quenching  rate  increases  with  vibrational  level  as  shown  by  Piper  (1982) 
and  Thomas  and  Kaufman  (1985)  the  rapid  decrease  in  population  indicates 
the  lowest  levels  are  important  for  the  production  of  0(^S). 


Fractional  Population 
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NjCA)  Vibrational  Population  at  150  km 
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Figure  6.6  Fractional  population  from  the  N2  Triplet  model  at 
150  km. 
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Using  the  population  levels  from  Figure  5.5  and  the  results  from  Piper  et  al. 
(1981)  and  Piper  (1982)  for  the  u=0,l  levels,  solving  for  the  effective 
branching  ratio  as  given  by  Equation  [5.4]  gives  a  value  of  0.312.  This  value 
does  not  agree  well  with  the  one  determined  from  the  yp-  fitting  studies 
performed  in  this  work.  This  indicates  that  the  higher  vibrational  levels 
are  required  to  fully  determine  the  0(^8)  production  from  the  quenching 
reaction. 

5.3.2  Dependance  of  the  Determined  Parameters  on  Other  Reactions 

o 

In  addition  to  the  major  reactions  for  the  production  of  the  5577A 
dayglow,  there  are  several  minor  sources  which  could  play  a  role  in  the 
determined  magnitude  of  the  reactions  from  above.  The  two  of  most 
importance  at  the  altitude  regions  of  the  Atmosphere  Explorer  observations 
used  in  this  study  are  the  dissociative  recombination  of  02"'’  (Reaction  [3.3]) 
and  the  ion  exchange  between  atomic  nitrogen  and  02'''  (Reaction  [3.12]). 
Both  of  these  reactions,  which  were  discussed  in  detail  in  Chapters  1  and  3, 
and  will  be  important  to  the  twilight  transition  period  (discussed  in  the 
following  section)  and  the  nightglow. 

These  reactions  are  used  to  determine  the  sensitivity  of  the 
determined  parameters  for  Reactions  [3.1]  and  [3.7].  The  two  parameters 
are  determined  for  different  combinations  of  the  rate  coefficient  for  Reaction 
[3.12]  and  the  branching  ratio  for  the  production  of  0(^8)  by  Reaction  [3.3]. 
For  the  branching  ratio,  in  addition  to  the  calculated  value  determined  as 
described  in  Chapter  3,  the  value  was  fixed  at  6%,  which  is  a  median  value 
of  the  branching  ratio  listed  in  several  studies  (Zipf  ,1970,  1980;  Kopp  et  al., 
1977;  Abreu  et  al.,  1983;  Guberman  1988).  The  rate  coefficient  for  Reaction 
[3.12]  is  not  well  known.  A  value  of  2.5  xlO  ^l  cm^  s'^  was  determined  by 
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Frederick  et  al.  (1976)  within  a  factor  of  two  due  to  uncertainty  in  the  actual 
atomic  nitrogen  density  profile.  For  this  study  the  value  was  set  at  2.5  xlO" 
cm^  s'l  and  then  multiplied  by  0,  0.5,  1.0,  and  2.0  to  remove  the  reaction 
from  the  model  completely  and  allow  the  rate  coefficient  to  vary  over  the 
limits  stated  by  Frederick  et  al.  (1976). 

The  two  parameters  determined  in  Section  5.2.3  were  found  for  each 
case  described  above.  The  results  of  the  fitting  were  similar  to  the  one 
shown  in  Figure  5.4  with  some  small  variations  in  the  values  previously 
determined.  The  results  are  summarized  in  Figure  5.6  for  the  eight  cases 
along  with  the  values  determined  by  the  previous  studies  for  comparison. 
These  results  will  be  summarized  below  along  with  the  rest  of  the  model 
analysis. 

5.3.3  Summary  of  Dayglow  Results 

The  photochemical  model  we  have  developed  is  able  to  reproduce  the 

o 

observed  daytime  5577A  thermospheric  emission.  Four  major  source 
reactions  were  incorporated  into  the  model  of  the  emission  profile.  The 
reactions  needed  for  the  model  are  photoelectron  impact  on  atomic  oxygen, 
dissociative  recombination  of  02'*',  quenching  of  N2(A3z^+)  by  atomic 

oxygen,  and  photo-dissociation  of  O2  which  becomes  more  important  below 
the  peak.  The  ion-exchange  reaction  between  02'’'  and  atomic  nitrogen 

o 

reaction  is  not  needed  to  explain  the  5577A  dayglow  emission  profile. 

In  this  section,  the  results  of  the  determination  of  the  magnitude  of 
the  impact  cross  section  for  photoelectrons  on  atomic  oxygen  and  the 
branching  ratio  for  the  quenching  of  N2(A)  will  be  summarized.  The 
vibrational  population  distribution  of  N2(A)  will  be  discussed  as  it  affects  the 
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production  of  0(^8)  by  the  quenching  reaction.  Some  comparisons  of  the 
photochemical  model  to  5577A  dayglow  observations  will  be  presented. 

The  impact  of  photoelectrons  on  atomic  oxygen  is  the  primary  source 
of  5577 A  dayglow  emission  in  the  thermosphere.  The  reaction  produces  the 
majority  of  the  OC^S)  even  at  the  peak  in  the  emission  layer.  As  can  be  seen 
from  Figure  5.4,  the  observed  profiles  of  5577 A  emission  is  best  reproduced 
by  the  impact  cross  sections  of  Henry  et  al.  (1969)  with  a  multiplying  factor 
of  1.0.  The  magnitude  of  the  impact  cross  sections  determined  by  this  study 
are  also  in  good  agreement  with  the  latest  theoretical  impact  cross  sections 
determined  by  Herrington  and  Burke  (1991,  Personal  Communication)  and 
are  about  a  factor  of  two  smaller  than  the  laboratory  cross  sections  of  Shyn 
et  al.  (1986). 

The  quenching  of  N2(A)  by  atomic  oxygen  is  an  important  source  of 
the  5577 A  emission  at  the  peak  in  the  layer.  While  the  photoelectron  impact 
is  still  the  largest  contributor,  this  reaction  is  needed  to  adequately  fill  in 
the  peak  in  the  emission  layer.  This  reaction  was  studied  for  changes  in 
the  contribution  to  the  volume  emission  by  fixing  the  total  rate  coefficient  to 
the  value  determined  in  Chapter  4  for  the  N2(A),  u=0  level.  The  value  for 
the  total  rate  coefficient  used  in  the  study  is  2.8xlO'H  cm^  s'l  which  falls 
between  2.1±0.4xl0'H  cm^  s*l  from  Piper  (1982)  and  3.5±0.6xl0’^l  cm^  s'l 
from  Thomas  and  Kaufman  (1985).  The  effective  branching  ratio  for  this 
reaction  determined  by  this  study  of  5577A  dayglow  observations  is  0.75±0.25 
for  the  assumed  total  reaction  rate  for  quenching.  The  effective  branching 
ratio  for  the  quenching  of  N2(A),  u=0  determined  by  these  simulations  tend 
to  support  the  results  from  Piper  (1982)  and  De  Souza  et.  al  (1985)  of  0.75±13. 
When  combined  with  the  total  rate  coefficient,  the  branching  ratio  gives  a 
production  rate  for  0(^8)  of  2.01x10*1^  cm^  s'^. 
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Figure  5.6  Results  of  the  tests  of  the  sensitivity  of  the  major 
reaction  parameters  determined  here  to  the  values  used 
in  the  minor  sources. 
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Differences  in  the  vibrational  population  may  account  for  the  discrepancy 
between  the  branching  ratios  determined  by  aeronomical  studies  and  by 
laboratory  experiment.  While  the  vibrational  population  in  laboratory 
studies  is  restricted  to  a  single  level,  the  dayglow  will  have  the  population 
distributed  over  all  vibrational  levels.  This  distribution  plays  an  important 
part  in  determining  the  effective  branching  ratio  for  the  quenching 
reaction.  Since  the  rate  coefficient  has  been  shown  to  change  with 
vibrational  level,  the  contribution  of  each  level  will  depend  on  the  fraction  of 
the  population  it  contains. 

The  branching  ratio  determined  by  Piper  (1982)  was  for  the  entire 
population  of  the  N2(A)  in  either  the  u  =  0  or  the  "0=1  level.  The  effective 
branching  ratio  from  this  study  for  the  OC^S)  dayglow  agrees  well  with  the 
values  of  Piper  (1982)  but  assumes  the  population  is  in  the  n=0  vibrational 
level  only.  In  the  dayglow,  the  vibrational  population  of  N2(A)  is  not  limited 
to  t’  '  lowest  two  levels  as  was  shown  in  Figure  5.5.  By  allowing  for  the 
distribution  of  the  vibrational  population  of  the  N2(A)  state  and  knowing  the 
dependance  of  the  branching  ratio  and  rate  coefficient  by  vibrational  level, 
the  effective  branching  ratio  discrepancy  can  be  resolved.  If  the  N2(A) 
vibrational  levels  responsible  for  the  production  of  the  0(^S)  in  the  dayglow 
are  mainly  the  lower  levels,  then  the  auroral  studies  of  the  effective 
branching  ratios  may  have  better  agreement  with  those  from  laboratory 
studies. 

Using  the  magnitude  for  the  impact  cross  section  and  the  effective 
rate  coefficient  for  the  N2(A)  quenching  reaction  determined  above,  the 
calculated  5577A  volume  emission  rate  from  the  model  can  be  compared  to 
the  inverted  emission  rate  profile  on  an  orbit  by  orbit  basis.  Examples  of 
these  comparisons  are  shown  in  Figures  5.7  and  5.8.  The  figures  show  the 
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relative  importance  of  the  source  reactions  for  the  production  of  OC^S)  for 
the  selected  conditions.  Figure  5.7  shows  AE-D  0846  downleg  which  starts 
at  3°  S  latitude  and  ends  at  about  13°  N  and  the  solar  zenith  angle  changes 
from  about  56°  to  about  65°  along  the  satellite  path.  Figure  5.7a  shows  the 
individual  reactions  used  in  the  model  and  the  total  emission  rate  from  the 
model  compared  to  the  inverted  emission  rate  along  the  satellite  orbit  while 
5.7b  compares  the  error  of  the  inverted  emission  rates  to  the  errors 
attributed  to  the  model.  The  area  between  the  two  lines  is  the  total  of  the 
model  ±20%.  Figure  5.8  is  for  AE-D  0912  downleg  which  travels  from  15°  S 
latitude  to  about  2°  N  while  the  solar  zenith  angle  changes  from  about  57°  to 
about  64°.  Figures  5.8a  and  5.8b  are  similar  to  5.7a  and  5.7b  respectively. 
Electron  impact  is  the  dominant  source  at  all  altitudes  between  250  km  and 
satellite  perigee  with  dissociative  recombination  a  secondary  source  at 
higher  altitudes  and  quenching  of  N2(A)  becoming  more  important  at  lower 
altitudes.  Near  the  peak  of  the  layer,  quenching  of  N2(A)  by  atomic  oxygen 
becomes  an  important  source  of  OC^S)  in  the  5577A  dayglow.  Further  down, 
below  satellite  perigee,  photodissociation  of  O2  (Reaction  [3.9])  quickly 

becomes  the  major  source  of  5577A  emission.  Though  Reactions  [3.3]  and 
[3.12]  become  more  important  in  the  twilight  and  nightglow,  they  are  not 
important  in  the  dayglow  calculations  in  the  photochemical  model 
developed  here. 

Even  though  the  model  is  run  in  equal  steps  in  time  and  the  inverted 
brightnesses  are  in  equal  steps  in  altitude,  good  agreement  is  found 
between  the  observations  and  the  model  predictions,  particularly  at  high 
altitudes.  At  low  altitudes,  where  the  satellite  is  near  perigee,  the  model 
Hoes  not  predict  the  emission  layer  as  well  especially  in  high  solar  zenith 
angle  cases.  This  is  because  the  assumption  of  small  horizontal  gradients 
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Figure  6.7  Results  of  the  model  compared  to  the  inverted  profile 
for  AE-D  0846D.  (a)  the  individual  reactions  and  total 
emission  rate  from  the  model  and  the  inverted  data,  (b) 
Shows  a  comparison  of  the  error  of  the  inverted  data  with 
the  model  ±20%  (area  between  lines) . 
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Figure  6.8  Results  of  the  model  compared  to  the  inverted  profile 
for  AE-D  0912D.  (a)  the  individual  reactions  and  total 
emission  rate  from  the  model  and  the  inverted  data,  (b) 
Shows  a  comparison  of  the  error  of  the  inverted  data  with 
the  model  ±20%  (area  between  lines) . 
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in  brightness  is  not  as  valid  near  perigee.  As  the  satellite  approaches 
perigee,  the  photometer  zenith  angle  will  undergo  a  larger  change  in  a  4 
km  step  than  when  at  400  km.  The  larger  photometer  zenith  angle  change 
means  that  horizontal  gradients  in  the  emission  will  not  be  completely 
corrected  for  in  the  inversion  technique.  Integrating  the  model  volume 
emission  rates  and  comparing  to  the  measured  brightness  will  not  avoid 
the  problem  of  horizontal  gradients.  Since  the  model  only  predicts  the 
emission  at  a  point  from  the  in-situ  measurements,  which  are  along  the 
satellite  track,  the  model  volume  emission  rates  are  not  vertical  profiles. 
Two  consecutive  measurements  of  brightness  are  treated  as  measurements 
from  a  vertical  profile  for  the  differentiation  to  determine  the  volume 
emission  rate.  Integrating  vertically  will  not  predict  the  brightness  seen  by 
the  photometer  looking  radially  out  at  any  given  point. 

5.4  5377A  Emission  in  the  Mid-latitude  Twilight 

At  midlatitudes,  the  upper  atmosphere  undergoes  a  rapid,  drastic 
period  of  change  twice  daily  known  as  twilight.  During  this  period,  the 
atmosphere  is  either  going  from  a  period  of  prolonged  sunlight  (day)  to  a 
period  of  dark  (night)  or  vice  versa.  This  transitional  state  makes  the  study 
of  the  twilight  both  interesting  and  difficult.  This  period  has  been  reviewed 
by  Bates  (1960),  Chamberlain  (1961),  Gadsen  (1967),  and  Hunten  (1967). 

Appendix  D  contains  a  discussion  of  the  twilight  in  general  and  the 
557 7 A  emission  in  the  twilight.  The  modeling  of  the  twilight  is  very  difficult 
and  the  inversion  of  the  observations  in  the  twilight  is  not  as  straight 
forward  as  in  the  dayglow.  The  photochemical  model  developed  for  the 
dayglow  is  able  to  transition  to  the  twilight  conditions.  As  the  solar  zenith 
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to  transition  to  the  twilight  conditions.  As  the  solar  zenith  angle  increases, 
the  various  reactions  that  are  driven  by  the  solar  flux  begin  to  decrease  in 
importance  while  those  driven  by  ion  chemistry  become  the  dominant 
sources.  This  is  illustrated  in  Figure  5.9  which  shows  the  shift  in  reaction 
importance  as  the  orbit  proceeds  into  twilight.  As  the  solar  zenith  angle 
increases  to  90°,  the  contribution  of  the  quenching  reaction  for  N2(A)  drops 
rapidly  as  does  the  electron  impact  on  atomic  oxygen.  The  two  reactions 
involving  02'*'  also  decrease  in  magnitude  but  at  a  much  smaller  rate  than 
the  solar  driven  reactions.  The  overall  results  of  the  change  are  shown  in 
Figure  5.10  which  shows  the  change  in  the  total  volume  emission  rate  as 
the  solar  zenith  angle  increases.  The  changes  illustrated  in  the  figure 
occur  in  less  than  six  minutes  of  time  along  the  satellite  path.  This  rapid 
change  in  chemistry  makes  inversion  and  modelling  of  the  twilight  very 
difficult. 

5.4.1  Model  Results  for  Twilight  Conditions 

The  photochemical  model  developed  and  used  in  the  dayglow  as 
described  above  should  be  able  to  also  model  the  twilight  5577A  emission. 
The  model  may  be  used  in  a  manner  similar  to  the  dayglow  cases. 
Elliptical  orbits  are  used  to  collect  the  in-situ  data  and  measurements  of  the 
brightness  above  the  satellite.  The  model  volume  emission  rates  are  then 
integrated  vertically  above  the  satellite  and  the  results  may  be  compared  to 
measured  brightnesses.  Also,  the  model  volume  emission  rates  will  be 
compared  to  the  results  of  a  circular  orbit  inverted  using  the  sequential 
estimation  method  as  described  in  Chapter  2  and  applied  to  Atmosphere 
Explorer  observations  (Abreu  et  al.,  1991).  Example  results  using  these  two 
methods  will  be  presented  and  briefly  discussed. 
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Figure  5.9  Model  source  reactions  for  5577A  emission  during 
twilight.  For  (a)  70°  (b)  80°  (c)  90°  solar  zenith  angle. 
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Figure  5.10  Comparison  of  total  emission  for  the  model  as 
satellite  proceeds  through  the  terminatior. 
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The  use  of  the  5577A  dayglow  model  in  twilight  is  complicated  by  the 
high  solar  zenith  angles  which  affects  the  calculations  of  the  various 
parameters  in  the  subsiderary  models.  The  complication  due  to  the 
assumption  of  the  change  in  background  emission  above  a  given  level 
thatwas  present  at  perigee  affects  the  twilight  also.  This  restriction  at 
perigee  is  due  to  the  rapid  change  in  photometer  zenith  angle.  As  the 
satellite  traverses  the  twilight  it  is  not  the  change  in  photometer  zenith 
angle  that  changes  so  quickly,  but  the  rapid  change  in  the  emitting  layer 
and  background.  The  use  of  this  method  of  inversion  is  not  capable  of 
handling  the  twilight  emission. 

An  example  of  the  problem  is  shown  in  Figure  5.11.  In  this  figure 
the  model  has  error  bars  of  20%.  While  those  of  the  measured  brightness  is 
set  at  20%  also,  it  could  be  higher  due  to  the  changes  discussed  above.  The 
general  features  of  the  observed  emission  is  well  represented  by  the  model. 
The  integrated  brightness  from  the  model  is  smaller  than  the  observed  but 
the  background  is  not  corrected  for  above  the  altitude  of  the  model. 

The  sequential  estimation  technique  is  better  suited  to  the  changes  in 
the  twilight.  Figure  5.12  shows  the  results  of  the  inversion  without  the 
subtraction  of  Rayleigh  scattered  light.  The  inverted  emission  is  shown  in 
Figure  5.13  after  the  subtraction  of  the  Rayleigh  scattering  and  Figure  5.14 
is  the  model  results.  As  can  be  seen,  the  model  can  reproduce  the  general 
features  of  the  transition  in  the  twilight. 

The  examples  and  results  presented  here  represent  preliminary 
results  of  the  investigation  into  the  twilight  emission  using  the  5577A 
dayglow  model.  These  figures  show  the  5577A  model  as  developed  here  is 
able  to  reasonably  reproduce  the  general  features  of  the  twilight  emission. 
However,  the  rapidly  changing  chemical  processes  in  the  twilight  are 
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Figure  5.11  Twilight  period  for  AE-E  1012D  from  integrated 
model  and  observed  brightness. 
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Figure  5.12  AE-C  Orbit  6646  inverted  5577A  emission  using 
sequential  estimation  without  subtraction  for  Rayleigh 
scattering. 
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Figure  5.13  AE-C  Orbit  6646  inverted  5577A  emission  using 
sequential  estimation  with  subtraction  for  Rayleigh 
scattering. 
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Figure  5.14  Model  5577A  emission  for  orbit  AE-C  6646. 
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difficult  to  model  and  to  observe.  Work  is  still  needed  to  study  the  twilight, 
both  in  the  model  for  the  rapidly  changing  twilight  conditions  and  in  the 
inversion  methods  to  be  used  in  the  twilight.  This  research  is  to  be 
continued  at  a  later  date  and  may  be  extended  to  include  other  emissions  as 
discussed  in  Appendix  D. 


CHAPTER  VI 


CONCLUSIONS 

The  01  emission  at  5577 A  has  a  history  that  goes  back  over  100  years 
since  its  discovery  in  the  aurora  and  nightglow.  This  emission  was  the 
central  part  of  the  study  presented  here.  As  was  pointed  out  in  Chapter  1, 
there  are  still  questions  about  the  various  processes  that  produce  the  0(^S) 
which  gives  rise  to  the  emission  even  though  it  has  been  studied  a  long 
time.  The  work  in  this  thesis  centered  on  answering  some  of  these 
questions  by  using  observations  of  the  emission  taken  by  the  Atmosphere 
Explorer  satellites. 

To  study  the  5577A  dayglow  emission,  a  photochemical  model  was 
developed  for  the  production  of  0(^8)  and  compared  with  the  observations. 
The  description  of  the  satellite  instruments  and  the  analysis  of  the  data  was 
given  in  Chapter  2  and  the  development  of  the  photochemical  model  was 
discussed  in  Chapter  3.  The  rest  of  the  thesis  describes  the  two  most 
important  reactions  to  the  dayglow  production  of  0(^8).  The  results  of  the 
study  are  given  in  the  next  section.  The  continuation  of  the  work  on  this 
project  is  discussed  in  the  last  section.  Also,  work  that  may  be  continued 
from  parts  of  this  project  that  were  only  briefly  started  here  will  be 
discussed. 
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6.1  Results 

The  photochemical  model  developed  here  is  able  to  reproduce  the 

o 

observed  daytime  5577A  thermospheric  emission.  Five  source  reactions  are 
need  in  the  modeling  of  the  emission  profile.  The  reactions  needed  for  the 
model  are  photoelectron  impact  on  atomic  oxygen,  dissociative 
recombination  of  02'*',  quenching  of  N2(A3Zy+)  by  atomic  oxygen,  ion- 
exchange  reaction  between  02‘*'  and  atomic  nitrogen,  and  photo-dissociation 
of  O2  which  becomes  more  important  below  the  peak.  The  ion-exchange 
reaction  between  02"^  and  atomic  nitrogen  is  not  a  major  reaction  needed  to 
explain  the  5577A  emission  profile  but  does  become  important  in  the 
twilight  period. 

The  impact  of  photoelectrons  on  atomic  oxygen  is  the  primary  source 
of  5577A  dayglow  emission  in  the  thermosphere.  The  reaction  produces  the 
majority  of  the  0(^8)  even  at  the  peak  in  the  emission  layer.  As  can  be  seen 
from  Figure  5.4,  the  observed  profiles  of  5577A  emission  is  best  reproduced 
by  the  impact  cross  sections  of  Henry  et  al.  (1969).  The  magnitude  of  the 
impact  cross  sections  determined  by  this  study  are  also  in  good  agreement 
with  the  latest  cross  sections  determined  by  Herrington  and  Burke  (1991, 
Personal  Communication)  and  about  a  factor  of  two  smaller  than  the  cross 
sections  of  Shyn  et  al.  (1986)  [  See  Figure  1.3]. 

The  quenching  of  N2(A)  by  atomic  oxygen  is  an  important  source  of 

o 

the  5577A  emission  at  the  peak  in  the  layer.  While  the  photoelectron  impact 
is  still  the  largest  contributor,  this  reaction  is  needed  to  adequately  fill  in 
the  peak  in  the  emission  layer.  This  reaction  was  studied  for  changes  in 
the  contribution  to  the  volume  emission  from  varying  the  branching  ratio 
for  0(^8)  production.  The  value  of  2.75  x  10'^^  cm3  g-l  foj.  the  total  rate 
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coefficient  determined  in  Chapter  4  agreed  well  with  the  values  of  Piper 
(1982)  and  Thomas  and  Kaufman  (1985)  was  used  in  this  study.  The 

o 

effective  branching  ratio  for  this  reaction  determined  by  this  study  of  5577A 
dayglow  observations  varies  from  approximately  0.65  to  about  0.75.  The 
effective  branching  ratio  for  the  quenching  of  N2(A)  determined  by  these 
simulations  tend  to  support  the  results  from  Piper  (1982)  and  De  Souza  et.  al 
(1985)  of  0.75±13  and  do  not  agree  with  those  of  earlier  studies  in  auroral 
cases. 

The  importance  of  the  vibrational  population  distribution  was 
demonstrated.  The  differences  in  the  vibrational  population  may  account 
for  the  discrepancy  between  the  branching  ratios  determined  by 
aeronomical  studies  and  by  laboratory  experiment.  While  the  vibrational 
population  in  laboratory  studies  is  restricted  to  a  single  level,  the  dayglow 
will  have  the  population  distributed  over  all  vibrational  levels.  This 
distribution  plays  an  important  part  in  determining  the  effective  branching 
ratio  for  the  quenching  reaction.  Since  the  rate  coefficient  has  been  shown 
to  change  with  vibrational  level,  the  contribution  of  each  level  will  depend 
on  the  fraction  of  the  population  it  contains. 

The  branching  ratio  determined  by  Piper  (1982)  was  for  the  entire 
population  of  the  N2(A)  in  either  the  u  =  0  or  the  u  =  1  level.  The  effective 
branching  ratio  from  this  study  for  the  0(^S)  dayglow  agrees  well  with  the 
values  of  Piper  (1982).  In  the  dayglow,  the  vibrational  population  of  N2(A)  is 
not  limited  to  the  lowest  two  levels  as  was  shown  in  Figure  5.5.  By  allowing 
for  the  distribution  of  the  vibrational  population  of  the  N2(A)  state  and  using 
the  the  branching  ratio  and  rate  coefficient  for  the  first  two  vibrational 
levels  from  Piper  (1982),  it  was  shown  that  the  higher  vibrational  levels  are 
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needed  for  about  50%  of  the  production  of  0(lS).  This  points  to  the  need  to 
extend  the  study  of  the  N2(A)  state  to  higher  vibrational  levels. 

The  production  of  0(^8)  in  the  dayglow  from  these  two  reactions 
determined  here  are  not  greatly  affected  by  the  other  reactions  in  the 
photochemical  model.  The  effects  of  changing  the  reaction  rate  of  Reaction 
[3.12]  and  the  branching  ratio  of  dissociative  recombination  of  O2  (Reaction 
[3.3])  were  investigated  in  Section  5.3.2.  This  study  showed  that  the 
variations  of  these  reactions  over  those  given  in  the  literature  will  not 
significantly  affect  the  results  for  the  major  reactions.  The  impact  cross 
section  does  not  shift  a  significant  amount  toward  those  determined  in  the 
laboratory.  Nor  does  the  branching  ratio  for  the  energy  transfer  reaction 
agree  with  the  values  determined  in  the  auroral  studies. 

The  possible  errors  due  to  the  determination  of  the  atomic  oxygen 
densities  by  the  satellite  measurements  was  also  investigated  for  the  effects 
on  the  determined  parameters.  The  possible  changes  due  to  errors  in  the 
atomic  oxygen  density  was  tested  by  assuming  an  error  of  ±30%. 
Increasing  the  oxygen  density  is  equivalent  to  measuring  the  density  low, 
which  is  the  most  likely  case  due  to  surface  recombination.  This  would 
make  the  best  fit  impact  cross  sections  even  smaller,  about  0.75  times  the 
theoretical  cross  sections.  This  makes  the  determined  branching  ratio  even 
larger  to  get  agreement  with  the  peak  in  the  emission  layer.  If  the  oxygen 
density  was  measured  high  by  30%,  the  impact  cross  section  and  the 
branching  ratio  will  be  smaller  but  the  change  is  not  sufficient  to  change 
the  overall  results. 

The  results  of  appl3nng  the  5577A  model  to  the  twilight  period 
demonstrated  that  the  model  is  able  to  reproduce  general  features.  The 
processes  in  the  twilight  are  difficult  to  model  and  to  observe  so 
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comparisons  are  not  exact  or  as  clear  as  in  the  dayglow.  The  5577A  model 
developed  here  is  able  to  reasonably  reproduce  the  twilight  emission  but 
needs  refinining.  The  ability  of  the  model  to  follow  the  changing  twilight 
conditions  was  shown  but  also  points  to  the  need  for  improvment  in 
inversion  methods  for  use  in  the  twilight  and  in  the  model. 

The  results  of  this  work  tend  to  support  a  set  of  values  for  the  various 
reactions  that  were  determined  in  several  separate  studies.  It  is  not 

o 

necessarily  the  final  solution  to  the  question  of  5577A  emission  in  the 
atmosphere.  The  efforts  here  provide  a  basis  for  the  starting  point  of  other 
investigations  that  may  give  these  answers.  These  future  investigations 
are  the  topic  of  the  final  section  of  this  chapter. 

6^  Further  Work 

The  questions  about  the  sources  for  the  greenline  emission  in  the  F- 
region  have  not  been  completly  answered.  The  future  work  on  the  5577A 
emission  based  on  this  study  falls  into  two  groups,  the  work  in  the  dayglow 
itself  and  work  to  expand  the  greenline  into  the  twilight  and  auroral 
regions  to  build  a  complete  model.  Other  possible  research  topics  on  related 
subjects  will  also  be  presented. 

In  the  dayglow,  the  major  work  should  concentrate  on  the  two  most 
important  reactions  as  discussed  here.  For  the  impact  of  photoelectrons  on 
atomic  oxygen,  the  discrepancy  between  the  laboratory  and  the  value  from 
this  study  and  theoretical  calculations  needs  to  be  addressed  with  detailed 
laboratory  investigations.  Improved  experimental  techniques  may  help  to 
solve  this  puzzle.  The  same  can  be  done  for  the  energy  transfer  reaction. 
The  rate  coefficients  and  branching  ratios  need  to  be  determined  by 
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vibrational  level  past  v=0,l  as  done  by  Piper  et  al.,  (1982).  This  can  best  be 
accomplished  in  the  laboratory  with  new  techniques  that  have  recently  been 
developed. 

Once  these  two  major  reactions  are  more  fully  understood,  a  better 
determination  of  the  contribution  by  the  minor  sources  can  be  studied.  This 
is  especially  true  for  the  ion  atom  exchange  reaction  between  02'''  and 
atomic  nitrogen  which  is  important  in  the  twilight  but  not  well  known.  The 
model  for  the  5577A  can  be  extended  into  the  auroral  regions  with 
modification  to  the  photoelectron  model.  This  can  be  done  as  the  model  is 
currently  used,  but  improved  results  from  experiments  will  make  the 
model  of  the  aurora  more  complete  since  the  same  reactions  as  the  dayglow 
are  thought  to  be  important  in  the  aurora.  When  these  steps  are  completed, 
the  model  for  0(^8)  emission  will  be  a  global  model  and  three  dimensional 
maps  of  the  5577A  emission  can  be  constructed. 

Other  possible  research  based  on  this  study  are  the  continued 
expansion  of  the  photochemical  model  to  other  emissions  which  has 
already  been  started  and  expand  the  use  of  the  vibrational  population  model 
for  molecular  nitrogen.  The  photochemical  model  can  be  converted  to 
calculate  other  emissions  based  on  the  same  profiles  as  the  5577A 
greenline.  The  nitrogen  vibrational  population  model  contains  many 
transitions  that  were  not  used  in  this  study  but  others  have  began  using 
them  for  various  studies.  The  triplet  model  used  here  has  a  companion 
model  for  the  low  lying  singlet  transions  of  molecular  nitrogen,  this  model 
can  be  converted  to  the  dayglow  also.  These  two  models  will  then  be  able  to 
describe  the  complete  set  of  nitrogen  transitions  in  the  atmosphere. 
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APPENDIX  A 


SOLAR  FLUX  MODEL 

The  subroutine  SSFLUX  listed  here  calculates  the  solar  EUV  and 
FUV  flux  in  the  range  1  to  1750  A  for  a  specified  level  of  solar  activity.  The 
calling  routine  supplies  a  scaling  switch,  ISC  ALE;  the  daily  10.7  cm  flux, 
F107;  its  81-day  centered  average,  F107A;  and  the  Hydrogen  Lyman-alpha 
flux,  FLYA.  The  subroutine  returns  the  longwave  and  shortwave 
boundaries,  WAVEl  and  WAVE2,  of  the  wavelength  bins  (A),  and  the  solar 
flux  in  each  bin  SFLUX  (photons  cm‘2  s'^).  The  dimension  of  these  arrays 
must  also  be  supplied  by  the  calling  routine  (currently  59).  One  of  four 
methods  are  used,  depending  on  the  value  of  ISCALE,  described  as  follows. 

If  ISCALE=0  the  flux  is  scaled  using  parameterization  methods 
based  on  F107  and  F107A.  For  ionizing  EUV,  Hinteregger's  contrast  ratio 
method  (Hinterregger  et  al.,  1981)  is  used,  based  on  the  Torr  and  Torr  (1985) 
bin  structure  for  reference  spectrum  SC#21REFW.  The  1026A  (H  LyP)  and 
335A  (FeXVI)  enhancement  ratios  are  calculated  from  Hinteregger's 
formula,  using  the  coefficients  which  reduce  to  the  reference  values  at 
F107=75.4,  F107A=71.5.  The  'best  fit'  coefficients  are  not  used  as  they 
produce  some  negative  values  at  low  solar  activity,  but  remain  in  an 
unused  data  statement  for  reference.  The  rest  of  the  spectrum  is  then 
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calculated  from  these  key  emissions  using  Hinteregger's  method.  Scaling 
factors  are  calculated  from  contrast  ratios  in  the  original  spectrum  data 
file.  For  FUV  in  the  1050-1350A  region,  50A  interval  averaging  is  done 
based  on  SC#21REFW,  and  scaling  factors  also  calculated.  For  Hydrogen 
Lyman  alpha  at  1216A,  which  is  treated  separately  as  an  individual  line, 
the  correlation  relationship  derived  from  SME  data  by  Tobiska  and  Barth 
(1989)  is  used.  For  the  Schumann-Runge  continuum,  the  Torr  et  al.  (1980) 

o 

50A  bin  structure  is  used  but  the  coefficients  have  been  adjusted  to  reflect 
SME  and  LASP  rocket  data  (Rottman,  1981  and  1988;  Mount  and  Rottman, 
1983  and  1985). 

If  ISCALE=1  linear  interpolation  between  high  and  low  activity 
spectra  is  done,  based  on  F107  alone,  and  assuming  the  low  activity 
spectrum  corresponds  to  F107=68  and  the  high  activitiy  spectrum  to 
F  107=243.  The  Hinteregger  SC#21REFW  and  F79050  spectra  as  binned  by 
Torr  &  Torr  (1985)  are  used  for  ionizing  EUV.  For  the  1050-1350A  region, 
the  SC#21REFW  spectrum  averaged  into  50A  intervals  is  used  for  low  solar 
activity;  the  high  activity  spectrum  was  obtained  by  scaling  this  spectrum 
using  the  contrast  ratios.  For  the  Hydrogen  Lyman  alpha  and  the 
Schumann-Runge  continuum,  linear  interpolation  amounts  to  the  same 
thing  as  the  aformentioned  parameterization. 

If  ISCALE=2,  the  flux  is  scaled  according  to  the  Tobiska  and  Barth 
method  (Tobiska  and  Barth,  1990).  The  procedure  involves  using  10.7  cm 
flux,  F107,  and  the  Hydrogen  Lyman  alpha  flux,  FLYA,  considered  as 
coronal/transition  region  and  chromospheric  emissions,  respectively. 
Chromospheric  fluxes  are  calculated  from  a  linear  correlation  relationship 
with  FLYA.  Coronal  fluxes  are  calculated  from  a  combination  of  F 107  and 
a  "modeled"  F107  based  on  FLYA.  When  no  value  for  FLYA  is  supplied,  it 
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is  calculated  from  F107  using  the  SME-derived  formula  as  above,  and  then 
used  to  scale  the  other  fluxes.  The  coronal  and  chromospheric  scaled 
fluxes  are  then  summed  to  obtain  the  solar  flux.  The  subroutine  takes  a 
value  of  FLYA  in  photons  cm'2  s'^,  but  converts  it  to  ergs  cm‘2  s'^  for 
internal  computations.  Since  the  wavelength  bins  used  here  extend  past 
those  used  by  Tobiska  &  Barth  for  both  long  and  short  wavelengths,  the 
linear  interpolation  method  is  used  for  those  bins  not  covered  by  the  Tobiska 
&  Barth  model. 

If  ISCALE=3,  the  Woods  and  Rottman  rocket  spectrum  for  10 
November,  1988  is  employed  from  300-1050A.  The  entire  spectmm  of  the 
solar  flux  is  calculated  by  the  contrast  ratio  method  of  Hinteregger 
discussed  above.  The  bins  of  the  woods  and  Rottman  spectrun  are  then 
inserted  in  the  calculated  spectrum. 

The  EUV  fluxes  between  50  and  250A  are  normalized  upwards  (after 
Richards  and  Torr,  1984)  for  both  the  linear  interpolation  and  the  contrast 
ratio  method.  The  normalization  coefficient  is  2  at  F107=68  and  reduces 
linearly  to  1  at  F107=243.  None  of  the  above  models  extends  shortwards  of 

o 

18A,  so  from  1-18  A  an  amalgam  of  sources  are  used  to  derive  an  estimated 
flux,  e.g.,  DeJager,  1964;  Smith  &  Gottlieb,  1974;  Manson,  1977;  Kreplin  et 
a/.,  1977;  Horan  &  Kreplin,  1981;  Wagner,  1988. 
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C  Definitions: 

C  ISCALE  =0  for  Hinteregger  contrast  ratio  method, 

C  =1  for  Hinteregger  linear  interpolation, 

C  =2  for  Tobiska  and  Barth  method 
C  F107  daily  10.7  cm  flux  (l.E-22  W  m-2  Hz-1) 

C  F107A  81-day  centered  average  10.7  cm  flux  (") 

C  FLYA  H  Lyman-alpha  flux  (photons  cm-2  s-1) 

C  WAVEl  longwave  bound  of  spectral  intervals  (Angstroms) 

C  WAVE2  shortwave  bound  of  intervals  (=  WAVEl  for  indiv.  lines) 

C  SFLUX  scaled  solar  flux  returned  by  subroutine  (photons  cm-2  s-1) 

C  LMAX  dimension  of  above  arrays,  currently  =  59 
C  WAVEL  =  WAVEl 
C  WAVES  =  WAVE2 
C  RFLUX  low  solar  activity  reference  flux 
C  XFLUX  high  solar  activity  flux 

C  SCALEl  scaling  factors  for  H  LyB-keyed  chromospheric  emissions 
C  SCALE2  scaling  factors  for  FeXVI-keyed  coronal  emissions 
C  LM  dimension  of  above  arrays,  currently  =  59 

C  SRA  'A  value  for  S-R  continuum  scaling  formula 

C  SRB  'B'  value  for  S-R  continuum  scaling  formula 

CBl  fit  coefficients  for  H  LyB 
C  B2  fit  coefficients  for  FeX^/I 
C  R1  enhancement  ratio  for  H  LyB 
C  R2  enhancement  ratio  for  FeXVI 

C  SFNORM  normalization  factor  for  scaling  flux  shortwards  of  2504 
C  CHINT  intercept  from  chromospheric  Ly-a  iinear  correlation 
C  CHSLP  slope  from  chromospheric  Ly-a  linear  correlation 
C  COINT  intercept  from  coronal  F10.7  linear  correlation 
C  COSLP  slope  from  coronal  F10.7  linear  correlation 
C  W1  weighting  factor  for  coronal  emissions 

C  W2  weighting  factor  for  coronal  emissions 

C  WCHR  weighting  factor  for  chromospheric  emissions 
C  WAVEAV  average  wvl.  for  converting  T&B  flux  from  ergs  to  photons 
C  FLYERG  value  of  Lyman-a  flux  in  erg  cm-2  s-1 
C  F107MD  modeled  F107  value  (uses  Ly-a  in  photons) 

C  FLUXCH  scaled  chromospheric  flux  value  based  on  Ly-a  (in  ergs) 

C  FLUXCO  scaled  coronal  flux  value  based  on  F107 
C 

SUBROUTINE  SSFLUX(ISCALE,F107,F107A,FLYA, WAVEl, WAVE2,SFLUX,LMAX) 
C 

PARAMETER  (LM=59) 

C 

DIMENSION  WAVEl(LMAX),  WAVE2(LMAX),  SFLUX(LMAX), 

+  WAVEL(LM),  WAVES(LM),  RFLUX(LM),  XFLUX(LM), 

+  SCALEl(LM),  SCALE2(LM),  SRA(8),  SRB(8),  Bl(3),  B2(3), 

+  CHINT(LM),  CHSLP(LM),  COINT(LM),  COSLP(LM),  Wl(LM), 

+  W2(LM),  WCHR(LM),  FLUXCH(LM),  FLUXCO(LM),  WAVEAV(LM), 

+  WARKLM) 

new  B's: 


O  Q  O  O 


199 


DATA  Bl/l.O,  0.0138,  0.005/,  B2/1.0,  0.59425,  0.3811/ 
old  B's,  commented  out: 

DATA  Bl/1.31,  0.01106,  0.00492/,  B2/-6.618,  0.66159,  0.38319/ 

DATA  WAVEIV  1750.00,  1700.00,  1650.00,  1600.00,  1550.00,  1500.00, 
+  1450.00, 1400.00, 1350.00, 1300.00, 1250.00, 1215.67, 

+  1200.00, 1150.00, 1100.00, 1050.00, 1031.91, 1025.72, 

+  1000.00,  977.02,  950.00,  900.00,  850.00,  800.00, 

+  789.36,  770.41,  765.15,  750.00,  703.31,  700.00, 

+  650.00,  629.73,  609.76,  600.00,  584.^3,  554.37, 

+  550.00,  500.00,  465.22,  450.00,  400.00,  368.07, 

+  350.00,  303.78,  303.31,  300.00,  284.15,  256.30, 

+  250.00,  200.00,  150.00,  100.00,  50.00,  32.00, 

+  23.00,  16.00,  8.00,  4.00,  2.00/ 

DATA  WAVES/ 1700.00, 1650.00, 1600.00, 1550.00, 1500.00, 1450.00, 
+  1400.00, 1350.00, 1300.00, 1250.00, 1200.00, 1215.67, 

+  1 150.00, 1100.00, 1050.00, 1000.00, 1031.91, 1025.72, 

+  950.00,  977.02,  900.00,  850.00,  800.00,  750.00, 

+  789.36,  770.41,  765.15,  700.00,  703.31,  650.00, 

+  600.00,  629.73,  609.76,  550.00,  584.33,  554.37, 

+  500.00,  450.00,  465.22,  400.00,  350.00,  368.07, 

+  300.00,  303.78,  303.31,  250.00,  284.15,  256.30, 

+  200.00,  150.00,  100.00,  50.00,  32.00,  23.00, 

+  16.00,  8.00,  4.00,  2.00,  1.0(V 

DATAWARl/  0.00,  0.00,  0.00,  0.00,  0.00,  0.00, 

+  0.00,  0.00,  0.00,  0.00,  0.00,  0.00, 

0.00,  0.00,  3.81,  5.00,  5.34,  6.90, 

+  2.44,  7.33,  8.62,  9.60,  4.54,  1.58, 

+  1.28,  0.53,  0.38,  0.34,  0.61,  0.55, 

+  0.38,  2.37,  1.07,  0.41,  1.74,  0.61, 

+  1.67,  0.95,  0.23,  0.58,  1.14,  0.86, 

+  3.32,  7.24,  1.31,  0.00,  0.00,  0.00, 

+  0.00,  0.00,  0.00,  0.00,  0.00,  0.00, 

+  0.00,  0.00,  0.00,  0.00,  0.00' 

DATA  RFLUX/  370.45,  203.69,  96.00,  69.71,  50.70,  26.67, 

+  17.21,  8.26,  12.86,  4.10,  5.20,  272.00, 

+  2.78,  0.70,  3.07,  3.64,  3.18,  4.38, 

+  1.78,  5.96,  4.22,  4.43,  1.93,  0.87, 

+  0.79,  0.24,  0.20,  0.17,  0.39,  0.22, 

+  0.17,  1.50,  0.45,  0.48,  1.58,  0.80, 

+  0.51,  0.31,  0.18,  0.39,  0.21,  0.74, 

+  0.87,  6.00,  0.24,  0.84,  0.10,  0.27, 

+  0.92,  1.84,  0.13,  0.38,  0.0215,  0.0067, 

+  l.E-3,  2.E-3,  l.E-5,  5.E-8,  l.E-KV 

DATAXFLUX/  400.00,  226.00,  107.50,  81.20,  58.70,  29.30, 

+  19.40,  9.90,  30.70,  11.20,  12.00,  438.00, 

+  6.50,  1.60,  6.40,  8.66,  9.04,  13.12, 

+  4.42,  13.18,  12.03,  13.29,  5.01,  2.18, 

+  1.59,  0.67,  0.43,  0.43,  0.72,  0.46, 

+  0.48,  3.02,  1.46,  1.02,  4.86,  1.59, 

+  1.57,  1.67,  036,  0.99,  2.20,  1.39, 

+  5.63,  11.28,  2.50,  4.14,  3.16,  0.59, 

+  3.70,  4.to,  0.3^1,  1.15,  0.18,  0.08, 
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+  2.5E-2,  5.E-2,  8.E^,  3.E-5,  5.E-7/ 

DATASCALEl/12573.0,  9220.0,  4820.0,  4694.0,  3227.0,  1090.0, 
+  922.00,  671.00,8009.80,3186.34,3033.78,  47555., 

+  1692.09,  405.95, 1516.20, 2731.70, 3314.57, 4375.00, 

+  1316.91, 3621.91, 3908.56, 4432.54, 1541.21,  531.73, 

+  364.83,  0.00,  116.00,  129.41,  162.48,  94.07, 

+  41.29,  709.50,  0.00,  268.47, 1561.05,  367.64, 

+  290.06,  184.36,  0.00,  86.15,  7.50,  0.00, 

+  0.00,2220.00,  0.00,  61.00,  0.00,  86.95, 

+  206.00,  135.89,  60.35,  157.12,  7.06,  0.75, 


0.00, 

0.00, 

0.00, 

0.00,  0.00/ 

DATASCALE2/  0.00,  0.00,  0.00,  0.00,  0.00, 

+ 

0.00, 

0.00, 

0.00, 

0.00,  0.00,  0.00, 

+ 

0.00, 

0.00, 

0.00, 

0.00,  0.00,  0.00, 

+ 

0.00, 

0.00, 

0.00, 

0.00,  0.00,  0.00, 

+ 

0.00, 

5.34, 

0.00, 

0.00,  0.00,  0.54, 

+ 

3.30, 

0.00, 

12.60, 

0.00,  0.00,  0.00, 

+ 

5.34, 

11.63, 

2.28, 

5.56,  24.93,  8.16, 

+ 

60.69, 

0.00, 

28.20, 

45.90,  40.80,  1.27, 

+ 

35.47, 

42.80, 

1.12, 

6.19,  1.26,  0.69, 

+ 

0.23, 

0.46, 

7.6E-3 

,  2.9E^,  4.8E-6/ 

DATASRA/  0.169,  0.128,  0.066,  0.064,  0.044,  0.015, 

+  0.013,  0.009/ 

DATASRB/  358.5,  195.3,  91.5,  65.6,  48.0,  25.7, 

+  16.3,  7.72/ 

DATA  CHINT/O.OOOOOE+00, 0.OOOOOE+00, 0.OOOOOE+00, 0.OOOOOE+00, 
+  O.OOOOOE+00, 0.OOOOOE+00, 0.OOOOOE+00, 0.OOOOOE+00, 

+  O.OOOOOE+00, 0.OOOOOE+00, 0.OOOOOE+00, 0.OOOOOE+00, 

+  O.OOOOOE+00, 0.OOOOOE+00, 0.OOOOOE+00, 0. 18223E-01, 

+  -0. 14763E-02, 0. 13765E-02, 0. 10044E-01, 0.48687E-01, 

+  0.79848E-02, 0. 16334E-02, 0.99890E-02, 0. 12362E-01, 

+  0.38317E-02, 0.OOOOOE+00, 0.22214E-02, 0.87883E-03, 

+  0.65437E-02, 0.20197E-02, 0.11897E-02, 0.25303E-01, 

+  O.OOOOOE+00, 0.69859 E-02, 0.15090E-02, 0.15679E-01, 

+  0.71813E-02, 0. 18873E-02, 0.OOOOOE+00, 0.44750E-02, 

+  0.43206E-03, 0.OOOOOE+00, 0.OOOOOE+00, 0.24958E+00, 

+  O.OOOOOE+00, 0. 17434E-03, 0.OOOOOE+00, 0. 1 1585E-01, 

+  0.11532E-01, 0.98362E-02, 0.24461E-02, 0.91696E-02, 

+  0.86457E-03, 0.OOOOOE+00, 0.OOOOOE+OO,  O.OOOOOE+00, 

+  O.OOOOOE+00, 0.OOOOOE+00, 0.OOOOOE+00/ 

DATA  CHSLP/O.OOOOOE+00,  O.OOOOOE+00, 0.OOOOOE+00,  O.OOOOOE+00, 
+  O.OOOOOE+00,O.OOOOOE+00,O.OOOOOE+00,O.OOOOOE+O0, 

+  O.OOOOOE+00,O.OOOOOE+00,O.OOOOOE+00,O.OOOOOE+O0, 

+  O.OOOOOE+00, 0.OOOOOE+00, 0.OOoOOE+00, 0.13097E-01, 

+  0.15670E-01, 0.20808E-01, 0.66142E-02, 0.18084E-01, 

+  0.20603E-01,  0.24699E-01, 0.91077E-02, 0.49109E-02, 

+  0.79057E-03,  O.OOOOOE+00, 0.73979E-03, 0.83584E-03, 

+  0. 1 1269E-02, 0.68429E-03, 0.34588E-03, 0.54965E-02, 

+  O.OOOOOE+00, 0.23091E-02, 0. 13033E-01, 0.32354E-02, 

+  0.27878E-02,  0. 17867E-02, 0.OOOOOE+00, 0.10032E-02, 

+  0.96861E-04, 0.OOOOOE+00, 0.OOOOOE+00, 0.35650E-01, 

+  O.OOOOOE+00, 0.96698E-03, 0.OOOOOE+00, 0. 16549E-02, 

+  0.44742E-02, 0.38025E-02, 0.23655E-02, 0.10185E-01, 

+  0.96072E-03,  O.OOOOOE+00, 0.OOOOOE+00, 0.OOOOOE+00, 


201 


+  O.OOOOOE+00,  O.OOOOOE+OO,  O.OOOOOE+00/ 

DATACOINT/O.OOOOOE+00,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 
+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+00,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO,  0.74083E-03, 

+  O.OOOOOE+00, 0. 13995E-02,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  0. 13215E-02,-0.52657E-03,  O.OOOOOE+OO, 

+  0.17231E-03, 0.OOOOOE+00,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  -0.57828E-02,-0.14201E-01, 0.44921E-02, 0.81510E-03, 

+  -0.38050E-01, 0.25868E-01.-0.74031E-01,  O.OOOOOE+OO, 

+  -0.49781E-01,-0.65154E-01,-0.93662E-01,  O.OOOOOE+OO, 

+  -0.59300E-01, 0.12726E-01, 0.32644E-02,-0.64290E-02, 

+  -0.17720E-01,-0.17972E-01,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO/ 

DATA  COSLP/O.OOOOOE+00,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 
+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO,  0.35613E-04, 

+  O.OOOOOE+OO,  0,67265E-04,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  0.7 1233E-05, 0.45645E-04,  O.OOOOOE+OO, 

+  0.20078E-03,  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  0.98801E-04, 0.23753E-03, 0.45868E-04, 0. 12901E-03, 

+  0.64519E-03, 0.20540E03, 0. 17408E-02,  O.OOOOOE+OO, 

+  0.9031  lE-03, 0. 17093E-02, 0. 13948E-02,  O.OOOOOE+OO, 

+  0. 15455E-02, 0.23847E-02, 0.87499E-04, 0.799 16E-03, 

+  0.34005E-03, 0.3 1881E-03,  O.OOOOOE+OO,  O.OOOOOE+OO, 

+  O.OOOOOE+OO,  O.OOOOOE+OO,  O.OOOOOE+OO/ 

DATA  Wl/  0.00000,  0.00000,  0.00000, 0.00000, 0.00000,  0.00000, 

+  0.00000,  0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 

+  0.00000,0.00000,0.00000,0.00000,0.00000,0.00000, 

+  0.00000,0.00000,0.00000,0.00000,0.00000,0.00000, 

+  0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 2.22570, 

+  5.40850,  0.00000, 0.00000, 0.00000, 0.00000,  0.00000, 

+  1.23810, 0.29122, 0.00000, 0.24172, 0.02037, 0.00000, 

+  0.00000,  0.00000, 0.00000, 0. 18000, 0.00000, 0.00000, 

+  0.33750, 1. 12500, 1. 12500, 0.80000, 0.24000, 0.00000, 

+  0.00000, 0.00000, 0.00000, 0.00000, 0.00000/ 

DATA  W2/  0.00000, 0.00000,  0.00000, 0.00000,  0.00000,  0.00000, 

+  0.00000,0.00000,0.00000,0.00000,0.00000,0.00000, 

+  0.00000,0.00000,0.00000,0.00000,0.00000,0.00000, 

+  0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 0.00000, 

+  0.00000, 1.80579, 0.00000, 0.00000, 0.00000, 0.74191, 

+  3.60570, 0.00000, 0.42239, 0.OOCOO,  0.00000, 0.00000, 

+  0.56275, 0.65524, 1.04710, 0.44890, 1.20180, 0.01272, 

+  1.13808, 0.00000, 1.17450, 2.07000, 1.05000, 0.00000, 

+  1.91250, 1.12500, 1.12500, 1.20000, 1.76000, 2.00000, 

+  0.00000, 0.00000, 0.00000, 0.00000, 0.OOOOC/ 

DATA  WCHR/  0.00000,  0.00000,  0.00000,  0.00000,  0.00000,  0.00000, 

+  0.00000,0.00000,0.00000,0.00000,0.00000,0.00000, 

+  0.00000, 0.00000,  0.00000, 0.04970, 2.06310, 1.07370, 
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+  0.98260, 1.00000,  1.51570, 1.54410, 1.91970, 0.00000, 

+  10.9300, 1.20380,  2.81790, 2.84590, 0.85290, 2.47300, 

+  1.20190, 0.64480, 0.52800, 1.65450, 0.52580, 0.58410, 

+  2.25100, 0.72810, 0.80550, 3.45310, 1.10184, 0.02540, 

+  1.72600, 0.58730, 1.00000, 1.00000, 1.00000, 1.00000, 

+  1.00000, 1.00000, 1.00000, 1.00000, 1.00000, 1.00000, 

+  0.00000, 0.00000, 0.00000, 0.00000, 0.0000(y 

DATAWAVEAV/  0.000,  0.000,  0.000,  0.000, 

+  0.000,  0.000,  0.000,  0.000, 

+  0.000,  0.000,  0.000,  0.000, 

+  0.000,  0.000,  0.000,1025.500, 

+  1031.900, 1025.700,  975.500,  977.020, 

+  925.500,  875.500,  825.500,  775.005, 

+  787.710,  770.410,  765.150,  725.500, 

+  703.360,  675.150,  626.930,  629.730, 

+  609.760,  576.985,  584.330,  554.370, 

+  525.000,  476.185,  465.220,  418.920, 

+  377.915,  368.070,  326.580,  303.780, 

+  303.310,  275.300,  284.150,  256.320, 

+  224.600,  174.340,  124.470,  75.255, 

+  39.620,  24.070,  0.000,  0.000, 

+  0.000,  0.000,  0.000/ 

c 

C  Hinteregger  Method: 

C 

IF  (ISCALE  .EQ.  0  .OR.  ISCALE  .GE.  3)  THEN 
R1  =  Bid)  +  B1(2)*(F107A-71.5)  +  Bl(3)*(F107-F107A+3.9) 
R2  =  B2(l)  +  B2(2)*(F107A-71.5)  +  B2(3)*(F107-F107A+3.9) 
DO  100L=1,LMAX 
IF  (L  .LT.  9)  THEN 
SFLUX(L)  =  SRA(L)  *  F107  +  SRB(L) 

ELSE 

IF{L.EQ.  12)  THEN 
SFLUX(L)  =  (F107  +  218.88)/  1.05453 
ELSE 

SFLUX(L)  =  (RFLUX(L)  +  ((R1-1.)*SCALE1(L) 

+  +  (R2-1.)*SCALE2(L))  /  1000.) 

ENDIF 

ENDIF 

100  CONTINUE 
ELSE 
C 

C  Linear  Interpolation  Method; 

C 

FRAT  =  (F107-68.)  /  (243.-68.) 

DO  200  L=1,LMAX 

SFLUX(L)  =  RFLUX(L)  +  (XFLUX(L)-RFLUX(L))  *  FRAT 
200  CONTINUE 
ENDIF 

C 

C  Normalize  flux  from  250  -  16  A  (after  Richards  &  Torr): 

C 

SFNORM  =  2.  -  1.  *  (F107-68.) /  (243.-68.) 

IF  (SFNORM  -LT.  1.0)  SFNORM  =  1.0 
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DO  300L=1,LMAX 

IF  (WAVEL(L)  .LT.  251.  .AND.  WAVES(L)  .GT.  15.) 

+  SFLUX(L)  =  SFLUX(L)  *  SFNORM 
300  CONTINUE 
C 

C  Tobiska  and  Barth  Method: 

C 

IF  (ISCALE.EQ.2)  THEN 
IF  (FLYALT.l.)THEN 
FLYA  =  (F107  +  218.88)  /  1.05453E-9 
F107MD  =  F107 
ELSE 

F107MD  =  (1.05453E-9*FLYA)  -  218.88 
SFLUX(12)  =  FLYA/1.E9 
ENDIF 

FLYERG  =  FLYA*  12398.*  1.6022E- 12/1215.67 
DO  500  J=  16,54 

FLUXCH(J)=CHINT(J)+(CHSLP(J)*WCHR(J)*FLYERG) 
FLUXCO(J)=COINT(J)+COSLP(J)*(W1(J)*F107MD+W2(J)*F107) 
SFLUX(J)  =  FLUXCH(J)  +  FLUXCO(J) 

SFLUX(J)  =  SFLUX(J)*WAVEAV(J)/(12398.*1.6022E-12*1.E9) 

500  CONTINUE 

ENDIF 

C 

C  Woods  and  Rottman  (10  Nov.  1988)  spectrum: 

C 

IF  (ISCALE  .EQ.  3)  THEN 
DO  550  L=  15,45 
SFLUX(L)  =  WARKL) 

550  CONTINUE 
ENDIF 
C 

C  Convert  from  gigaphotons  to  photons: 

C 

DO  600  L=1,LMAX 
WAVEl(L)  =  WAVEL(L) 

WAVE2(L)  =  WAVES(L) 

IF  (SFLUX(L)  .LT.  0.0)  SFLUX(L)  =  0.0 
SFLUX(L)  =  SFLUX(L)  *  1.E9 
600  CONTINUE 
C 

RETURN 

END 
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APPENDIX  B 


INTERNATIONAL  REFERENCE  IONOSPHERE 

The  Internatiional  Reference  Ionosphere  (IRI)  model  was  described 
briefly  in  Chapter  3.  There  the  use  of  the  Booker  approach  of  building 
profiles  was  noted  (Booker,  1977)  along  with  the  use  of  the  Epstein  functions 
to  construct  layers.  This  apendix  will  describe  these  topics  in  further 
details. 

The  IRI  model  uses  the  first  three  Epstein  functions  to  represent 
altitudnal  variations  as  well  as  special  variations  in  latitude  and  diurnal 
features.  These  three  functions  are  defined  as: 

EPS.,(h;HX,SC)  =  ln(l  +  e’‘)  [B.l] 

EPSo(h;HX,SC)=]/(l  +  e-’‘)  [B.2] 

EPS,(h;HX,SC)  =  e’'/(l+e’‘f  [B.3] 


with: 
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and 


X  = 


h-HX 

SC 


EPS„.=-^EPS,  . 
dx 


[B.4] 


[B.5] 


EPS-1  describes  a  transition,  EPSO  a  step,  and  EPSl  a  peak  at  h=HX  with 
width  SC.  These  functions  are  used  in  combinations  to  construct  the 
profiles  of  the  model.  The  step  function  can  be  used  as  a  filter  function 
taking  a  value  of  0  below  a  certain  altitude  and  a  value  of  1  above.  Rawer 
(1987)  explained  how  three  such  filter  functions  can  be  used  to  construct  a 
jiont  analytical  representation  of  the  electron  density  profile.  The  three 
Epstein  functions  are  shown  in  Figure  B.l. 

The  Booker  function  is  a  single  analytical  function  that  is  capable  of 
representing  simultaneously  the  principle  features  of  the  D-,  E-,  F1-,  and 
F2-Regions  of  the  ionosphere.  In  the  Booker  method  (Booker,  1977),  several 
subsections  were  used  to  construct  a  skeleton  profile  where  each  subsection 
is  approximated  by  a  straight  line.  The  derivative  of  this  profile  function 
can  be  representd  by  a  sum  of  Epstein  functions,  stepping  from  subsection 
to  subsection.  Integrating  from  a  height  ho  to  height  h  gives  the  Booker 
function, 

y{h)-y(h„)  =  (h  -h.)Dy,  +  |;(Dy,.,  -Dy,)SC,AEPS  (B.61 

3=1 


where: 


Figure  B.l  Plot  of  the  Epstein  Functions  used  in  the 
International  Reference  Ionosphere  model  (from  Bilitza, 
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AEPS  =  EPS_j(h;HXj,SCj)  -  EPS.,(ho;HXj.SC  J 


[B.7] 


and 

M  =  the  number  of  subsections 

Dy  =  the  constant  gradients  of  y  in  the  subsections 

HX  =  the  subsection  boundaries 

SC  =  the  thickness  of  the  transition  zones  between  the  subsections. 

The  selection  of  SC  is  critical  to  the  sucess  of  the  function 
representing  an  actual  profile.  If  SC  is  small,  the  model  profile  will  closly 
reproduce  the  skeleton  profile.  If  it  is  too  large,  there  can  be  interference 
between  adjacent  layers.  The  best  results  are  achieved  when 


—  <SC<—  , 
20  10 


[B.8] 


where  D  is  the  height  range  of  the  smaller  of  the  two  subsections.  Also,  the 
choice  of  ho  should  be  in  a  region  where  the  skeleton  and  the  actual  profile 
are  close,  ie.,  outside  of  any  highly  variable  transitions.  These  constraint 
make  it  difficult  to  represent  layer  profiles. 

Layers  are  constructed  by  using  the  Epstein  transition  functions  as 
introduced  by  Rawer  (1984).  By  combining  the  Epstein  function  with  a 
linear  term,  the  layer  function  is, 

L(h;h^,HX,SC)  =  AEPS-(h-h,„)EPSo(h„;HX,SC)/SC  [B.9] 


where 


AEPS  =  EPS.,(h;HX,SC)  -  EPS.i(h„,;HX,SC). 


P.IO] 
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The  layer  peaks  at  h=hm-  Figure  B.2  shows  examples  of  layers  for  a  variety 
of  values  of  HX  and  SC. 
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APPENDIX  C 


PHOTOELECTRON  MODEL 


The  model  of  the  photoelectons  used  in  this  work  is  essentially  the 
sam  as  described  in  Solomon  (1987),  Solomon  et  al.  (1988),  and  Solomon  and 
Abreu  (1989).  Here,  a  discription  of  the  electron  transport  equation  is  given 
and  the  impact  cross  sections  for  excitation  and  ionization  are  listed.  The 
model  is  discussed  in  general  terms  in  Chapter  3. 

The  electron  transport  calculations  used  in  the  model  use  the 
continuity  equations  of  Banks  et  al.  (1974)  for  the  flux  along  a  magnetic  field 
line.  Neglecting  gravity,  parallel  electric  fields  and  magnetic  field 
convergence,  thes  transport  equations  are. 


^  =  -7  X  k  {^a  +  Pc  -  Z  ^kPe  ^  . 

as  p  T  '  P  k  2p  p 


[C.l] 


dO' 

ds 


X n. (o;  +  )<!>■  +  ^  ^  ■ 

p  k  p  k  ^p  p 


[C.2] 


where: 


q*(E,s)  =  5^nJs)5^{R<I)-(E',s)  +  [l-R]cD"(E',s)} 

k  l.f>K 


[C.3] 
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q-(E.s)  =  Xn,(s)X{R<I>*(E’.s)  +  [l-R]cD-(E',s)},  [C.4] 

k  j.oE 

where: 

R  =  p;,(E>‘(E'^E)  [C.5] 

and 

s=distance  along  a  magnetic  field  line  (positive  upward) 

O 

<J>*  = — =  electron  fiux  upward  along  s 

O. 

O  =  — -  =electron  flux  downward  along  s 

n^(s)=neutral  species  number  density 
Pg(E)  =elastic  backscatter  probability 
Og(E)  =  elastic  cross  section 
q(E,s)=secondary  electron  production  rate 
q*  =electron  production  rate  due  to  cascading 
Pgj  =inelastic  backscatter  probability 
o^,  =inelastic  cross  section 

J 

j=process  index 

k=neutral  species  index 

p=mean  pitch  angle  cosine  (assumed  0.577). 

The  equations  are  solved  by  implicit  integration  with  boundary  conditions  of 
the  fluxes  are  equal  and  near  zero  at  the  bottom  and  specified  by  an  incident 
electron  flux  at  the  top. 

To  calculate  the  electron  fluxes  at  a  given  point  in  the  atmosphere, 
the  electrons  are  binned  according  to  energy  in  variable  sized  energy  bins. 
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Starting  with  the  highest  energy  electrons,  the  secondary  and  degraded 
primary  electrons  are  computed  and  added  to  the  electron  flux  at  each 
altitude  and  each  bin  of  lower  energy.  The  excitation  and  ionization  due  to 
this  energy  of  electrons  is  also  found.  Then  the  next  highest  energy  bin  is 
calculated  and  the  results  added  to  those  below  it  in  energy.  The  degraded 
fluxes  are  filled  in  the  neighboring  bins  by  a  weighted  distribution  method 
similar  to  that  of  Swartz  (1985)  at  allow  for  the  differences  between  electron 
energy  and  bin  mean  energy. 

The  cross  sections  for  excitation  and  ionization  of  the  neutral  species 
were  calculated  using  the  analytic  fit  parameters  of  Green  and  Stolarski 
(1972)  and  Jackman  et  al.  (1977).  For  the  excitation  parameters  of  O,  N2, 
and  02,  The  Green  and  Stolarski  formula  was  used. 


where  =  6.51  x  lO'l"^  cm"'^  eV'^,  W  is  the  excitation  energ>,  E  is  the  incident 
electron  energy,  and  Aq,  v,Q  are  the  fitting  parameters.  The  parameters 

for  the  excitation  cross  sections  are  given  in  Table  C.l. 

For  ionization  of  atomic  oxygen,  the  Green  and  Stolarski  formula 
was  used  for  the  cross  section  and  the  Opal  et  al.  (1971)  description  of  the 
secondary  electron  production  spectrum  was  employed, 

“■(E.)  =  q„A  I  ^(l-e-)’dw 


[C.7] 
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Table  C.l  Parameters  for  the  Excitation  Cross  Sections  used  by 
the  Photoelectron  model. 
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S{Ep,E,)  = 


itE 


i+(e«/e)  ’ 


[C.8] 


where  I  is  the  ionization  potential,  Ep  is  the  energy  of  the  primary  electron, 
Es  the  energy  of  the  secondary,  e=EpAV,  Eis  an  empirically  derived  energy 
degradation  coefficient  (15.0  eV  for  atomic  oxygen),  and  A,  Q,  7,  v  are  the 
fitting  parameters.  The  atomic  oxygen  ionization  parameters  are  given  in 
Table  C.2. 

For  the  molecules,  the  Jackman  et  al.  formula  for  the  differential 
cross  sections  were  used. 


S,(E,.E.) 


A(E,)r^(Ep) 
(E.-T„(E,))%r“(Ep)  ’ 


[C.9] 


where: 


A(E)  = 


^  K  ^ 

v% 


In 


I  ^ 


[C.IO] 


EsEp 

Ep  +  r, 


[C.ll] 


To(eJ  =  t, 


Ep-Tb  ■ 


[C.121 
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These  parameters  are  given  in  Table  C.3.  Figure  C.l  shows  the  electron 
impact  elastic  cross  sections  and  backscatter  ratios  used  in  the  two  stream 
model  (After  Solomon,  1987). 


217 


iiaiiSii 

11 

II 

9S 

ii 

«-i  f-i  O)  CO  o 

^  S3  S  S  S  s  g3 

q 

^  » 

18.75 

8 

23.60 

Ills 

lA 

oo' 

CO  a 
CO  c 

F 

0  CO 

o'  eo 

H  F-t 

CO 

CO 

fH 

13.8 

13.8 

«  N  N  00  ^  <o  o 

^  s;  s?  S5 

o 

t- 

CO  (i 
fH  ^ 

fJ  e 
CO  c 

i 

9 

47.20 

80.00 

■■ 

OOOOl 

C 

S' 

tx 

vuwi. 

O'OOOI 

1000.0 

O'OOOI 

1000.0 

1000.0 

T. 

ar  Oxy 
1.86 
1.86 
1.86 
1.86 
1.86 
1.86 

1.86 

S 

2; 

w 

fH  F" 

■1 

: 

4.71 

4.71 

4.71 

TTOUl.l 

B1 

8 

d 

<s 

3 

u 

"o 

0.00 

nnn 

\S 

>  d 

000 

0.00 

0.00 

4'  "“ssssss 

^  o  o  d  d  d  d 

000 

s 

0.00 

non 

00*0 

OO'O 

88 
o  d 

gpmi 

3.76 

1.74 

1  74 

1.74 

rH 

1.74 

1.74 

ipllii 

8 

d 

000 

000 

OO’O 

8 

d 

88 
d  d 

gliili 

A 

d 

1990 

f:: 

CO 

d 

0.371 

0.530 

1 

1 12.10 

1 16.10 

1 16.90 

1 18.20 

1 20.00 

1  23.00 

B 

S« 

A 

1 

1 

1  16.73 

uo 

C"* 

ad 

F-< 

1 

38 

11 

c 

m 


0.001 


0.01 


10 


0.1  1 
Energy  (keV) 


Solid  line— O  Dashed  lioe — Oj  Dotted  line — N: 


Filfure  C.l  Elastic  Cross  sections  and  Backscatter  Ratios  used 
in  the  Photoelectron  model  (from  Solomon,  1987). 
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APPENDIX  D 


INTRODUCTION  TO  TWILIGHT 

Twilight  is  defined  by  Chamberlain  (1961)  as  the  time  when  sunlight 
is  shining  on  an  emitting  layer  of  the  atmosphere  from  below.  This  is 
distinctive  from  the  daytime  when  the  sun  is  shining  on  the  emitting  layer 
from  above  and  at  night  when  the  emitting  layer  receives  no  direct  sunlight 
at  all.  Ground  based  twilight  observations  may  be  used  as  a  substitute  for  a 
series  of  rocket  flights  which  will  give  much  better  continuity  of  coverage 
(Hunten,  1967). 

The  twilight  has  several  differ;.iil  periods  all  of  which  are  based  on 
the  relative  position  of  the  sun  to  the  emitting  region  of  the  atmosphere. 
These  divisions  are  based  on  the  instant  when  the  geocentric  zenith  angles 
of  the  central  point  of  the  solar  disk  in  increments  of  about  6  degrees  as 
indicated  in  Table  D.l.  The  times  of  these  periods  are  tabulated  annually  in 
the  Astronomical  Almanac  (from  US  Government  Printing  Office, 
Washington  and  HM  Stationary  Office,  London). 

The  simplest  figure  for  twilight  is  given  in  Figure  D.l  where  the 
atmosphere  is  singly  scattering.  The  more  realistic  picture  is  given  in 
Figure  D.2  where  complications  arise  from  multiple  scattering,  adsorption, 
and  refraction  by  the  atmosphere  and  by  reflections  from  the 


Figure  D.l  The  simplest  interpretation  of  twilight  observations 
by  a  singly  scattering  Rayleigh  atmosphere. 
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Figure  D.2  Schematic  indications  of  the  interfering  processes 
that  can  occur  in  the  twilight  including  multiple 
scattering. 
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surface  and  clouds.  These  topics  are  reviewed  by  Gadsden  (1967)  and 
discussed  in  detail  by  Rozenberg  (1966). 

Observations  of  the  daytime  airglow  began  about  1961,  study  of  the 
twilight  started  about  1938  (Kvifte,  1969).  The  study  of  the  emissions  during 
twilight  can  provide  valuable  information  about  the  atmospheric  structure, 
dynamics  and  chemistry.  The  period  of  twilight  covers  changes  from  when 
the  daytime  airglow  which  is  dominated  by  solar  and  photoelectron  driven 
processes  into  the  nighttime  where  the  emissions  are  provided  by  various 
chemiluminescent  processes.  During  this  period  the  total  sky  brightness 
diminishes  rapidly  and  the  weaker  radiation  become  easier  to  observe.  The 
change  in  the  atmosphere  can  be  used  to  study  either  the  dayglow  or 
nightglow  by  treating  it  as  perturbation  of  a  quasi-steady  state  (Gadsden, 
1967). 

The  twilight  emissions  have  been  subdivided  by  Hunten  (1967)  into 
groups  that  are  separated  by  altitude.  The  lowest  is  from  O2  bands  at  30  to 
50  km,  with  some  emission  of  OH  in  the  same  region.  The  alkali  metals  of 
Na,  K,  Li  and  Ca'*'  form  the  next  group.  The  first  three  are  in  a  layer  at 
about  90  km  with  Ca'*’  being  somewhat  higher.  The  emissions  from  these 
metals  is  due  to  resonance  scattering.  The  third  group  is  the  F-region 
emissions  from  forbidden  transitions  of  atomic  oxygen  and  nitrogen,  N2‘*' 
and  orthohelium. 

The  excitation  processes  for  the  emission  in  the  twilight  is  expected  to 
be  those  in  the  dayglow  only  at  a  much  smaller  rate.  This  is  especially  true 
for  those  processes  that  are  directly  driven  by  solar  and  photoelectron 
processes.  As  the  solar  zenith  angles  increase,  the  path  of  the  light 
through  the  atmosphere  lengthens.  This  screening  by  the  lower 
atmosphere  to  the  radiation  changes  the  effective  radius  of  the  earth  as 
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shown  in  Figure  D.3.  Due  to  increased  attenuation,  from  adsorption  and 
scattering,  the  flux  will  decrease  and  photoelectron  production  will  drop  for 
a  given  point  in  the  atmosphere. 

The  atomic  oxygen  lines  observed  in  the  twilight  glow  at  the  F-region 
altitudes  are  the  red  (6300A)  and  green  (5577A)  forbidden  lines  for  Ol^D)  and 
0(^8)  and  a  weak  4368A  line  from  the  allowed  transition  of  (43p-33S).  At  F- 
regions,  the  red  line  is  greater  in  intensity  than  the  green  line  emission. 
Both  lines  have  the  same  excitation  mechanisms  as  in  the  dayglow. 
Photodissociation  of  O2  in  the  Schumann-Runge  continuum  produces  the 
bulk  of  the  0(^D).  The  0(iS)  is  mainly  produced  by  the  photoelectron 
excitation  but  as  twilight  progresses,  dissociative  recombination  is  the 
responsible  mechanism  which  will  continue  into  the  nightglow. 

e> 

At  twilight,  the  6300A  01  emission  has  an  enhancement  over  the 
nightglow  (Chamberlain,  1961;  Gadsden,  1967;  Hunten,  1967,  Kvifle,  1969). 
The  5577A  greenline  has  been  reported  to  have  an  enhancement,  much 
weaker  than  that  for  the  red  line,  of  a  factor  of  two  (Dufay  and  Dufay,  1948). 
It  was  not  observed  for  many  years,  until  Megill  (1960)  using  a  birefringent 
filter  substantiated  the  observations.  The  expected  photon  emission  from 
fluorescence  of  2972A  radiation  is  not  expected  to  contribute  significantly  to 
the  twilight  5577A  emission  due  to  the  low  screening  height  for  2972A 
radiation  at  about  50  km  due  to  ozone  absorption  and  the  altitude  of  the  peak 
in  atomic  oxygen  near  100  km  (Gadsden,  1967).  This  excludes  resonance 
fluorescence  as  a  source  of  the  late  twilight  enhancement.  The  most  likely 
source  is  the  dissociative  recombination  of  02"^,  (Reaction  [3.3]),  which 
would  also  provide  the  enhancement  for  the  red  line. 

Frederick  et  al.  (1976)  suggested  that  Reaction  [3.12],  the  ion-atom 
interchange  between  atomic  nitrogen  and  02'^,  can  explain  the  structure  of 
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Zenith 


Figure  D.3  Geometry  of  twilight  observations  from  the  ground 
The  observer  at  O  measures  the  emission  at  P.  The 
geometrical  shadow  height  is  Zg  and  the  true  screening 
height  is  Z]. 
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the  5577A  emission  in  the  twilight.  The  behavior  of  the  twilight  bbllk 
emission  has  been  studied  by  Schaeffer  and  Feldman  (1972),  Hays  and 
Sharp  (1973),  and  Kopp  et  al.  (1977).  Kopp  et  al.  (1977)  demonstrated  that  the 
5577A  emission  is  asymmetric  in  the  morning  and  evening  terminators. 
This  is  shown  in  Figure  D.4  and  they  concluded  the  differences  can  be 
explained  by  using  Reaction  [3.12]  with  a  rate  coefficient  of  2.5  x  10' cm3  §- 
1,  which  supports  the  results  of  Frederick  et  al.  (1976). 
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ALTrruOE  (Km) 


Figfure  D.4  Contour  of  the  5577 A  emission  rate  as  a  function  of 
altitude  and  SZA.  (After  Kopp  et  al,  1977). 
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